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1. INTRODUCTION 
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AnA rh^ ^ result of the major significance of flight velocity on iet noise 
fUalt ^PP^''^";.'|;^°"?'stency between jet noise trends observed in actual 
9 ht and in flight simulation facilities, NASA Lewis Research Laboratory 
decided to embark on a study of the problems of understanding the generation 
propagation, and measurement of exhaust-generated noise from jets fn flight ’ 
The work described in this report represents a contracted stud/of !he funda- 

no?se f%m^a*^^ controlling the generation, propagation, and measurement of 
noise from a conical jet in forward flight. 
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More specifically, the primary objective of this program was to obtain 
a fundamental study of the problems Involved in either theoretically predict- 
ing or empirically extrapolating the jet noise field (from both subsonic and 
supersonic exhausts) received on the ground from an aircraft flying overhead 
in typical takeoff and landing approach patterns. Such prediction or extra- 
polation would be made from static ground tests both with and without 
relative velocity effects on the nozzle configuration. The ultimate goal of 
this program was to make detailed recommendations based on the studies 
conducted under the following four specific tasks, for a possible follow-on 
program. 


Task 1: Literature Survey and Preliminary Investigation 

Task 2: Acoustic Propagation Effects 

Task 3; Source Alteration Effects 

Task k: Investigation of Verification Techniques 

The first task involved a compilation and assessment of the existing 
literature on jet noise flight effects and related subjects. This was com- 
pleted early in the program, and a literature survey report was submitted for 
NASA's internal use. The remainder of the program comprised a basic study of 
the generation and radiation of noise from jets in flight. In addition, the 
problems of measurement of noise from aircraft in flight were studied in some 
detai 1 . 

The intent of this contract was that it act as a program definition 
phase for more comprehensive work to follow. This intent has been met, and 
in several areas, new data and theory are presented which lead to definitive 
conclusions regarding the generation of noise of jets in flight. 

This report is compiled in three major sections which describe the 
technical work accomplished in this contract. A synopsis of the technical 
work described in these three sections is given in Section 5* 

Contributors to the work presented in this report are Robert H. Burrin, 
Christopher L. Morfey, Philip J. Morris, C. Benton Reid, H. K. Tanna, 

Brian J. Tester, and M. Clay Whiffen. 
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2. INFLIGHT SIMULATION EXPERIMENTS ON JET NOISE 


2.1 OBSERVED EFFECTS OF FORWARD MOTION ON JET NOISE 


During the initial period of intensive research on jet noise and the 
development of jet noise suppressors in the 1950's and early 1960's, the 
effect of aircraft motion on the source (the "relative velocity" effect) was 
considered, and the reduction in noise level at peak polar angle was duly 
noted. Attempts were made to describe this reduction by using the relative 
velocity rather than the absolute jet velocity in noise estimation proce- 
dures. It was observed, however, that this did not fully agree with the 
measured results, and a further term was introduced by Greatrex (ref. 5) 
and Coles (ref. 6) in an attempt to allow for an increase in noise genera- 
tion due to elongation of the jet mixing region which was known to take 
place in flight. During this period, some theoretical work was carried out 
by Ffowcs Williams which included the effects of aircraft motion (ref. 7); 
however, this does not appear to have been widely used. 

It is probably correct that, although many flight tests were made and 
compared with static results, most of these were of a developmental nature 
rather than being oriented toward engine research purposes. Hence, only 
peak levels were required, and the more fundamental aspects of comparing the 
directivities and spectral shapes in detail were rarely dealt with. In 
these early investigations, the effects of ground reflection were certainly 
not well considered, and atmospheric attenuation data were ill-defined. 

Also, the standard frequency analysis system was at that time in octave 
bands, which would obscure important differences that may exist between 
static and flight spectra. 

The corrections usually applied to static acoustic test data to theo- 
retically convert the results to the moving source or inflight results 
relative to a stationary observer are (a) a Doppler shift in frequency, (b) 
a change in filter band level due to change in Doppler-shifted center fre- 
quency, and (c) a so-called dynamic effect or convective amplification due 
to source motion relative to the fixed observer. 

Corrections for effects (a) and (b) have been given by Mangiarotty and 
Turner (ref. 8). The source motion effect (often referred to as the dynamic 
effect m the past) was derived by Ffowcs Williams (ref. 7). The effects of 
relative velocity on the source strengths and on radiation efficiency through 
flow/acoustic interactions have not been determined. In the SAE prediction 
method (ref. 9), the noise generated by a jet in motion is simply taken to 
be equal to that of a static jet with exhaust velocity equal to the relative 
velocity of the moving jet (jet velocity minus forward speed of the jet). 

Recent Measurements and Observations 

A considerable amount of noise measurements, both from model-scale 
configurations operated in the inflight simulation mode and from full-scale 
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flight tests, has been published, but the majority of these studies have 
been concentrated on various suppressor nozzle configurations [Coles, et al y 
(ref. 10); Brausch (ref. 11); Burley and Karabinus (ref. 12); Hoch and 
Hawkins (ref. 13); Von Glahn, et aly (ref. l4); Burley and Johns (ref. 15); 
Burley, et aly (ref. l6); Chamberlain (ref. 17); Von Glahn and Goodykoontz 
(ref. l8); Burley and Head (ref. 19); Gubkina and Mel’nikov (ref. 20); 

Von Glahn, et alj (ref. 21); Brooks and Woodrow (ref. 22); Bushel 1 (ref. l); 
Cocking and Bryce (ref. 23)]. Data on the basic conical (or convergent) 
nozzle are, therefore, rather scarce. Furthermore, there is considerable 
disagreement between the results of various investigations, and the situation 
is quite inconclusive at the present time. The use of the corrections given 
above does not explain the differences measured between inflight and static 
tests. Several factors can be suggested which play very important roles in 
producing widely varying and misleading conclusions in such experimental 
programs. In particular, these are internal or extraneous noise sources in 
the rig or the engine under test; geometry of the configuration under test 
(which must influence the flow properties in the shear layer of the Jet); 
Improper account of the effects of environmental factors; and finally, 
inaccurate and/or inadequate techniques for the acquisition and analysis of 
the noise data. 

The discovery of the importance of internal core noise or "tailpipe” 
noise under static conditions on many of the test vehicles gave an explana- 
tion of why the expected relative velocity effects may not be observed in 
flight. This was that the tailpipe noise sources would dominate the in- 
flight noise and there is no obvious reason why noise from these sources 
should fall with forward velocity. As a result, many of the earlier measure- 
ments carried out to study the effects of forward motion on jet noise, from 
both conical and suppression nozzles, have been invalidated. 

However, further investigation has shown that in the jet velocity 
regimes where internal noise would not be thought to dominate even in flight, 
the changes of the measured noise characteristics (in particular, the field 
shape) with jet and aircraft velocity, are not consistent. Even the shock- 
associated noise level at 90° to the jet axis is found to be increased in 
flight, contrary to the theoretical suggestion that dynamic effects or con- 
vective amplification due to aircraft motion should be zero or negligible at 
this angle. Therefore, it is clear that further understanding of the 
characteristics of jet noise is required under flight conditions. 

Inspection of the more recent full-scale flight results [Hoch and 
Hawkins (ref. 13), Brooks and Woodrow (ref. 22), and Bushel 1 (ref. 1)] and 
comparison with static test data reveal that, although a "relative velocity" 
reduction in turbulent mixing noise is observed close to the downstream jet 
axis in flight, no general reduction is observed around the complete field. 
Indeed, the observations range from small reductions in the rear arc to 
actual amplifications of the noise in the forward quadrant; the level in the 
forward arc rises with aircraft speed instead of falling as a simple rela- 
tive velocity effect would predict. In contrast, results from the inflight 
simulation experiments at model scale conducted by NGTE (ref. 23) and NASA- 
Lewis (ref. 21), indicate that noise reductions with forward velocity are 
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present around the complete field, and in general, the magnitude of the 
relative velocity effect increases as the observer moves from the forward 
arc towards the downstream jet axis. This fundamental discrepancy between 
the inflight simulation results and the full-scale flight experiments is 
illustrated qualitatively in Figure 2.1. The changes in OASPL with forward 
motion can also be displayed conveniently as a relative velocity (jet 
velocity minus forward speed) dependence. This is shown in Figure 2.2, 
where the inflight effects are expressed quantitatively in terms of the 
relative velocity exponent, m. 

To summarize, it is clear that there are significant differences be- 
tween the inflight effects as observed from flight tests and those from 
model-scale inflight simulation experiments. These discrepancies need to be 
resolved before the flight simulation techniques can be used as a long-term 
research tool to obtain a fundamental understanding of forward velocity 
effects on jet noise. A significant effort in the present program has been 
directed toward this problem, as will become apparent in the later parts of 
this report. 


2.2 FLIGHT SIMULATION TECHNIQUES 


In the preceding section, published inflight effects on jet noise, 
derived from facilities which simulate forward motion, were briefly intro- 
duced without mentioning any specific simulation concept. The purpose of 
this section is to present brief descriptions of various flight simulation 
methods that are either in use or gaining increasing interest at the present 
time. But before this is done, it is perhaps worthwhile to summarize the 
problems involved with actual flight tests. 

2.2.1 Full-Scale Flight Testing 

Full-scale flight testing is a direct method giving the required noise 
measurements with no intermediate steps or recourse to theory. However, the 
cost associated with any comprehensive flight test program is prohibitive, 
and there are several disadvantages which are difficult to overcome: 

(1) The internal noise sources are almost always significant in a full- 
scale engine, and their contribution to the total sound field must be 
established and, if necessary, adequately suppressed to levels well below 
the unsilenced jet noise so that no significant contamination occurs when 
the magnitudes of turbulent mixing noise decrease in flight. 

(2) The envelope or test window of jet operating conditions (exhaust 
velocity and temperature) and flight speed is usually restricted In a flight 
test program. 

(3) The sound generated by the jet exhaust flow has to propagate 
through a complex inhomogeneous environment on its way to the ground-based 
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observer, and many of the phenomena associated with propagation of sound 
through a real atmosphere have yet to be understood and properly quantified. 

(4) The acquisition of flyover noise data and the correction procedure 
necessary to yield free-field data is extremely complicated. It involves 
statistical averaging of nonstationary acoustic data, accurate determination 
of aircraft position and speed, proper account of ground reflection, and 
accurate atmospheric absorption corrections. These aspects of flyover data 
acquisition are discussed in detail in Appendix 2A of this report. 

In view of the above difficulties associated with flight testing, 
several methods for investigating inflight effects on jet noise have evolved 
in recent years. The four most promising types of flight simulation facili- 
ties are: 

(1) jet mounted on a ground based moving vehicle 

(2) jet mounted at the end of a rotating arm (a spinning rig), 

(3) stationary jet in a large-scale acoust i ca 1 ly-treated wind 
tunnel, and 

(4) stationary jet immersed in a larger surrounding jet (a free-jet 
faci 1 i ty) . 

Each of these has its own merits as well as facility oriented problems, and 
these are discussed briefly below. 


2.2.2 Ground-Based Moving Vehicle 

Investigation of inflight effects on jet noise by mounting the engine 
configuration on a ground-based vehicle, which can be translated at 
constant speed, appears attractive at first sight since the relative motion 
between a moving source and a stationary observer can be simulated exactly. 
Two examples of such a facility are the high speed train and the rocket 
powered sled. However, these facilities are accompanied by problems which 
are similar to those encountered in aircraft flyover noise measurements. 

2.2.3 Spinning Rig 

The use of a spinning rig, where the model -scale jet exhaust is 
mounted at the end of a rigid arm that rotates at constant speed, also pro- 
vides an exact relative motion simulation. The main advantages of this 
method are that the costs involved are much smaller than flyover testing 
and the atmospheric propagation effects are diminished by an order of 
magnitude since the propagation distances are smaller. However, there are 
severe technical problems associated in operating such a facility. In 
particular, because the angular position of the moving jet source relative 
to fixed microphones is changing rapidly, the problems associated with 
noise data acquisition and reduction are very complicated. It should also 
be noted that even if the contamination from internal noise sources and 
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rig-generated self noise can be reduced to acceptable levels* the undesirable 
effects of centrifugal forces in producing distortion of the jet flow have 
yet to be established. 

2. 2.4 Large-Scale Anechoic Wind Tunnel 

The step from static testing in relatively idealized test facilities to 
full-scale flight testing is large. Therefore, a very useful intermediate 
step can be provided by operating the jet configuration under examination in 
a wind tunnel, in which the effect of forward motion is produced by the 
tunnel air flow. 

A comprehensive review of the use of wind tunnels for acbustic and 
aerodynamic measurements in the inflight simulation mode, and current re- 
search aimed at the design and operation of large wind tunnels, has been 
published recently by AGARD (ref. 24, 25 , and 26). 

Aerodynamic measurements on a moving configuration do not appear very 
practical, but can be achieved relatively easily in a wind tunnel environ- 
ment. The acoustic measurements can be achieved by making the test section 
anechoic, although the costs required to modify any large-scale tunnel, such 
as the NASA-Ames 12.2 m x 24.4 m facility, to obtain acceptable background 
noise levels and to minimize wall reflections, are high. In addition, 
operating costs are considerable for comprehensive evaluation of inflight 
effects on jet noise. 

Simulation of aircraft motion in such a manner is attractive in that a 
model rig free from internal noise can be set up and tested in a controlled 
anechoic environment. The jet operating conditions and the tunnel speed 
can be controlled independently. In addition, the measurements are carried 
out In essentially a static (steady state) environment. However, the 
quality of acoustic measurements is highly dependent upon having a quiet 
means of propulsion for the tunnel and in developing microphones which are 
Insensitive to the air flow across them. The noise levels induced by air 
flow over the microphones placed within the tunnel flow are usually high, 
and therefore preclude any meaningful measurements at reasonably high 
tunnel velocities of interest. 

A minor disadvantage of this technique is the lack of a true relative 
velocity effect between source and observer, unless the microphone is 
traversed past the nozzle at the moving stream velocity. However, the 
effects of forward motion on jet structure as well as changes in noise 
characteristics for an observer moving with the jet can be established. The 
translation of acoustic data from a wind tunnel facility to the fly-over 
situation can be achieved far more readily within the bounds of existing 
theoretical expertise. 

2.2.5 Free-Jet Anechoic Wind Tunnel 

This flight simulation method removes many of the limitations of con- 
ventional wind tunnels for acoustic testing. The aircraft forward motion in 
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a free-jet or open-jet facility is simulated by immersing the primary jet 
exhaust configuration in a larger surrounding jet. The microphones are 
placed in the far field outside the free jet flow, and the complete test 
section is enclosed in an anechoic chamber to provide a free-field environ- 
ment. In addition to retaining other advantages of a large-scale wind 
tunnel, a carefully designed free-jet facility is also capable of providing 
low background noise levels. Furthermore, the cost associated with opera- 
ting these facilities is minimal In comparison to the large wind tunnels. 

However, since the sound generated by the model jet has to propagate 
through the shear layer of the free jet before it reaches the microphones, 
this facility does introduce additional problems which have to be overcome. 
These are (l) refraction by mean velocity gradients, and (2) scattering by 
turbulence in the free-jet mixing layer. The refraction phenomenon gives 
rise to a change in the amplitude of the incident wave as well as a modifi- 
cation in the raypath angle. Corrections for these amplitude and angle 
changes therefore need to be established, and these can be obtained theo- 
retically and/or experimentally. Detailed description of the correction 
procedure, derived during the present program, is given in Section 3-1 -3. As 
regards the effects of turbulence scattering in a free-jet facility, all 
qualitative evidence available to date suggests that this is not an impor 
tant effect. However, further work is required to produce quantitative 
results on this phenomenon. This scattering effect is discussed in some 
detail in the later parts of this report. 

Having presented an overview of the various flight simulation facili- 
ties that can be utilized to study the effects of forward motion on jet 
noise, the inflight simulation experiments that were conducted in the 
present program are described in the following subsections. The experiments 
were conducted in the Lockheed Anechoic Free-Jet Facility, which is 
described next. 


2.3 ANECHOIC FREE-JET FACILITY 


2.3.1 Facility Description 

An existing anechoic chamber was modified utilizing as many other 
existing pieces of hardware as possible. Prior to the construction of the 
full scale facility, a one-fifth scale model was built to confirm various 
aerodynamic performance concepts and to aid in the design of the free-jet 
working section as well as the shape of the collector. 

The basic anechoic room was 3-^ m (long) by 3.^ m (wide) by 5.2 m (high) 
between wedge tips. The interior was lined with fiberglas anechoic wedges 
which provided a 99% echo-free environment at all frequencies above 100 
Hertz. The chamber was completely isolated from the rest of the acoustics 
laboratory since it was mounted on massive springs. A spr i ng-tens i oned 
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cable floor, suspended from the walls, provided easy access to the interior 
of the chamber for instrumentation calibration and test section changes. 

A planview schematic of the complete facility is shown in Figure 2.3. 
Starting from the left, air was drawn into the intake, through the honeycomb 
and screens to the contraction, across the anechoic room to the collector, 
through the diffuser, the two right angle corners with turning vanes, through 
the duct silencers to the transition section. The facility was powered by the 
jet ejector, whose exhaust and entrainment flows were diffused through the 
17.1 m long muffler/diffuser section on the right of Figure 2.3. 

The i ntake/cont ract ion section was an existing part of a research smoke 
tunnel, and the ejector/d i f fuser section formed the major part of a now de*- 
activated inlet noise absorption test facility adjacent to the research 
complex. The ready availability of this power source prompted its i ncorpo- 
ration into the free jet facility design rather than a fan. 

Because of the high noise levels generated by the 8.6 cm diameter jet 
ejector, being operated at pressure ratios up to 8 in order to induce flows 
through the working section of up to 75 m/s, a significant amount of acoustic 
treatment was necessary between the ejector and the anechoic room. The design 
criterion adopted was to assume that the total noise generated by the ejector, 
propagated upstream to the collector and then radiated spherically to the 
microphone measuring stations within the room. The amount of attenuation 
provided by the duct was required to be sufficient to reduce the ejector noise 
at the microphones to levels at least 10 dB below the lowest noise levels to 
be measured in the test program over the entire frequency range of interest. 

Of course, this acoustic treatment had to be traded off against aerodynamic 
performance, the aim being to achieve at least 60 m/s free jet velocity. 

To achieve this noise criterion, the tunnel between the collector inside 
the room and the ejector intake was constructed from 1.27 cm plywood inner and 
outer walls separated by wooden studs. The cavities in the floor, walls, and 
roof were all filled with a total of 21,000 kg. of dry sand to minimize the 
flanking (or structure borne) transmission of sound from the ejector. The 
duct was lined throughout with 15*2 cm. of 32 kg/m^ density polyurethane foam. 
Acoustically treated turning vanes were installed in both right angle bends 
and 2.1 m long, low aerodynamic loss, lAC (Industrial Acoustics Company) quiet 
duct silencers were installed downstream of the second turn. 

The intake shown in Figure 2.4 is approximately 2.03 ni x 2.84 m and the 
contraction provides a working section area 0.76 m x 1.07 m, the latter dimen- 
sions being vertical. In order to perform adequate directivity measurements 
of noise generated by the model jets placed inside the free Jet, a minimum 
working section length of 2.9 m was required. This allowed measurements to be 
made in the range 30° ~ 90° from the downstream jet axis. The microphone arc, 
centered at the nozzle exit plane, had a radius of 2.74 m (54 nozzle diameters 
for 5.08 cm. diameter nozzles) and was placed outside the free Jet flow. 


10 




C LU 

o 2: 
</> 


11 



The free jet centerline was oriented at an angle of 10° from the closest 
room wall to allow for free jet spread without having to trim existing room 
wedges. The tunnel diffuser entrance was 1.12 m x 1.1»2 m and diffused to 
1.52 m X 1.83 m in a length of almost 6.10 m. The distance between the turns 
was about 3*05 m and the total length of lined duct from collector to ejector 
intake was about 15-5 m. The total installation was aligned axially using a 
low power laser. 

The air supply to the jet ejector originated from the main 2.07 x 10^ 

N/m compressor which supplied air to all research center facilities. Ejector 
air was controlled by a 15.2 cm. Fisher valve, with a 5-1 cm. Annin valve in 
parallel for low free jet velocities and fine tuning at higher velocities. 

The ejector ducting was 25.^ cm. pipe over the majority of its length but re- 
duced to 20.8 cm just prior to the pylon mount. The ejector air supply 
ducting and diffuser are shown in Figure 2 . 5 . 

For each of installation and minimum blockage (and therefore minimum flow 
disturbance) in the working section, the air supply ducting for the model jet 
was installed axially in the intake/contraction section rather than through a 
swept pylon mounted on the anechoic room wall. The ducting was designed to 
avoid any flow separation within the accelerating free jet flow in the con- 
traction section, a total ly welded construction being adopted for this 
purpose. The ducting was aligned by using a low power laser, placed at the 
end of the collector/diffuser and aimed along the free jet centerline, thus 
ensuring that the model jet flow would exhaust axially in the free stream. 

The cold air supply was regulated by a 5*1 cm. Annin automatically controlled 
valve. A Lockheed-built muffler, and an lAC (Industrial Acoustics Company) 

PRV2 muffler, were connected in series to reduce upstream valve noise to suffi- 
ciently low levels so as not to contaminate any measurements to be made in the 
room. The lAC muffler is shown in Figures 2.3 and 2.k. Downstream of the 
muffler was a 20.3 cm. diameter plenum 6.1 m. long, followed by a 1.8 m length 
of 15.2 cm pipe, and finally, approximately O .91 m of 10.2 cm pipe faired with 
the 15*2 cm. section to preclude any outer flow separations. The nozzles were 
attached to the end of the 10.2 cm. pipe, setting the exit plane approximately 
30 cm. beyond the end of the free jet contraction. 

Four nozzles were specially made for use in this facility, and in the 
wind tunnel (described in paragraph k.2.2): A convergent nozzle and three 

convergent-divergent nozzles designed by the method of characteristics, to 
operate n^inally at Mach numbers of 1.1|, I. 7 , and 2.0. They were turned 
from aluminum bar stock, the external shape being a straight taper from 11.4 
cm. 0/D to 5-1 cm. 0/D over a length of 27.9 cm. to mate smoothly with the 
faired 10.2 cm. I/D air supply pipe, once again to preclude the possibility 
of flow separation. The nozzles were attached with six countersunk screws 
around the circumference at the 10.2 cm. end. The four nozzles are shown in 
Figure 2.6, together with the cross-sectional drawing of the M = 1.7 nozzle. 

Nozzle performance was evaluated by a simple shadowgraph system. 
Unfortunately, the facility used in these tests could not be operated at the 
'ressure ratio required for the Mach 2.0 nozzle, but the M = 1.4 and I .7 
zzles were tested to determine the pressure ratios for optimum shock-free 
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Figure 2.5 Anechoic free-jet facility ejector diffuser 
and air supply ducting. 





Figure 2.6 5.08 Cm. diameter jet nozzle configuration. 



performance. Figure 2.7 shows a series of three shadowgraph photographs for 
each of these nozzles, indicating that the nozzle Mach numbers for optimum 
performance were 1.37 and 1.67 at pressure ratios of 3-06 and 4.73> 
respectively. 

2.3«2 Aerodynamic Performance Evaluation Tests 

2 . 3 . 2 . 1 Flow visualization tests. Initial tests were performed in the 
anechoic chamber using a smoke generator and a tufted wand to ascertain 
visually several important characteristics of the free jet flow in the work- 
ing section and flow induced in the surrounding areas of the room. 

These tests indicated that the Jet appeared to be stable throughout its 
length. There was no discernable snaking or deflection due to the proximity 
of the closest wall or the ceiling. The jet had no tendency whatsoever to 
attach to either of these areas, presumably due to the effectiveness of the 
uneven surface of the anechoic wedges. Stagnation on the collector bell- 
mouth was regularly disposed around the periphery, and no separation was 
detected inside the diffuser. 

In the chamber itself, outside of the free jet flow, the induced veloc- 
ities were quite small, everywhere a good deal less than ] 0 % of the tunnel 
velocity. The highest velocities, approaching 10^ of the tunnel velocity, 
were noted in the standing vortex around the collector bellmouth. The only 
microphone likely to be affected by this flow was located beneath the 
collector at 30° to the downstream axis. 

2. 3 . 2 . 2 Flow surveys using a specially built pressure probe rake. A 
series of initial calibration tests were carried out using a specially de- 
signed pressure probe rake. This rake could be positioned at any axial 
station between the free jet exit plane (at the end of the contraction) and 
the collector bellmouth. It was an easy task to manually perform vertical 
traverses with the rake at selected axial stations. 

The rake itself was made from a 1.22 m length of 7 cm x .12 cm stream- 
lined steel tube. It was supported by brackets at its extremities and 
drilled at 2.5^ cm intervals along its entire length to accommodate 3-18 mm 
probes. In this way, the pressure probes could be arrayed differently for 
each axial station as necessary. Up to 16 total pressure probes and 5 static 
pressure probes were utilized during the test surveys. A second order curve 
was fitted to the readings of the five static pressure probes, distributed 
along the rake, and the room static pressure at the edges of the free jet to 
determine the static pressure at each of the total pressure probe locations. 

Surveys were made at the exit plane of the contraction and then at O. 56 , 
1.12, 1.67» 3nd 2.2^ m downstream. The surveys at each station were made 
from the upper to the lower extremities of the free jet, and the spacing of 
the measurements was chosen accordingly. The rake, positioned at 0.56 m 
axially, is shown in Figure 2.8. In addition, piezometer pressures were 
recorded in the settling chamber of the intake prior to the contraction! arid 
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10 2 upstream of the exit plane. (The static measurement points were 

s^e man^^:: “he ^Ith tSenuLi ",:;^hs"'of"rn“^“h: ™“““ent‘s 
pSo! ^“^“Ta'^:n^t\^=h™ea1™- :rp“ 

manometer display is shown in Figure 2.9. tyP'cal print of the 

a three-dimensional plot of the potential core edae 
(defined here as the 98 ^ velocity point) for a tunnel velocity ofT l/s 

DotenH^lT" Pressure at the jet exit-plane was computed for the 

potential core. The mean excluded the boundary layer areas of the contrac- 
tion section and the centrally-positioned modei jet nozzTe aVr suppric^ 
r?i which*^^? 'ncluded in the averaging is shown cross hatched in Figure 
probes ° ^ positions of the total and static pressure 

tunnercoLmonff ^^'tten for on-the-spot evaluation of 

tunnel conditions from piezometer, temperature and barometric data and 

Xc?iv r differential necessary for a reqi^red tu’nnel 

velocity to facilitate test condition setup. 

eiecto^®wItah;’nf°"'" °^''^'''^tion of free jet nominal velocity. 

l-ecEor pEEsEEre “iiE 

2 . 3 . 2. 3 Conclusion. The aerodynamic performance evaluation tests 
carried out showed that the facility could be confidently used for model 

werJ no aircraft motion on jet noise There 

were no problems associated with free jet stability, collector/diffuser 

performance, flow separation in the intake/contraction or excessive air flow 
circulation in the room. Continuous tests could be carried out at free Jet 
velocities up to 70 m/s. Although free jet velocities up to 75 m/s were 
possible, limitations of available supply air mass flow precluded the 
simultaneous operation of the 5.O8 cm model jets. 


2 . 3.3 Acoustic Performance Evaluation Tests 

accurlcv°nfrh‘° acoustic design criteria and to ensure the 

facilitv measurements, the anechoic free-jet 

P-for»ance evalpation tests, and the 

2.3.3. 1 Anechoic quality Of facility. The fi rst series of cal i brat ion 

ensire'^that^^thl^'’^'^ evaluate the anechoic quality of the facility and to 
warin proposed microphone distance of 2.74 m (54 nozzle diameters) 

exit driver unit placed at the nozzle 

tensity vs. distance plots were obtained with a traversing microphone 
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Figure 2*11 Positions of total and static 
pressure probes for free-jet 
exit plane calibration. 
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arrangement. The microphone was traversed along three different directions 
(30®, 60® and 90® to the downstream jet axis), and for each of these tra- 
verses, measurements were made of the SPL fall-off as a function of distance, 
both for a pure tone noise source and for a one-third octave filtered white 
noise source. Along each of these traverses, the results for discrete 
frequency excitation were slightly inferior to the results for filtered one- 
third octave excitation. However, the results along the three directions 
from the sound source were essentially similar. A typical set of the inten- 
sity-distance plots at various frequencies is presented in Figure 2.13^ It 
can be seen from this figure that the anechoic quality of the facility was 
acceptable down to a frequency of approximately 200 Hz, and that at a 
distance of 2.7^ ni, the microphone was in the acoustic far field at all 
frequencies (above 200 Hz) of interest. 

2.3*3«2 Background and electronic instrumentation noise. The back- 
ground (or ambient noise in the anechoic room was measured with a 12.7 mm. 

BsK microphone (Type 4133) > which had, a dynamic range of 32 to 160 dB. The 
microphone was located in the center of the room, directly under the free-Jet 
test section at a distance of approximately 3 m from the test section center 
line, which corresponded roughly to the 60° measurement position in the jet 
noise test program. The resulting background-pl us- i nstrumentation noise 
spectra at various free-jet (or tunnel) velocities, Vj, are presented in 
Figure 2. 14. 

It should be noted that the total noise floor in the subsequent experi- 
^mental program also contains a contribution from the complete instrumentation 
system used in the data acquisition and reduction process. The magnitude of 
this Instrumentation (electronic) noise contamination unfortunately varies 
from one spectrum to another, depending upon the settings selected. The 
results from the test program (raw data) should therefore be scrutinized 
carefully, and all spectra or parts of spectra, which are obviously seen to 
be contaminated by the total background and measuring instrumentation system 
noise, should be rejected in the analysis of the results. This important 
aspect is elaborated further in Section 2.4, 

2.3«3«3 Big internal noise tests. The '‘acoustic cleanliness" or 
internal noise tests were conducted to produce evidence for the lack of any 
significant contamination, in the measurement arena, from noise sources 
which are normally present upstream of the nozzle exit plane. In the first 
instance, the tests described below were conducted with only the lAC muffler 
in position. It was found that the internal noise contamination was rather 
high, and further silencing of upstream valve noise was necessary. This was 
achieved successfully by connecting the Lockheed-built muffler (used in 
other research programs) in series with the lAC muffler. The tests de- 
scribed below refer to this configuration, which was used in the subsequent 
jet noise experimental program. 

The internal noise tests, with zero tunnel velocity, were carried out 
with the basic instrumentation system, consisting of 12.7 mm. BSK micro- 
phones, a 1/3“Octave analyzer and a level recorder, so that the instrumenta- 
tion noise was kept to a minimum. In order to establish the magnitude of 
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In conclusion, it can be safely stated that the internally generated 
noise from the facility was not significant at least down to Vj/L = 0 32 
11 data for Vj/ao>0.32 obtained from this jet noise rig represented true 
turbulent mixing noise, unaffected by internal noise. tL lod veJdci t^Sata 
would, however, be influenced by background and instrumentation noise Is 
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2 . 3 .^ Point Source Tests 

In theory, a tremendous amount of information can be gained from ooint 
source tests. In particular, in the anechoic free-jet facimr PoTnt 
source tests can be very useful in both cal ibration and basic research test- 
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results and to provide experimental values or confirm theoretical values for 
the amplitude and angle corrections. However, a fairly major effort is re- 
quired to conduct some of the more useful tests and no provision was made 
for extensive point source testing in this prooram. 

Some preliminary point source tests were conducted with the goal of 
determining basic source directivity in the flow, amplitude of source reflec- 
tion at the free-jet shear layer, the far-field amplitude and directivity, 
and the effect of the free-jet mixing region turbulence on scattering of 
sound from the source. The "point" source is shown in Figure 2.17. The 
sound was generated by a commercially available 100-watt audio driver unit 
and propagated through an inverse conical horn section through the right 
angle bend to the 0.6^t cm diameter opening. Output levels of approximately 
115 dB at 1 m from the opening could be obtained from 1 KHz to 5 KHz. The 
source exit point was located as shown in Figure 2.l8 at an axial location of 
20.3 cm. from the free-jet exit plane. 

The first test conducted was a neai field radial traverse for the 
purpose of determining the order of magnitude of internal reflections from 
the shear layer at 90°. Because of differences between the free-jet results 
and inflight data, especially at 90°, it was hypothesized that free-jet 
corrections at 90°, which assume no internal reflections for waves normal to 
the mixing layer, might be incorrect due either to free-jet divergence or 
finite axial velocity gradients in the mixing region. 

The traverse mechanism is shown in Figure 2.17* The microphone was 
equipped with a probe with a right angle bend. The sound signal was sensed 
through static pressure ports. Microphone position was detected by a poten- 
tiometer which drove the X-axis of an X-Y recorder. A typical output is 
shown in Figure 2.19 for four frequencies. The measured amplitude vs. 
distance from the source shown in Figure 2.19 is compared with inverse 
square law for discrete frequency excitation. It is obvious that near the 
source, direct radiation is dominating the sound pressure signal and an in- 
verse square law behavior is observed. However, at a distance some 25 cm. 
from the source, reflections from the exposed free-jet lip have a consider- 
able influence. This problem was not given further consideration since, as 
stated in section 2.3-3. it did not exist for broadband noise and also it 
did not preclude determination of internal reflections from the free-jet 
shear layer, which was the primary objective of this particular test. 

However, because of the strong reflections near the lip line, the source was 
positioned at 10.2 cm. from the lip line in order to change the angle of the 
unwanted reflection and increase the source amplitude at the mixing layer 
interface. A typical result is shown in Figure 2.20 for 2500 and 3150 Hz. 
The solid line is for no-flow, the dashed line for 30.5 m/s and the dotted 
line for 6l.O m/s. The source amplitude remained unchanged during these 
tests. While at 61.0 m/s there were minor perturbations near the free-jet 
lip line, these perturbations were linked to a measurement problem rather 
than a reflection. In the mixing region, the probe microphone sensed high 
broadband turbulent pressure fluctuations, and as a result of non-normal un- 
steady impingement on the static ports, self-noise was generated. The 
ability of the narrow band filters to remove this broadband component of 
noise was inadequate in the mixing regions; thus the measurement was a sum 
of discrete plus random. 
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Figure 2.17 Sketch of point source and microphone traversing 
mechanism located in free jet. 
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Figure 2.20 Effect of free-jet velocity and source frequency on 
sound transmission through free-jet shear layer. 





If, however, the data were examined on either side of the mixing 
region, some general observations could be made. First, no detectable in- 
ternal reflections were observed. On the outside, some small variations 
in level, admittedly influenced by a reflected wave field, were detected. 

Due to the surface reflection problem, no specific conclusion could be de- 
rived for this region. Thus, the only significant conclusion gleaned from 
this test was that internal reflections at the 90° angle were not a 
problem, as expected from the vortex sheet model. 

Another potential problem with free-jet forward flight effects testing 
is the scattering of sound, generated within the potential flow region, as 
it propagates through turbulent eddies in the mixing layer. This problem 
has been analyzed by Mani (ref. 27), who looked at the change in directivity 
for a highly directive source. He showed considerable diffusion in the far 
field sound pressure directivity for frequencies as low as ^o^s/^o“^ where 
in this case r^ was the radius of the source region. The magnitude of dif- 
fusive scattering was controlled by the parameter DR where R was the far- 
field radius and D was a diffusion constant which is defined as 1/4 of the 
mean square angular deviation of the ray, per unit length of the ray path, 
due to scattering. The diffusion constant is shown, for high frequency, to 
be approximately represented by 


DR = U^R/a^) f(o)oL/ao) 

where is the turbulent Mach number and fCoiQl/a^) is inversely related to 

the length scale of turbulence. 

For values of DR = 10 ^ and 5x10 ^ at frequencies t^o'^s^^o ” ^ 
calculated diffusion is very significant and in some cases a very peaky angu* 
lar distribution is diffused to the extent that it nearly appears to have 
originated from a point monopole. Thus, the diffusion or turbulent scatter- 
ing problem is one which deserves attention in free-jet acoustic studies. 

Another approach at examining the turbulent scattering problem con- 
siders spectral broadening of an original discrete frequency signal. As the 
discrete tone wavefront propagates through the turbulent eddies, the speed of 
propagation of the wave varies. This random variation in propagation speed 
along the wavefront results in a narrow-band random modulation in the ampli- 
tude and phase of the signal which manifests itself as spectral broadening. 
Many studies have been conducted In this area of research [Rudd (ref. 28), 
Brown and Clifford (ref. 29)]; however, experimental data on free-jet facili- 
ties have not been available until recently. Foley and Paterson (ref. 30) 
discuss a test in the United Aircraft anechoic free-jet facility. They dis- 
play a test result for a 24,000 Hz signal propagating through the shear 
layer of their free-jet facility for 6l.O m/s tunnel velocity. Spectral 
broadening in this instance appears to be fairly significant although no 
quantification is given. They state that the effect is not significant 
until the wavelength is less than the Integral scale of turbulence. 


Because of source limitations it was not possible to conduct scattering 
tests above 5000 Hz. The source directivity was essentially spherical at 
that frequency, thus an examination of angular diffusion was not possible. 
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However, an examination of the frequency spectrum in the far-field was under- 
taken with and without flow. The spectrum was examined through a 2 Hz narrow 
band filter. A typical result at 5000 Hz and 52.5^ angle to the jet axis is 
shown in Figure 2.21, first for no flow and then at 6l.O m/s tunnel velocity, 
for a source positioned 10.2 cm above the lip line and 20.3 cm downstream. 
At 20,3 cm. in this facility, the integral scale of turbulence corresponded 
approximately to wavelengths of 1.3 cm to 3*8 cm or frequencies of 8000 Hz 
to 2^,000 Hz at the position where the ray propagated through the mixing 
layer. Examination of the two narrow band spectra shows that there is no 
detectable difference and in fact further examination shows that the spectral 
shape is that of the 2 Hz filter. 

Thus, it can be confirmed that at wavelengths greater than the scale of 
turbulence, spectral broadening is not a problem in free-Jet facilities. 
Further work is required, however, for higher frequencies. 


2.4 EXPERIMENTAL PROGRAM AND PROCEDURE 

The influence of forward motion on unheated Jet exhaust noise was 
examined in the present program by measuring the turbulent mixing noise in 
the far field from three 5*08 cm diameter nozzles, operated in the anechoic 
free-jet facility. The free jet to nozzle area ratio was 400. The nozzles 
employed were the M = 1 convergent nozzle, and the M = 1.4 and 1.7 convergent- 
divergent nozzles. In order to avoid the contamination of data from shock- 
associated noise, these latter were operated at their design pressure ratios 
only. Nine microphones were placed on an arc of radius 2.74 m (R/d = 54) , 
centered at the nozzle exit plane. The measurement angles were 0m = 3O°(7i°) 
90° relative to the downstream jet axis. 

The experimental program was carefully planned to yield results at (1) 
constant jet efflux velocity (Vj/ao) with varying tunnel velocity (Vj/a^) , 
and (2) constant tunnel velocity with varying primary jet velocity. Forty 
test conditions (combinations of Vj/aQ and Vj/ao) were initially chosen. The 
nominal values of jet exit velocity and tunnel velocity were 

Vj/a^ = 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.17 and 1.345 (8 values); 

Vj/a^ = minimum, 0.05, 0.10, 0.15 and 0.20 (5 values). 

The minimum tunnel velocity, Vjq, is defined as the tunnel velocity which 
is provided by the ejector action of the primary jet in the absence of any 
additional velocity supplied by the facility ejector itself. Its magnitude, 
therefore, increases as Vj/3q increases. 

For each combination of Vj/a^ and Vj/aQ, the static pressure ratio 
Pj^/pf was determined by using the analysis developed in Appendix 2B. The 
values of various parameters defining the test conditions are tabulated in 
detail at the beginning of Appendix 2C. 
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Figure 2.21 Effect of propagation through shear layer 
on spectrum shape of discrete tone source 
of 5000 Hz discrete frequency for narrow 
band (2 Hz filter) spectra measured In 
far-field at 52.5° to jet axis. 
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2.4.1 Test Procedure 

The quantities to be measured in setting up the test conditions are 
defined in the sketch shown in Figure 2.22. In order to set the experiment 
at the desired combination of Vj/a^ and Vj/a^ (for the convergent nozzle) or 
Mj and Vj/aQ (for the two con-div nozzles), the following sequence was 
adopted: 

Tunnel : 

1 . Sped fy Vj/ao 

2. Read p^, T^ 

3. Calculate 8 q (= aj) 

4. Calculate Vj 

5. Compute (p| “ pj) using tunnel calibration computer program 
on HP-65 

6. Set tunnel at (p^ - Pj) to obtain desired Vj/a^^ 


Jet : 

7. Specify or Mj 

8. Compute Pp^/py using relationships derived in Appendix 2B 

9. Read py 

10. Calculate p^^ 

11. Set jet at pp^ to obtain desired Vj/a^ or Mj 

12. Check if (p^ " Py) has altered, and adjust if necessary. 

During each test, all pressure and temperature measurements were re- 
corded; subsequent to the test program, these measurements were used to 
compute the exact jet and tunnel operating conditions, which are tabulated 
in Appendix 2C, for every test conducted. 


2.4.2 Data Acquisition 

In t6e present experimental program, nine 12.7 mm, BSK microphones Type 
4133 with FET cathode followers Type 2619 were mounted on the microphone arc 
at 7i° intervals from 30^ to 90^ to the downstream jet axis. The responses 
were recorded simultaneously on a l4-channel Honeywell FM tape recorder at 
305 cm/s. In order to obtain the 1/3“OCtave spectrum from 200 Hz to 40 KHz 
using a Hewlett-Packard real time l/3“Octave audio spectrum analyzer, the 
tape speed was reduced to 76.25 cm/s on playback. The l/3“Octave levels 
were then recorded on an incremental tape recorder for subsequent detailed 
analysis using a data reduction program developed for use on the Univac 4l8 
digital computer. This program incorporated the rrucrophone frequency 
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Figure 2.22 Measurements for setting test conditions. 
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response corrections and atmospheric absorption corrections. The results, 
were finally displayed in the form of tabulated one-third octave band sound 
pressure levels over the frequency range from 200 Hz to 40 KHz> and the 
overall levels were also computed over this frequency range. 

2.4.3 Data Qual i ty 

All acoustic data acquired in the present experimental program were 
carefully scrutinized in order to ensure that the results which are selected 
for detailed analysis would not be contaminated to any significant extent by 
contributions from any other noise sources. The data from the entire test 
program are given In Appendix 2C for future reference. In using these 
data, the following points must be considered at all times: 

(1) There is no significant contamination from rig internal noise 
sources, as established in Section 2.3.3- 

(2) The two major sources of extraneous noise are the tunnel background 
noise and the measuring instrumentation electronic noise. As expected, the 
influence of these sources is maximum for low Vj/aQ, high Vj/aQ test condi- 
tions. In general, when present, the effect of background noise is noticea- 
ble at low frequencies, whereas the presence of instrumentation noise 
results in a slight lift in the spectrum levels at the high frequency end. 

(3) Although the supersonic nozzle calibration tests (shadowgraph tests) 
revealed no significant shock pattern in the jet flow at design pressure 
ratios for all nozzles, close examination of the acoustic data indicates that 
in the case of the M = 1.4 nozzle, the spectra are slightly contaminated by 
shock-associated noise at the higher frequencies. This indicates slight 
imperfections in the machining of the nozzle contour. 

(4) At high tunnel velocities, the spectra measured by the 6^ = 30° 
microphone are affected by the aii flow generated noise over this microphone 
as described in Section 2. 3-2.1. 

(5) The spectra measured at 0m = 75° contain a slight ripple due to 
room reflection, which was not detected in the earlier acoustic calibration 
tests. However, this effect is consistently present, and therefore it does 
not influence the final conclusions derived from the data analysis. 

The results which are presented in the rest of this report are essen- 
tially free from any of the above limitations; hence, any conclusions or 
scaling laws derived from them are considered to be accurate. 

One further point regarding data presentation needs to be made here. 
Although the measured spectra contain minor variations in level about the 
mean variation from one 1/3-octave band to the next, the spectra presented 
in this report are shown as continuous spectra, purely for convenience. A 
typical example, showing the actual measured 1/3-octave band levels and the 
corresponding continuous spectrum, is presented in Figure 2.23. It can be 
observed that the presentation of data in this fashion takes no unjustifi- 
able liberty, and indeed it simplifies the spectrum plot considerably. 
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especially when it is necessary to plot several closely spaced spectra on 
the same figure. 


2.5 MEASURED (OR UNCORRECTED) RESULTS 

The turbulent mixing noise results from the i nf 1 i ght s imulat i on experi- 
ments, conducted in the anechoic free-jet facility, are first presented 
without applying any facility corrections. These corrections were mentioned 
briefly earlier, and their magnitudes as well as the correction procedure 
are derived in the present program in Section 3-1 • The corresponding cor- 
rected results will be presented and discussed in the next section. The 
major effects of forward motion on the directivity and spectral characteris- 
tics of turbulent mixing noise, derived from the uncorrected results are 
discussed below. 


2.5.1 Overall SPL Results 


The influence of tunnel velocity (Vj/a_) on the directivity of overall 
SPL at various jet efflux velocities in the wide range from Vj/a^ = 

0.1* to 1 . 31*5 is shown systematically in Figure 2. 24. It should be noted 
that the angle is the measurement angle relative to the downstream Jet 
axis, and not the emission angle. For a fixed measurement angle 9^,, the 
angle at which sound was emitted from the jet, 6g, varies as the tunnel 
velocity is varied. Conversely, in order to examine the changes in noise 
radiated at a fixed emission angle Og with Increasing tunnel velocity, the 
measurement angle 0 must be varied appropriately. The relationships be- 
tween these two angles form a part of the correction procedure that will be 
described in section 3-1 • 


The directivities shown in Figure 2. 24 contain remarkably 1 i ttle 

fixed value of jet velocity Vj/a^, the directivities at 


I n- 


scsttcr. For a * vvj.wn... w. j — , o u . , -ru* 

creasing tunnel velocities appear virtually parallel at first sight. This 
observation agrees qualitatively with the results of the free-jet experi- 
ments conducted at NASA-Lewis (ref. 21). However, a closer examination also 
reveals that, in general, the reductions in OASPL with forward velocity are 
slightly angular dependent, being a little higher as the measurement angle 
approaches the jet axis. 


In order to provide a quantitative description of the reductions in 
OASPL at all measurement angles considered, the results were correlated on 
the relative velocity (VreL“^J”^T^ basis, as follows. 


The overall intensity of turbulent mixing noise at measurement angle 
can be written to scale according to 


I (6m) 


y m y n 

''rel ''j 


( 2 - 1 ) 


where the exponents m and n are functions of e„,. The corresponding OASPL is 
therefore given by 
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®m 

(a) Vj/a^ = 0.4 



(b) Vj/ao = 0.5 


Figure 2.24 Effect of tunnel velocity on the directivity of measured 
(uncorrected) overell SPL at various jet exit velocities 
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Figure 2.24 Continued. 
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( 2 - 2 ) 


OASPL(e„) 

m 


10 logio 


fw"* uni 

^''rel ’'jJ 


the reduction In OASPL fro. .inl.u. 
city kVjq) to any other tunnel velocity (V^) can be written as 


AOASPL(e^) = 10 logjo 


Vj-Vto 

Vj-Vt 


m 


(2-3) 


o^'^KuredTAskL ‘"f '1,°'^ 

'0 log,„ ((Vj-Vto)/V„eJ function of the velocity parameter 

the vI^ocnrp“a^'m«e"'!n ngvrells"::^ 

::s: rsciJu-fstn^kiTr i 

except at 6 =30° where as mpnr' ^ ounds of the experimental accuracy, 
• L m * wnere, as nientioned previouslv fh^ ^ 

microphone produced undesirable pressure fluctualio^s in 


tr™:'!'';„'^! ''r'f''”" — >ty exponent 


measurement angle 0 is Dresenf'pW exponent m with 

figure are the values obtained f »-k 2*26. Also shown in this 

Von Glahn, et at (ref 21) It^sho experiments reported by 

p-jcei -f-;-- 

p.lr rrl-ret^rmerThl 

be expected in view of the differences in the two facNilies Ld " 

2 . 5.2 l/3“0ctave Spectral Results 

three';::a::Ji:::nTangier"' 

range are shown in Figures 2.27, i 28 and 2 M ' velocity 

sented in these figures is suff clenr,r„ro:-a' bf'- 

recalled that some of the snertm do* '? Provide a complete picture. It Is 
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of frequency. On the other"hand'"'^tt^[s^"’ fairly independent 

other hand, ,t is misleading to derive any conclusions 
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Relative velocity scaling of measured 
(uncorrected) overall SPL reductions 
3t various measurement angles. 
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Figure 2.25 Continued. 
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Figure 2.25 Continued. 
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Figure 2.26 Relative velocity exponent from measured 
(uncorrected) overall SPL results. 



1 /3-OCTAVE SPL 



Figure 2.27 Effect of tunnel velocity on measured 
(uncorrected) spectra at low jet 
velocity (Vj/a^ = 0.6). 
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at this stage from the measured spectra at angles other than 90° to the jet 
axis, since the emission angle corresponding to a fixed measurement angle 
could vary significantly with varying tunnel velocity. The spectra at 
0m =60 and 37-5° are presented here purely for future reference. As a 
general observation, however, we may note in passing that at small measure- 
ment angles, the reduction in sound pressure levels with increasing tunnel 
velocity is frequency dependent, and the magnitudes of these reductions 
decrease progressively as the observed frequency increases beyond the peak 
frequency. 

It remains to be seen how the spectral characteristics observed above 
from the uncorrected results are modified by the facility corrections. This 
aspect is examined in the following section. 


2.6 CORRECTED RESULTS 


In the previous section, the effects of forward motion on the direc- 
tivity and spectral characteristics of jet mixing noise were examined in a 
a preliminary manner where the measured results were not corrected for any 
of the complications introduced by testing the model jet in a larger co- 
flowing stream, with the microphones placed outside in a stationary medium. 
An adequate understanding and accurate quantification of these facility 
corrections are vital to the success of a free-jet facility in simulating 
forward velocity effects on jet noise. Considerable effort in the present 
program was directed towards this aspect of inflight simulation, and the 
complete details of various phenomena involved, together with a systematic 
data correction procedure, are given in section 3-1 of this report. In 
essence, the correction procedure takes proper account of source distribu- 
tion effects in a jet flow, the downstream convection of sound waves by the 
tunnel flow, and the refraction of sound caused by the free-jet shear layer. 
The refraction effect causes a change in the ray path angle and a change in 
the sound pressure amplitude as the sound waves propagate through the shear 
layer. The correction procedure incorporates all these effects in a 
realistic manner, and using the measured or uncorrected results at fixed 
measurement angles, 0^, it yields results corrected to constant emission 
angles, 0^, for an observer moving with the jet. 

In the present program, the measured results at three jet exit 
velocities (Vj/aQ = 0.6, 0.9 and 1.3^5) were subjected to the correction pro- 
cedure, and the corresponding corrected results at fixed emission angles are 
presented in this section. The spectral results are examined first. 

2.6.1 Corrected 1/3 -Octave Spectra 

The corrected one-thi rd octave spectra corresponding to the measured 
spectra discussed in the previous section are presented In Figures 2.30, 

• 31, and 2.32, for Vj/aQ = 0.6, 0.9, and 1.3^5, respectively. In order to 
obtain the corrected spectrum at some emission angle, 6 , the measured 
spectra have to go through an interpolation procedure, both in frequency and 
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1 /3-OCTAVE SPL - dB 
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Figure 2.31 Effect of tunnel velocity on corrected spectra 
at medium jet velocity (Vj/Sq = 0.9)- 
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Figure 2.31 Concluded. 
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In anale as wi H be described in section 3.1. At low frequencies and low 
Q Jh ^procedure requires measured values at angles 9^ smaller than those 
'uUd ^'Experimental program, and hence. 

have a lower limiting frequency; it is not possible to obta n corrected 
levels at frequencies lower than this limiting frequency. In the spectra 
p^sented in the figures, the lower limiting frequency, when present, .s 
identified by a vertical line in most cases. 

The general features exhibited by the corrected spectra at constant 
mission TngUs (6^) -e quaUtative,, similar “ 

„»asured spectra he nSi' (or'elf^^tL) correc- 

c/vnp nusntltStlVCdlfiCrcnCSS* \ • •Ui — 

tions^at all tunnel velocities and frequencies are found to be ; 

Thr reductions in 1/3~octave levels with increasing tunnel velocity are vir 
I^al if identical to those noted previously at 9^ = 90°. although the absolute 
levels ft all Eel velocities are a little higher due to the source oca- 
tion corrections. In contrast, at smaller emission angles, the shear 1 y 
refraction corrections do become significant. Their magnitu . 

tSe tunned ^^'ocity increases. The consequence of this effect is that the 
reductions obtained from corrected spectra at constant Og ^ ^ '' . 

the rLuctions observed in uncorrected spectra at the same value of constant 
T tn other words, at low emission angles, the relative velocity effect is 
laraeffof the corrected results than for the corresponding uncorrected 
results at the same values of measurement angles. The reductions in correct 

polL Uet. with tunnel velocity were scaled ‘^7''“7ird nn Ire 
basis, as done previously for the uncorrected results, and the findings 

described in the following section. 


2.6.2 Corrected Overall SPL Results 

Th« ovd-rall sound pressure levels for the corrected spectra were ob- 

frequency limit imposed by the correction procedure, as mentioned above. 

The effect of tunnel velocity Vy/a^ on the directivity 

^ ^tl r?s'L»?rd he'rht constant 

^i ter I„rthfr:?IIe;:epresents a true 

irs.nlrrrrili^iltirt^irsrtl^uctr :r:t; Ilissl. anqles con- 
sidered hert. The magnitudes of the reductions increase slightly as 
observer moves from 0e = 9O° towards the downstream jet axis. 

In order to obtain the relative velocity exponent m for 
results the OASPL reductions at all emission angles were plo e ag 

the velocity parameter 10 logio ^ aLuEf Lta^sufected to 

iSe Ifciluy^orreEs wf" limi ied'wi thi n the scope of the present program 
Ife fumblrof^Safa'oints in the figure, through which th s o 
exponent m is plotted, is unfortunately limited. Neverthele . 
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(b) Vj/ag = 0,9 


Figure 2.33 


Effect of tunnel velocity 
overall SPL at various jet 


on the directivity of corrected 
exit velocities. 
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Figure 2.}k 


Relative velocity scaling of corrected overall 
SPL reductions at various emission angles. 
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results are considered to be fairly accurate. At 90° to the jet axis, the 
overall intensity is proportional to 5-5 powers of the relative velocity, a 
result which is closely in agreement with the scaling laws derived from 
jet flow measurements in Section 

Finally, the variation of resulting relative velocity exponent m with 
emission angle 0^ is plotted in Figure 2,35* For comparison, the exponent 
values obtained in Section 2.5 from uncorrected data are also shown in this 
figure. In order to compare the forward velocity effects on jet noise ob- 
tained from the present experiments with the results of other investigations, 
the latter are also included in the figure. In particular, the published 
inflight effects studies considered here are (1) the free-jet experiments 
conducted by NASA-Lewis (ref. 21), (2) the wind tunnel experiments conducted 
by NGTE (ref. 23), and (3) the full-scale flight results reported by Rolls- 
Royce (ref. l). Figure 2.35 thus provides a simplified overview of the 
inflight effects on jet noise on the relative velocity basis from all types 
of facilities. The major facts and implications can now be discussed in some 
detai 1 . 


(1) The first major observation can be derived by comparing the rela- 
tive velocity exponents obtained in the present experiments from the 
corrected data and the uncorrected data. The facility corrections, designed 
to convert results from a free-jet simulation to a corresponding wind tunnel 
simulation, are negligible at 90° to the jet axis, but produce a significant 
difference at lower angles in the rearward arc. 

(2) In the NGTE wind tunnel experiments, the procedure for deriving 
the relative velocity exponent m included the so-called dynamic effect 
correction (i.e., an amplification due to relative motion between tunnel 
flow and stationary jet), which was not included so far in the present 
data manipulations. In order to obtain a true comparison with the NGTE 
results, this correction was also applied to the present results by 

us i ng 


AOASPL(0e) 

wi th 
dynamic 
cor rect ion 


1 + Mjcos0g 
1 +MjqCOS0^ 

wi thout 
dynamic 
correction 


A0ASPL(6e) - lO'^ogiQ 


(2-4) 


where 


AOASPL(Oe) 

wi thout 
dynamic 
correction 


10 logio 


Vj ~ Vto 
Vj-Vt 


m 


(2-5) 


The overall SPL reductions with dynamic corrections are therefore lower than 
the overall SPL reductions without dynamic corrections. They were correlated 


69 





REFERENCE [23J 
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PRESENT (CORRECTED 
WITHOUT DYNAMIC EFFECT) 

PRESENT (CORRECTED 
WITH DYNAMIC EFFECT) 


STATIC 


REFERENCE [21 ] 
(UNCORRECTED) 


/^ULL-SCA^ 

^Zflight 



EMISSION ANGLE 


- DEGREES 


MEASUREMENT ANGLE 9^ - DEGREES (FOR UNCORRECTED RESULTS) 


Figure 2.35 Comparison of relative velocity exponents from various 
inflight simulation experiments and flight tests. 



on the relative velocity basis at each emission angle 0^, and the resulting 
relative velocity exponent values are given by the broken curve in Figure 
Q seen that the effect of this dynamic correction is zero at 

0e 90 * and the magnitude of the effect increases as 6^ decreases. 

Comparison of these modified exponents with the corresponding NGTE 
results in Figure 2.35 shows that the agreement is good at all emission 
angles. The small differences are perhaps to be expected for two reasons: 
(a; the present results were corrected to take proper account of realistic 
source location, whereas no such correction was applied to the NGTE results 
(where the sources were assumed to be located at the nozzle exit plane)- 
(b) the microphones in the NGTE tests were placed at a sideline distance'of 
twenty-two nozzle diameters from the jet axis, which is not large enough to 
provide results that approach true far-field radiation characteristics! 


The good agreement between the results of the present free-jet experi- 
ments, carefully corrected to yield results which would be expected from a 
corresponding wind tunnel simulation, and the results from the NGTE wind 

which Utilized a large test section to nozzle area ratio 
(- 5000 ;, provides three important implications: (a) the effect of turbu- 

lence scattering in a free-jet mixing layer, which was not included in the 
data correction procedure, does not appear to be significant; for if it was 
significant, then the agreement between the free-jet experiments and the wind 
tunnel experiments, where the scattering problem does not arise, would not be 
as good as we see here; (b) the free jet to nozzle area ratio of hOO used in 
the present experiments appears to be adequate for accurate evaluation of in- 
flight effects on jet noise characteristics; (c) a free-jet facility is capa- 
ble of simulating many of the inflight effects on jet noise, providing ade- 
quate and accurate facility corrections are applied to the measured results. 


^3) Finally, the comparison of relative velocity exponents from in- 
flight simulation experiments with the envelope of exponents obtained from 
full-scale flight results is rather discouraging. At low angles to the jet 
exhaust, it appears that the flight simulation results are in reasonable 
agreement with the actual flight results. At larger angles, however, the 
reductions in noise levels with forward motion observed in the model simula- 
tion experiments are larger than those measured in flight tests. At 0 =90°, 
there is little or no change in the flyover noise levels, whereas the free- 
jet experiments indicate significant noise reductions which scale on the 5 5 
powers of the relative velocity. 


The obvious question to be raised at this stage is, what is the reason 
tor this discrepancy between the flight simulation results and the full- 
scale a i rcraft flyover results? The first possibility stems from the fact 
I present experiments, although the jet velocity regime was simu- 

lated adequately, the effect of jet heating was not considered; in contrast, 
all flyover results are obtained from heated jet flows. Hence, an exact 
caparison for matching jet operating conditions should be carried out after 
obtaining results from heated flight simulation experiments. Having stated 
this, however, the acoustic scaling laws from heated jets, derived in 
Appendix ifC, suggest that although this possibility will bridge the gap 
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between model simulation and flight results considerably, it might not be 
able to provide the complete explanation. 

In attempting to resolve this discrepancy further, it is readily admit- 
ted that there are subtle differences in the two sets of experiments: in 

the flight simulation tests, there is no relative motion between the source 
and the observer, whereas in the flyover tests, the noise source is being 
convected relative to a stationary observer. However, further considerations 
also show that none of the features associated with source motion (for 
example, eddy convection velocity effects, source acceleration effects, etc.) 
will affect the results at 90° to the direction of motion. It must therefore 
be concluded that the flight results considered so far may not represent pure 
turbulent mixing noise, and are likely to be contaminated by engine internal 
or other noise sources. 


2 . 6.3 Prediction of Stat i c-to-Fl i ght Noise Reduction 

We have previously alluded to the point that although the concept of 
relative velocity exponent is very convenient in semi -empi ri cal correlation 
schemes, it does not provide any physical explanation of static to flight 
noise reductions, except at 90° to the jet axis where It is a measure of the 
changes in equivalent source strength. Hence, an attempt was made to 
correlate the measured {with facility corrections) OASPL reductions at all 
emission angles in terms of a theoretical result for convective amplification 
due to eddy convection, derived in Section 3- 

The starting point for the theoretical correlation formula is that the 
Intensity at 0e = 9O° is proportional to The intensity at any other 

emission angle 6g is then given by adding the eddy convective amplifica- 
tion and jet motion amplification (resulting from the relative motion 
between tunnel flow and stationary jet) to this basic dependence, which is 
the source alteration effect. The overall SPL can therefore be expressed as 


1 + MjCOs6g I I 


0ASPL(6e) 


10 log 


10 


(Vj-Vj)5-5 (CAvj_v^) 


where CA is the eddy convective amplification given by 

(1 - M.cose-)^ 

CA s — ; , ( 2 - 

{(1 -M(.cos6g)^ + ' 

and the subscript Vj - Vj denotes the velocity at which the eddy convection 
Mach number Me is to be evaluated. This convect i ve ampl i f i cat ion result is 
derived from geometric acoustics in Section 3, and is valid outside the 
zone of silence. The reduction in OASPL from minimum tunnel velocity, Vyo. 
to any other tunnel velocity, Vj, can be written as 
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AOASPL(0e) 


[|Vj-Vto 

5.5 

CAVj - \IjQ 

1 +M-j-cos0g 


1 Vj - Vt 


CAVj - Vt 

1 + M-pQCOS0g 



(2-8) 


In using equations (2-7) and (2-8), the eddy convection velocity was assumed 
to be 0.67 times the applicable relative velocity, and a was taken to be 0.3, 
Although these values are consistent with the values obtained from static jet 
flow measurements by several investigators, their magnitudes for Jets in a co- 
flowing stream need to be derived experimentally. 


For all measured (with facility corrections) OASPL reductions, the corre- 
sponding theoretical OASPL reductions were calculated, and the comparison of 
results at all emission angles 0^ (including 90°) greater than the cone of 
silence angles where 0^ is given by 


is shown 


0^ = cos 


in Figure 2.36. 


1 


1 + 0.67 




(2-9) 


The calculated noise reductions are in close agreement with the measured 
(with facility corrections) noise reductions, the majority of the points 
agreeing within ±1/2 dB. 


The good agreement observed here provides an important indirect impli- 
cation. The calculation formula used here Is basically similar to the formula 
used by NGTE; although there are significant physical differences between the 
two, the numerical differences over the range of Vj/aQ considered here are not 
significant. The fact that both the wind tunnel results and the free-jet 
results correlate well with the prediction formulae suggests that the corrected 
results from the present free-jet experiments are compatible with the results 
obtained from the wind tunnel experiments. This provides further Indirect 
evidence to our previous conclusion that a free-jet facility is very suitable 
for a fundamental study of forward motion effects on jet noise. 


2.7 SUMMARY AND CONCLUSIONS 

The effects of forward motion on the characteristics (both directivity 
and spectral) of turbulent mixing noise from jet exhausts were examined by 
conducting inflight simulation experiments in the Lockheed anechoic free-jet 
facility. The highlights of the experimental program and the major con- 
clusions are as follows: 

( 1 ) The facility was capable of providing tunnel velocities up to 75 m/s 
with a test section of 0.76 m x 1.07 m. Prior to the jet noise experiments, 
the facility was subjected to detailed aerodynamic and acoustic performance 
evaluation tests. 
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Figure 2.36 Comparison of measured (with facility corrections) 
OASPL reductions with calculated OASPL reductions. 
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(2) The aerodynamic performance tests revealed that there were no 
problems associated with free-jet stability, and the air flow circulation 
velocities in the anechoic chamber were acceptable. 

(3) The acoustic performance tests established that (a) the facility 

was anechoic for frequencies down to 200 Hz, (b) the facility background noise 
in the measurement arena was low, and (c) the jet noise results would not be 
contaminated by rig internal noise at least down to Vj/ao = 0.32. 

(4) Limited calibration experiments using a point source established 
that in the frequency range from 1 KHz to 5 KHz, no detectable internal re- 
flections in the free-jet test section were observed for tunnel velocities 
up to 6l.O m/s. It was also confirmed that at these frequencies, the effect 
of turbulence scattering in a free-jet facility was not significant. These 
point source experiments need to be extended to higher frequencies. 

(5) In the jet noise experimental program, the free-jet to nozzle area 
ratio”^ was 400. The nozzles (diameter = 5-08 cm.) employed were M = 1 conver- 
gent nozzle, and M = 1.4, 1.7 convergent-divergent nozzles. The measurement 
angles were 6^ = 30^ (7i°) 90^ relative to the downstream jet axis, and the 
microphones were placed at 54 nozzle diameters. Acoustic measurements were 
conducted at forty different test conditions (combinations of Mj/a^ and 
Mj/a^); eight values of jet exit velocity were chosen in the range 

0.4 < Vj/aQ < 1 *3^5> and the tunnel velocity was varied up to Vj/aQ = 0.2. All 
acoustic data were carefully scrutinized prior to detailed analysis, and all 
data likely to be contaminated by extraneous noise sources (background/ 
instrumentation/shock noise) were not utilized. 

(6) The uncorrected results (i.e., data to which no facility correc- 
tions were applied) were found to be in good agreement with the results from 
published free-jet experiments. 

(7) The measured results were subjected to a systematic data correction 

procedure, which was derived in the present program, and which converted the 
results from a free-jet facility to the corresponding results that would be 
obtained in a large-scale wind tunnel simulation. In essence, the correction 
procedure took proper account of source distribution effects in a jet flow, 
the downstream convection of sound waves by the tunnel flow, and the refrac- 
tion of sound caused by the free-jet shear layer. Using the measured or 
uncorrected results at fixed measurement angles » ®m» the procedure finally 
yielded results corrected to constant emission angles, for an observer 

moving with the jet nozzle. 

(8) In general, the magnitudes of the facility corrections were negli- 
gible at 90^ to the jet axis, but produced a significant effect at low 
angles in the rearward arc. 

(9) At 6^ = 90^, the effect of tunnel velocity on the corrected spectra 
was virtually independent of frequency, and the spectra at various tunnel 
velocities were nearly parallel. This implied constant reduction in equiva- 
lent source strength at all frequencies. At smaller angles to the downstream 


+5ee Footnote on page 76. 
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jet axis, the noise reductions with forward velocity were slightly frequency 
dependent; at 0^-40°, the magnitudes of these reduct ions decreased slightly, 
but consistently, as the frequency increased. 

(10) The corrected overall SPL results indicated that the magnitudes 

of the inflight noise reductions increased as the observer moved from Og = 90° 
toward the downstream Jet axis. The reductions at all emission angles were 
scaled on the relative velocity (Jet velocity minus tunnel velocity) basis. 

At 0g = 9O°, the overall intensity was found to be proportional to 5.5 powers 
of the relative velocity, a result which agreed closely with the theoretical 
velocity dependence, I that is derived purely from source alteration 

considerations in Appendix 4C. 

(11) Comparison of the relative velocity exponents m, obtained from 
the present corrected results, with the exponent values obtained in the NGTE 
wind tunnel experiments (where the test section to nozzle area ratio was 
5000), showed that the agreement was good at all emission angles.' From this 
comparison, it is concluded that (a) the effect of scattering caused by 
turbulence in the free~Jet mixing layer did not appear to be significant; 

(b) the free Jet to nozzle area ratio of ItOO used in the present experiments 
was adequate for accurate evaluation of inflight effects on Jet noise, and 

(c) a free-Jet facility is capable of simulating most of the inflight 
effects on Jet noise, providing adequate and accurate facility corrections 
are applied to the measured results. 

(12) The comparison of results from inflight simulation experiments 
with full-scale flight results shows that although the agreement is reason- 
able at low angles to the jet exhaust, there are significant discrepancies 
at larger angles. At Bg = 90° , there is little or no change in the flyover 
noise levels, whereas the free-Jet experiments indicate significant noise 
reductions. All attempts to explain this discrepancy in terms of source 
motion (relative to a fixed observer) effects lead to the conclusion that 
the flyover results at 90° will not be affected. It is therefore concluded 
that the flight results considered here may not represent pure turbulent 
mixing noise, but are likely to be contaminated significantly by engine 
internal or other noise sources. It is recommended that the "acoustic 
cleanliness" of all flight results be examined thoroughly, both for existing 
flight data as well as any future flight test results. 

( 13 ) Finally, a theoretical formula for the prediction of inflight 
noise reductions at emission angles outside the so-called zone of silence 
was presented, and the correlation between calculated and measured OASPL 
reductions was obtained to an accuracy of ±1/2 dB. 


~^The area ratio of 400 in the present free-jet experiments is based on the 
nozzle exit area for the 5.08 am. diameter primary jet, and is not based on 
the cross-sectional area of the 10.2 am. upstream air supply pipe. This is 
justified, since there is no flow separation over the outer surface of the jet 
nozzle. For a suppressor nozzle configuration, however, it may not be possible 
to consider the area ratio in this manner. In particular, for highly segmented 
nozzles, the area ratio should not be based on the effective area at the 
nozzle exit plane. In such cases, it is necessary to examine the spread rate 
of the jet flow in order to ensure correct aerodynamic simulation. 
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3, ACOUSTIC PROPAGATION THEORY 


exhaul; ““''o" °" je< 

theory is crystaUiL " E^a'^n^aroSf;'’ ■■?'?, Propagation 
to specify the features required of fl^^h? ™ °f fhght BmulaUon-, that is, 
relate acoustic measuremi^s rea ' fj J™ '» 

those taken under flight t^^t cLmL " facilities with 

of thl^^nuSstnfni^M';.;;””;-'': idealized models 

basic features ^^n be a^; r^' V"" '“'allzed models, 

in a straightforward way then the relations can be established 

realistic, in a step-byltin J.sM ^ Progressively more 

of each real effect can be assessed fn'; 

known?^earsound'p^opagijIon^ uHdeV'??' '"^»^ence of all 

ditions (e.g. turbulencHraE.!! ^ simulated flight con- 
However, when the investigation is^severelv UmirnT^^ clearly desirable, 
cial constraints, a choice should be made 1 t- temporal and/or finan- 

"overview, study’and a^ 'tn ^ In™';,;! “r," "PP^fPrily superficial 
more important aspects of the Drohi.™ Th^ ' focused on one or two of the 
chapter is of the l^tL; typ^. investigation described in this 

tion i;d'rf‘i;g"h;-;i™te;?;n' »f = fipht test condi- 
tion features of two (static) pract 1 "y r''® acoustic propaga- 

Simulation facility are descrihe^ ? realizations of the ideal flight 

acoustic propagation through the f ree- iershea^^r^ ’ free- jet facility, 

able detail with the aid of aeomptriy Jsyer is studied in consider- 

with a semi -empi ri cal model of axial theory. This, in conjunction 

leads to a simple corr^Son pro^rr h'at 

acoustic data to estimated fMnht a converts free-jet measured 

work is described in Appendice; 3A threuglT;;, PsP«'s “f this 

on the so-caii;d^no^-“;;u';;icinK;;;;?oi;y"th;; ir';hr- r 

:if^rd;;;;;;;iro; ;;;;a“tr"“Vv 

out of the primary jet flow into the rad? t' ^ ' *■ Propagates 

description of thLe effects ar^ nr realistic 

Li 1 lay equation; the inte^preL^Sn of ^°«^tions to the 

forward motion effects. sLe numerical extended to include 

presented to illustrate how f low-acoustirfl/fir^ equation are 
of silence can be calculated to a onr>H • effects outside the cone 

analytic expression. The f low-acoustic^iS^e^a«ioreff completely 

-the" oo^pt; 

relative to the stat ion:^^^o;s^rveru;^r ?figt 
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3.1 


ACOUSTIC THEORY FOR FLIGHT SIMULATION FACILITIES 


3.1.1 The Ideal Flight (IF) Condition 

The Ideal Flight condition is shown in Figure 3.1. The acoustic sour ce- 
the aiJ^aft exhaust ~ moves with a uniform speed, Va, along a straight and 
eve! night path in an infinite atmosphere, i .e. the ground ° 

nections^and scattering from the ground surface are absent, 
aircraft itself is assumed to be acoustically transparent.) The speed 
aircraft = o of the atmosphere are uniform everywhere and 

sound and mean ^ 3 °’rest (apart from the motion due to the sound 

lav.'useU) L. atmosphere is also ideal in the sense that the sound wave 
is not attenuated except through spherical divergence. 

In Finure 3 1 the aircraft is shown in two positions; the first, in 

(or:L^’;or“;n':r"tS:jre:s;r:^i^‘h;stor:'oftL'’:" 

rLur'u::, :Lf"that putse bas travelled the dUtance Rr to t e oh- 

server at the speed of sound, aQ. At this reception time the observer/ 
aircraft angle and distance are denoted i|^ and R^. 

The relation between Gg and .1; can be derived by working in terms of the 
pulse travel time, t, where 


t = Rr/a< 


(3-1) 


as indicated in Figure 3-2. The result 


cot ijj = 


Va + aQCOsSe _ ^A ^°^^e 


3o sinSg 


is equivalent to that given by Ribner (ref. 30. ' -e- 


cosij; ■ Ma {1 ~ MaS i n^ip 


cos0„ - 


® {1 - MfiSin^ili}^ - MaCosiJ; 


(3-2) 


(3-3) 


cos6g = cosi); {1 


Mnsin2,(i}i - MflSin^iJ^ 


(3-0 


occurs when the reception angle is equal to 90 . 


The relation between emission and reception distances Rpi 
from the geometry sketched in Figure 3.2 


fol lows 
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AIRCRAFT POSITION AT 



79 



RECEPTION TIME, t = tg + t 

(t = Rr/ag) 



Figure 3.2b Ideal flight condition geometry at reception time 
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.3b Variation of emi ssion/reception radius ratio, 

with reception angle, , for several aircraft Mach 
numbers Ma“^A^^o* 




= {1 + cosOg} ^ 

or equivalently, according to Ribner (ref. 31 ) 


(3-5) 


- (1 - MA^in^i/;) V(1 +M^cos9g) . (3-6) 

This ratio is evaluated in Figure 3-3b for the range of parameter values 
used m Figure 3 . 3 a. 

In the next section this Ideal Flight condition is transformed to a 
frame of reference in which both the aircraft and observer are at rest. 


3 . 1.2 The Ideal Wind Tunnel (IWT) 

The transformation of the Ideal Flight condition, described in the pre- 
vious section, to a frame of reference in which both acoustic source and 
observer or microphone are at rest, is carried out in two steps. Starting 
with the IF condition at reception time, as shown in Figure 3.1|a, the first 
step is to set the microphone in motion so that it occupies the same 
position at reception time but is moving with the same speed and direction 
as t e aircraft, as shown in Figure 3-4b. The moving microphone measures 
the same acoustic pressure amplitudes as the stationary one but over a modi- 
fied time period; a pulse detected over a time period AIq by the stationary 
microphone is detected by the moving microphone over the modified time 
period At = Ato Dq where 

Do = {1 + coseg/ag} \ (3-7) 

Hence, in general, the microphones measure the same mean square pressures and 
the same proportional bandwidth mean square pressures — at different fre- 
quencies related by the Doppler factor Dq. That is, if the stationary 
microphone measures a pressure power spectral density (psd) P at frequency 
o)q and the moving microphone measures a pressure psd P at frequency w then 

‘^o-^o('^o) “ u).P((d) ( 3 - 8 ) 

where 


Wq - OlDg. (3-g) 

The second step is simply a Galilean coordinate transformation to the frame 
of reference that moves with the aircraft and microphone. This does not 
change the stationary-moving microphone measurement relations given above. 
However, it does mean that the atmospheric medium now moves with a uniform 
speed Vyy In this new reference frame, as shown in Figure 3-^c. This is 
referred to here as the "Ideal Wind Tunnel" (IWT): "Ideal" because the 

"cross section" is infinite and the flow velocity Is absolutely uniform and 
parallel (equal in magnitude to the uniform aircraft speed). Also, of 
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^2 


s 



(a) Ideal flight condition 
at reception time. 


(b) Ideal flight condition at 

reception time with co-moving 
microphone (velocity - V^) so 
that it occupies position 0 
at the reception time 
t = te + 


(c) Ideal flight condition at 

reception time in a co-moving 
reference frame such that both 
source and microphone are at 
rest; now atmospheric medium 
moves, with uniform velocity 
V/\ - "Ideal Wind Tunnel." 


Figure 3.^ The three stages in the transformation from the 
ideal flight condition to the ideal wind tunnel. 



course, the ideal properties of the atmosphere in th^ if 

over to this new reference frame. "" <=°"‘^^tion carry 

In the IWT the sound pulse emitted by the source at nncin c 
d ] ono the rav oafh ^ n tu j* * position So propa* 


3.1.3 


Fa^?t^'e^!"''nIt^r' Simulation 

i-acilities. Data Correction to IWT Conditions 


- tSFi - a 

the source. Threresenci of tha^ ^ a™ ^ ■'i^tance from 

f.cientiy ciose'oS^^^r^Lele,; th ca? "''“h".: If? ^f“-“ 

“fiietrrfii: if fifiiiit^ff-'ff fiir"';'-' 
ntfef;:i?"“ ^ochte^iaiffififMf? iziv^^tr 

that fa“;:i ffg“„"f?fiaif [ffeffif"" tef lect/scatter sound 

phone locations In the real wfnH^r ? direct radiation at micro- 

direct radiation Lerg^es trl^^^^ 'r facility, the 

as It propagates through the free-jet shear llajer°" Thf effects 

levels at the microphone locations'^nurQiHe. resulting radiation 

from the (required) direct radiation evicr* * ^ ^ differ considerably 
reflection In either typrof Jlc n^^ ^"^ 'sterna} 

fact. It will be argued^that usuallv^condi r°*^ thought to be Important, in 
tion within the free-jet shear layer can h 'a"* sound propaga- 

Then. by definition there a^rJ ^^ if acoustics, 
fraction processes can be calculaf h transmission and re- 

models. Turblllfla ffaff.ffi ?? ”' ''' "*" «“l>Iished anslyticsl 

limited expeffleffal'esif rfseff I fl'sflf • 

unimportant effect in the present cofJeiJ." ^^9gest that this is an 

on thf"Lf?mfm sff^cf-micropW disjlffl(f) '•estriction can arise 
of this distance is norfoffibTf tSs ill I ‘^^^^J^tent scaling 

of one hundred jet nozzle diameters Sr^ofrmne'f^ed 
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but only twenty diameters might be possible in a real wind funnel due to 
its finite cross-section. Since the jet mixing noise source is distr buted 
axially, reduction in the mi crophone-jet nozzle separation will result in 
a progressive distortion in the measured directivity pattern relative to 
that measured at large or true far-field distances. 

In what follows the emphasis is upon the free-jet facility and how 
theoretical and empirical models have been 

that convert measured data to estimated, true far-field IVfT data. The 
actual corrected or estimated IWT data is given in Section 2 along with the 
basic frc6“jet measured data. 

A detailed description of the models and correction procedures is given 
in Appendices 3A - 3F. In Appendix 3A theoretical expressions ’ 

based upon geometric acoustics, for the radiation from various types of 
source distributions (simple, quadrupole; volume acceleration, volume dis 
placement; at rest, convected; compact, noncompact) when the 
Ltion is immersed in a parallel, stratified, sheared flow. Appendix 3B 
deals with the more general case of wave acoustics propagation a" 
nite, parallel, stratified sheared flow, mainly in connection with the 
Li 1 ley equation solutions described in section 3-2. Thus, Appendices 3A and 
3B are not exclusively concerned with the correction procedure and "»ode1s, 
but they provide a general introduction and lay the foundation for the theo 
retical work in this section. 

The mean flow model used in Appendices 3A and 3B is stratified; that is, 
the mean flow properties vary only with the transverse coordinate. Thus, a 
ray may propagate a significant distance in the flow direction, but the 
axial lean velocity, for example, is not allied to vary in that d-rection 
as it does in a real jet flow. Appendix 3C describes an investigation in 
which that restriction is relaxed, but It is shown that axial variations 
and gradients have a negligible influence on the refraction of sound rays by 
the free-jet shear layer. That is, a 8tvat^f^ed f^oM model ^s adequate for 
present purposes and the ray paths may be drawn as if the shear layer has 
been replaced by a vortex sheet at the free-jet lip- line. On the 
hand, the theory of Appendices 3A - 3C rests entirely on the assumption that 
the sound propagation is governed by the laws of ge^tric acoustics (GA). 

A qualitative justification of that assumption is given in Appendix 3C; 
detailed quantitative information on the accuracy of GA for infinite strati- 
fied sheared flow models of the primary jet are being obtained from Lilley 
equation solutions under an on-going USAF/DOT contract. A complete 
Lilley equation study should be undertaken of radiation from quadrupole 
sources at realistic locations in a primary jet that is surrounded by a 
(finite cross section) free-jet, in order to assess the accuracy of the 
present gecxnetric acoustics approach. 

With the results of Appendix 3C the rest of the GA correction procedure 
is straightforward. In Appendix 3D, a simple but adequate model ?,ven 
for the effective axial location of the jet noise sources as a function of 
frequency or Strouhal number. In Appendix 3E the GA energy conservation law 
is used to relate measured acoustic pressures outside the free-jet with the 
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(required) acoustic pressures inside the uniform, 
the free-jet- The actual correction procedure is 


potential core region of 
then describee^ in Appendix 


3.2 THEORETICAL FLIGHT ALTERATION EFFECTS 
ON FLOW-ACOUSTIC INTERACTIONS 


thP ^ alteration effects on flow-acoustic interactions in 

’'""t can be calculated outside the cone of silence 

w?th '' APP«"dix 3A. For example, equation 

(3A-15) with equation (3A-26) can be written as » H 


Vo^^o) 


oc 


DoD 


2 

s , rel 


/D 


5 

m, re I 


if 


(volume acceleration source, 
denoted in subsequent analysis 
by V = 1 ) 


( 3 - 10 ) 




< COS0 


V(as + Vs^rel) 


(3-11) 


where 


Vs.rel “ Vg - Va 

®s,rel “ ^ “ l^s,rel cos0g/aQ 

•^m.rel “ " (^c " ^a) cos0g/ao}2 + {(v^^ 

* “s.ral 

t>o “ (1 + Va '^osOg/aQ}"^ 



(3-12) 

(3-13) 


(3-U) 


Is the mean ax I a 1 fluid veloci ty at the typical source position; that 

nr^cc^* of the IF (or IWT) proportional bandwidth mean square 

bv a volume acceleration quadrupole distribution is altered 

by forward motion in three ways, through: 


(a) the dynamic effect (Dq factor), 

(b) the eddy convection alteration effect (Vc^Vc ~ '^a) , and 

(c) the flow-acoustic alteration effect (Dg^-Ds ^el) 


31) un?- ' established features of Ribner's (ref. 

analoav !'T extension of the Lighthill jet mixing noise 

analogy model. It should be noted that the so-called dynamic effect can also 
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appear In the Ideal Wind Tunnel directivity expression where the observer is 
stationary with respect to the nozzle. (Here Vj has been replaced by y . 

In general, flow-acoustic alteration effects are described by numerical 
solutions to Li 1 ley's equation which, in effect, replace the geometric 
acoustics analytic result (given above, for example); essentially this means 
that the flow-acoustic interaction factor Dg rel replaced by the so- 
called flow factor . as outlined in Appendix 3B. Lilley s equation 

solutions are regarded as more realistic than GA results in describing flow- 
acoustic interactions although in many practical cases differences between 
the two are small. This is indeed the case, outside the cone of silence, 
when axial source locations are chosen according to the model described in 
Appendix 3D, and this will be illustrated with examples given below. 

The GA result given by equation (3-10) refers to a volume acceleration 
source distribution; the corresponding volume displacement result is 


“o^o 


(o)q) “ Dq 


(volume displacement 
source denoted by v = 3) 


( 3 - 15 ) 


In a recent paper (ref. 32) Tester and Morfey have presented arguments in 
favor of the displacement type of source, at least for modeling jet mixing 
noise with Lilley's equation. Thus numerical solutions given before are 
confined to this type of source [the flow factor has the superscript 
(v) =(3)] although throughout Appendices A and B and elsewhere in this 
section analytical results are given for both types of source since the 
volume acceleration source has the same frequency dependence as that of 
classical Lighthill quadrupole source distribution and therefore remains of 
considerable interest. 


To summarize lmc uli i • --i c / \ 

solutions to that equation define a value for the flow factor rp^lV^o^ 
F(u>), the ratio of the far-field radiation intensity to its value without 
the mean flow (e.g. Vg=VA or Vg = VA = 0) the source strength being held 
constant, 
context. 


the utilization of Lilley equation results, the numerical 


The flow factor Fpel^o) reference frame 


(\)) 2V + 3 

“ ''rel(“o) V^m,rel 


( 3 - 16 ) 


and F((jj) in the IWT reference frame contexts 

wP (o)) 




m 


( 3 - 17 ) 


Thus, although the proportional bandwidth mean square pressures are identi- 
cal. the modified Doppler factors representing convective amplification due 
to source motion are different in each case and the dynamic effect factor is 
excluded from the IWT expression; hence the flow factors are not i dent i cal . 
The IWT modified Doppler factor D^ = io^/u is given by equation (3A-19; or 
(3B-l) (in which the transverse wavenumber scaling velocities, 

Ve 2 =Ve 3 "'^et> neglected) 
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“m - (“-'>1 «cl)“ + k? = ,„20^ ' (3-, 8) 

The wavenumber kj is first defined in the IF frame 

•^1 = cos0g (3-19) 

and IS invariant under the Galilean transformation to the IWT frame. Thus 
can be written as 

Dm - {(l - Dq V(-2 cos0g/aQ)2 + Dq (VgiCos0g/aQ) ^ (3-20) 

The appropriate modified Doppler factor for the IF f rame, rel » ' s simply 
^m.rel ^m'^^o 

[{1 - (Vgj -Vq) cos0g/aQ}2 + {Vgj cos0g/ajj}2]^ (3-21) 

where now the eddy convection velocity appears in a "relative to observer 

form [as it does in the Ribner (ref. 31 ) or Ffowcs Williams (ref. 7) 
resu I ts J . \ / 

In the following section the relation between the two different flow 
factors IS given and the choice of reference frame (in which the Li 1 lev 
equation has been solved) is discussed. ^ 

3 . 2.1 Choice of Reference Frame for Li 1 ley Equation Solutions 

the 3B the relation between the Li 1 ley equation flow factor in 

the ideal flight and wind tunnel reference frames, Fpel, F is given as 


^ ^rel^‘*’o) (volume displacement source) 


F<‘'U) = 0| 


(volume acceleration source) 


( 3 - 22 ) 


( 3 - 23 ) 


where 


Dq “ = 0 + V/^cos0g/aQ)“l. 

of t'lle modified by the presence 

of the co-flowing stream (see Figure 3.5) of velocity V = V*. In the IF 

reference frame, Frei(o)o) changes in response to the reduction in the mean 

ablve'^o^ T ''j Vj -V^ = Vj,rel. With the simple relation given 

a^H fl^ f reference frame to the other is straightforward 

the [p f can be generated in either reference frame. Here 

the IF frame is chosen since the flow profile is identical to that in the 
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'^a ~ '^A “ '^a,rel 


V 


a 


(a) Ideal wind tunnel . 


(b) Ideal flight. 


igure 3-5 Typical mean flow profiles for the Li 1 ley equation 
in each reference frame: (a) ideal wind tunnel, 

(b) i dea 1 flight. 



(ref. 33), which also incUra ?c^^ JelutEra^ 

comparisons with measured static jet noise di recti. i EydaEE 


3*2.2 Definition of Paramete 


rs 


However, the results given in referf^nr^>< (oo\ 

H!^tiit'HolL?£;i“” 


The frequency parameter is the 


modified Strouhal number where 


^m ~ ^ ^m , re I 


^ ~ ^o'^/'/j.rel 


=“r,/27r) 


(3-24) 

(3-25) 


•^m.rel = iol,re] t ’ + “? (Va , rel /^o) " ] + (af-a^) ^el ‘^osGe/a^) 2 (3-26) 


^s.rel - {1 + (Vg-V/^) cos0g/aQ}“^ 

ai = 0.2, 01 = 0.4 
(Vs -Va) = 0.6 (V 3 -Va). 


( 3 - 27 ) 
(3-28) 
( 3 - 29 ) 

(3l-2()°(?rwld™ ?r from equation 

velo^tieE V J ^ "““"■Pfio" fot the wavenumber scalEnV 

ocities Vgi, Vgj are proportional to V, „i; the values of m ft . 
obtained by Szewczyk and Morfev (ref 'fri^ cne values of aj. 0i were 

?hTi;Eo de-:E:i^^:tr“t 3^. 

d-::u'E);aEE- - ;-tf-vorEhf niq7E“Ea'se'w”?:;s^^ :;rqEa^nti- 

place of Dn, .el in equ^t iori 3 - 24 r^"^hr.!"i"l!!f^^ could be used in 


lace of D,,,el in equation 3-24r T re:act va u?^^^^ 

3eris"“EeEd"E:„“:aEE°:L^Ee^ -- - -'e 


prrmcLTy 


COSGa 


COS0, 


■’s ,rel 


(3-30) 
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In the following numerical results, error function veloc. ty prof i les 
have Len used which correspond, nominally, to actual profiles at the 
locations of% and 1. diameters in the intial mixing region of an ax i symmetric 
iet In addition, an "inverse exponential squared" profile has been used to 
represent the transition and fully developed regions. Each profile is de 
fined by the ratio of its vorticity thickness, 
at which the mean velocity is 66.3'<:' "f ' ts center 


and the radius, 
[ i ne value. 




3 . 2.3 Numerical Results 

At the nominal one diameter axial location (■Sjri = .2e), Li 1 ley 
flow fLtor remits a,e shown in Figure 3-6 for - 60 “ and mod. f ,ed Strouhal 
numbers S = 1 0 2 0 and ^i.O. According to the simple law described in 

ApSndix 3D. th^ one diameter location is equal to the effective 

i =q radiation. Thus, radiation at Strouhal numbers less than this 

tend to orTginate further downstream. The results in ^ ‘ ^ 

clearly that radiation at Strouhal numbers greater than S„, 
geometric acoustics limit by less than 1/2 dB. 

At the nominal four diameter axial location for 6s = 60° and S^, = 0.5 , 1 -0 
and 2 0 the flow factor results shown in Figure 3-7 obey a similar trend, 
^is axial location corresponds to the effective one for S„ a 1 .0 rad , at ,on 

"d^ ^'r::rof“!nte“^tnrt^:ie %::n;“s'tha: .he ra.e 

:xiar:ou:ce"i:crtton ra":uo'n'wrtiltrot,har:;m6e? is 

r-' " ^hire',: ^i^tljripe^d^ r?: In Ihe-' 

ui^ey equation flow factor, fU,' (« o) "6ich can be closely predicted by the 
change in the GA limit. That limit is given by 


:( 3 ) 

rel 


(ojq) 


GA 


- O' 


tji 

s , rel 


(3-31) 


its derivation is briefly described in reference (32). Th.s feature ,s also 
e^hi^tled b)the results shown in Figu.e 3-8 for a ful ly-developed mean 
velocity profile. In this case, a different Sl.oohal number is util z , 

Sm', where 


Sm' - (Dm, rel Lq) (^ '‘l)/'^a,rel 
= (2 rj/d) (Vj , rel , rel ^ • 


(3-32) 


Since 2 r^ is approximately equal to (V, for this profile, is roughly the 

local modified Strouhal number. 

(3-33) 


(^m,rel ^o^ , re 1 
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2.0 
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LILLEY 

EQUATION 

SOLUTIONS 



i,rel '^a " '^A 


''a.rel/ao 


Variation of the flow-acoustic interactions factor, , with 
the relative centerline velocity, Vg^rel. <"rom Li 1 ley''lquat ion 
solutions and from geometric acoustics theory. Profile: One 

diameter axial location, = 0.26; V-iVi, Source emission 

angle, 9 = 60°. ® 




, re 1 ^ 


(3) 

Figure 3-7 Variation of the flow-acoustic interactions factor, , with 

the relative centerline velocity, from Lilley equation 

solutions and geometric acoustics theory. Pvofite: Four ^ 

diameter axial location^ " 1-05; V 3 - Vj. Source emission 

anglcy 0 ^ = 60^. 


9 ^ 






2 


T» -i, 


O 0.25 
□ 0.50 
A 1-00 


Figure 3*8 Variation of the flow-acoustic interactions factor, f Lj , with 
the relative centerline velocity, Vg^rel. from Li 1 ley'^^uat ion 
solutions and geometric acoustics theory. Profile: Fully 

developed^ 'S^^/r = 1.82; Va?^Vj. Souroe emission angley Qo = 6 




i.e, one based on the local shear layer thickness and the local centerline 
mean velocity. This takes a value of order unity so that the difference be- 
tween the Li 1 ley equation flow factors and that given by the GA limit for 
all values less than unity are small and the slopes are nearly identical. 

The main influence of the wavenormal angle is to modify the differ- 
ences between Li 1 ley equation and GA flow factors as illustrated in Figure 
3.9 for 63 = 30 ^; the slopes, however, remain very similar. The variation of 
00 with 63 is given in Table 3-1* 


(''a,rel/ao)/6s° 

0 

30 

60 

120 

.25 

29.6 

40.0 

62.3 

122.7 

.50 

39.7 

46.6 

64.2 

126.0 

.75 

46.4 

51.4 

65.9 

130.2 

1.00 

51.3 

55.2 

67.4 

135.6 

1.25 

55.2 

58.3 

68.7 

143.1 


Table 3*1 Emission angle values, as a function of 
source emission or wavenormal angle 63 and 

ya,rel^^o ^s,rel 

(cos 0 g = COS 05/(1 coseg/ao)). 


3 . 2,4 The Geometric Acoustics Result for Data Correlation 
Outside the Cone of Silence 

It is clear from the preceding results (including Table 3.1) that for 
moderate changes in V 3 or A, at constant Gg or G^, the GA limit can be 
used, to a good approximation, to calculate estimated forward motion altera- 
tion effects on flow-acoustic interactions. From equation (3B-12) (Appendix 
3B) and equation (3*31) it follows that changes in directivity at constant 
Strouhal number in the IWT reference frame can be estimated from the 
expression 


D 


6 

s , rel 


Do/D 


9 

m, rel 


(3-34) 


where 


Ds,re1 = 1 - (Vs-Va) coseg/ao 
= (1 + Va coseg/ao)”^ 

and D^ Pgj can be evaluated with equation (3-26). An alternative form for 
rel’ Section 2 for data correlation purposes, is 

^m,rel ~ (*^s,re1 ^'^c, rel ^ (3~35) 

(a = 0.3) . 
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There the Doppler factor I -s c a i ( 1 a led with 


Equation (3"36) based uprn? i 
velocity, and the nean flow c h- : 

Vg , are identical. Also, it imi ■ 
changes with forward r!r?iiofi s!iK.<. ii 


V 


J , re 1 


(3-36) 


ossuinption that the eddy convection 
j! thv rftective source location, 

iia: the .dry -onvection velocity 

i i ^ ; -n as 


'^J , -- 0 . 6 / V , i . 


(3-37) 


Neither equation (3^36) nor ecjuatoM 1 / i /mc boscu on any solid evidence, 

experimental or theoretical; they oj - j ' e^iir^pnons which are conven-* 
ient to use and appe<jf o) r^e w Om:, 

The expression (3 ‘3M above h;W) i i , and the results compared 

with overall direclivit> aiteratio!! i^: ; :iear iwec data as shown in Figure 
2.2. The excellent ay tv-cnsei 1 1 - -ieqraiv: in-efittve to extend the 

present i nves t i ya L t on tr> t hr- low ^ aiyl. w nside the cone of silence 

and to non- i sothe t ilia ! jet^v. 


CONhi USlONb 

Detailed conclusions lor a.r. Ii vd Miv : la i n - s t i ga t i on areas have been 

given in the appropriate subsurtions aoi , iLo)vnr! i t : they can be summarized 
as foil ows : 

(1) It is necessary to dt:tou^ * d,^a 1 i /tol, cvnfn[)tual models of flight 

and simulated fliglit cc^nditioris ir '>r !ef i ) hlgidlght the important acoustic 
features and i nter re 1 a t i onsh i ps sjkI spto i fy i [jo basic acoustic features 
required of flight s i mu 1 a t i on fucil Mts In r f:: ‘real Flight (IF) condition 

emission and reception angles and distances inu-d re distinguished; here their 
functional relations have also bem def oroJ an i evaluated. In the Ideal Wind 
Tunnel (IWT), IF data can be acqrjired i* the distance travelled by the sound 

wave in the wavenormal direct inn eiid [he ..aven o ;’:i nigle are both held 

constant. Then only a Doppler i r egnen: ■ hiiL ine dr r o be applied to convert 
data from the IWT to the IF condition. Aitiplitede ' c or rec t i ons“ are not re- 
quired for the mean sgiiare i>rvrss'Me n -. ■<>r I ■ * e|H>r t i ona 1 bandwidth mean 

square pressure. 

(2) In a real wind tunnel Miuh? e -o ility, unvjanted reflec- 

tions and scattering from the oucei huuMuiary r cn] i ;.n r,f the uniform tunnel 
flow are probably unimportant- Hewever , (l>e area ratio effect and the 

axial distribution of the tuirbulen! jet raxing rn/tst source should be inves- 
tigated. The finite ar'ea ratio effect drould be ixaMoned with the aid of 

Li 1 ley equation solutions for the coirspK In firs, field. The axial source 
distribution effect can be niodeleii eiiip i i ' u 1 ; c n. v cid distance and angle 
corrections for the f’leasured data. Otlit-i wise i nt ideal and real wind tunnel 
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are acoustically equivalent if ’’atmospheric’' attenuation effects in the real 
wind tunnel can be estimated accurately. 

(3) Acoustic propagation in the flight simulation facility 

has been studied in cons i de^rab 1 e detail in tlie geometric acoustics approxima- 
tion. The approximation can be Justified, / be ; ^ , for the present 

application by an analysis of the Li 1 ley 6^,c^Uf^Lion in relation to the 
equation governing geometric acoustirs opaga f i on . However, the accuracy of 

Geometric Acoustics (GA) should be assessed in future work with a quantita- 
tive study based upon Li I ley equation solutions. 

The position and angle of ray path^ tfirouqli tlie free-Jet shear layer 
have been calculated; the mean flow vjithin the shear layer was allowed to 
vary axially in a realistic way. The results are nearly identical to those 
that are calculated with a vortex sheet model of the free*-jet flow, over a 

wide range of parameters. However, the amplitude of the sound along a ray 

path is governed by the ray tube gec>metry mhl tlie law of energy conservation 
in a ray tube, not the transmission of a sound ray by a vortex sheet. Hence, 
angle and distance corrections that convert free-Jet data to estimated Ideal 
Wind Tunnel data can be based on a standard vortex sheet model while trans- 
mission-amplitude corrections sliould lie calculated with the analytic formula 
derived here. 

The complete correction procedure also uses an empirical model of the 
turbulent mixing noise source axial distribution (as a function of Strouhal 
number) to calculate additional angle arrd distarKse cor r ec t i ons , 

(k) The influence of forward motion on f 1 ow~’acc^us 1 1 c interactions 
(within the primary Jet flow alone) lias beer^ examined in the Ideal Flight 
condition through a re- i nterpretat ion of the "static” numerical solutions to 
the Li 1 ley equation and wltii geometric acoustics. It has been demonstrated 
that outside the cone of silence the I illey equation radiation levels change 

with forward motion in a way that can be calculated, to a good approximation, 

with the analytic, GA result. When this GA f 1 ow^accMjs t i c alteration effect 
was combined with the standard, tlieoret i c a I , dynamic and convective amplifi- 
cation alteration effects, the predicted 'llrcctivity alteration with forward 
motion was found to be in excellent agreement witii measured free-jet data- 
outs I de the cone of silence. This successful correlation of directivity 
data Justifies an extension of the present wot k to include forward arc 
angles, angles inside the cone of silence and non i sot lie rma I primary jet flow. 


99 



4 . 


EFFECT OF FORWARD VELOCITY ON THE 
STRUCTURE OF A TURBULENT JET 


4.1 INTRODUCTION 


In this chapter of the report, the changes in structure of a round jet 
caused by variation of free stream velocity will be examined. The work to 
be described is in three sections. The first section is an experimental 
investigation, and turbulence measurements of a round jet in a wind tunnel 
are presented. These measurements are used to provide scaling laws for 
various characteristic lengths and velocities as functions of velocity 
ratio and jet Mach number. The second section is a numerical study of the 
development of a turbulent jet in a moving stream. The results of this 
analysis are compared with experimental measurements. The third section 
makes use of the scaling laws derived from the experiments to predict the 
changes in radiated noise with changes in velocity ratio brought about by 
alteration of the turbulent jet structure. Several appendices associated 
with this chapter are also included. These provide a complete tabulation 
of the velocity measurements and a discussion of the data processing tech- 
niques used to analyze the laser velocimeter measurements. A simple 
analysis is also available to describe the decay of the jet center- 
line velocity. Finally, an original method, based on the velocity at 
the dividing streamline, of noise radiation sea ling with changes in f orwa r d 
velocity is given. 

Before the investigations are described in detail, it is necessary to 
discuss to what extent the experiments performed simulate the case of a 
full-scale jet moving through stationary air at a uniform forward velocity. 
There are two questions which need to be considered. Firstly, whether the 
stationary jet in a moving stream simulates a moving jet in a stationary 
stream from the point of view of turbulent structure? Secondly, can the 
measurements made at model-size be used to predict full-size structure? 

The answer to the first question is based simply on a change in 
reference system. A control volume is considered which surrounds the dis- 
turbed flow region around the jet and the axes are fixed relative to the 
jet exit. So long as the boundary conditions in the jet reference frame 
are unchanged, the only difference that can exist is due to changes in the 
equations which govern the fluid motion due to motion of the reference 
frame. So long as the motion of this reference frame is given by a uniform 
linear velocity, the forces on a fluid particle remain unchanged. Clearly, 
this would not be true if the frame of reference fixed with the jet were 
moving with some angular velocity or a uniform acceleration. Thus, uniform 
translation of the axes fixed with respect to the jet will not change the 
fluid motion in that reference frame. It was noted, however, that the 
boundary conditions must not be changed by motion of this reference frame. 
This leads to two points. Firstly, the flow is not simulated by a jet into 
stationary air whose jet exit velocity is equal to the difference in 
velocity between the original jet and free stream since this transformation 
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has not moved the jet nozzle itself with the appropriate uniform motion of 
the reference axes. Secondly, the jet development for a jet moving In an 
infinite quiescent medium is only simulated by a stationary jet in a moving 
medium if the moving medium is also infinite in extent. This condition was 
referred to as the “ideal wind tunnel*' in section 3- Clearly the real wind 
tunnel will not satisfy this criterion completely due to the presence of 
the tunnel walls. This change in boundary conditions manifests itself in a 
deceleration of the medium and hence the creation of an adverse pressure 
gradient. However, the degree to which the real wind tunnel changes the 
flow is a function of the relative dimensions of the tunnel cross-section 
and the jet. The magnitude of the free-stream pressure gradient will be 
considered in the experimental study given in section 4.2. It should also 
be noted that a reference frame fixed with respect to the jet is the obvious 
choice for a description of the fluid motion since, once the operating con- 
ditions have been set, the motion is stationary in this frame. Clearly, the 
motion is not stationary for a frame of reference moving with the surround- 
ing air. 

The second question to be discussed is whether the model size experi- 
ments may be scaled to full-size? It has already been noted that the 
equations for the fluid motion are unaltered by the choice of reference 
frame, subject to the conditions given above. Thus, in an ideal case the 
two situations may be scaled with the appropriate use of nondimensional 
quantities such as the Reynolds number and Mach number, etc. However, the 
situation is not ideal and a number of factors may influence the scaling. 

If the question is confined to scaling of a cold model scale air jet to 
full scale, the changes in the fluid composition or nozzle geometry which 
would be encountered in a real jet engine do not apply. However, certain 
conditions are likely to not scale. The wind tunnel free-stream turbulence 
has a particular amplitude and length scale. This turbulence will influence 
the jet development tending to enhance the mixing process if the amplitude 
is large and the scales' small. If the atmosphere through which the aircraft 
is assumed to be flying were completely at rest and the aircraft structure 
did not create any turbulence in the flow surrounding the jet, there would 
be an obvious difference between the wind tunnel simulation and the flight 
case. However, the atmosphere Is not still and the aircraft surfaces may 
influence the air surrounding the jet, though relative levels and scales of 
the atmospheric turbulence are unlikely to be duplicated in the wind tunnel. 
Another feature which is unlikely to scale exactly is the turbulence at the 
jet exit, being a function of the flow configuration upstream of the jet 
exit. The Initial turbulence level is known to influence the jet 
development. 

It can be seen that those factors which prevent identical scaling are 
in each case due to non-scaling of the boundary conditions. The degree to 
which these changes may influence the flow in this situation needs examina- 
tion; however, in gross terms their influence is likely to be overshadowed 
by the alterations expected when a real jet engine is considered. It is to 
be concluded then that the turbulent structure of a full-scale jet with 
forward velocity moving through stationary air may be well simulated by a 
wind tunnel model jet experiment. 
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In the next section an experimental study of the effect of a secondary 
stream on the development of a turbulent jet will be described* 


4.2 EXPERIMENTAL INVESTIGATION OF THE EFFECT OF A SECONDARY 
STREAM ON THE STRUCTURE OF A TURBULENT JET 


In the following section an experimental investigation will be described. 
Measurements of turbulent velocities were made, using a laser velocimeter. In 
a turbulent jet mounted in a low speed wind tunnel. The variation of the jet 
development with the change in velocity ratio and jet exit Mach number were 
examined. Before describing the experimental study, previous investigations 
of two-stream mixing problems will be discussed. 


4.2.1 Summary of Previous Work 

There have been few experimental investigations of the mixing of two 
streams of finite velocity. This has been mostly due to the difficulties In 
constructing adequate experimental facilities. Those experiments which have 
been conducted can be divided into five classes. These are: 


(i) two-dimensional mixing layer, 

(ii) the two-dimensional or plane jet in a moving stream, 

(ill) the round jet in a wind tunnel, 

(Iv) the round jet in a free jet facility, and 
(v) coaxial round jets. 

4. 2. 1.1 The two-dimensional mixing layer. The two-dimensional mixing 
layer is perhaps the simplest two-stream flow to be examined. Analytically, 
convenient similarity solutions for the mean velocity have been obtained 
using an eddy viscosity assumption and this has prompted several experlmert" 
ta 1 investigations. Sabin (ref. 37) used a water channel with movable side 
walls to examine the effect of pressure gradient on the development of the 
mixing layer. Velocity ratios between .47 and .66 were obtained for very 
low flow velocities. Measurements of mean velocity were made and it was 
concluded that the dimensionless velocity profiles in a similarity coordi- 
nate could be expressed as functions of the velocity ratio alone. Yule (ref. 
38 ) made an extensive study of the mixing layer. Two velocity ratios, .3 and 
.61 were examined. Measurements were made of mean axial velocity, three 
components of the velocity fluctuations, shear stress and longitudinal and 
lateral space correlations. The variation of the measured spreading parame’* 
ter agreed with the measurements of Miles and Shih (ref. 39) and was given by 


2 - = (1 - 
^o (1 +A)i ’ 


(4-1) 


where o is the spreading parameter defined by Gortler (ref. 69) and Oq is the 
value of o for A, the velocity ratio between the two streams, equal to zero. 
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It should be noted that the data of Sabin (ref. 37) gave a relationship of 
the form, 


£_ = (1 - A) 

0Q ( 1 + a) 


(4-2) 


Yule (ref. 38) also noted that the peak value of the nond i mens iona 1 i zed 
total turbulence intensity, ^^/(Vg-Vb)^, increased with A. Jones et dl 
(ref. 40) also made spatial cross-correlations and calculated the variation 
of convection velocity across the layer for a primary stream velocity of 
30.5 m/s and a secondary velocity of 9.14 m/s. Measurements of intermit- 
tency were also made. Brown and Roshko (ref. Al) made measurements using 
air, fiitrogen, and helium. Use of these various gases allowed the effects 
of velocity ratio and density ratio to be examined separately. It was con- 
cluded that a relationship of the form of equation (^-2) could be used to 
describe the rate of growth of the shear layer. I t wi 1 1 be shown in section 
that the behavior of this spreading parameter for small values of X is 
important to the prediction of noise radiation. Brown and Roshko (ref. ^1) 
gave a very careful discussion of possible descriptions of this behavior; 
however, the scatter in the data did not enable one relationship to be pre- 
ferred over any other. Brown and Roshko (ref. ^1) also showed that the 
change in the growth of the mixing layer due to density differences in two 
streams of different gases was less than the change due to density differ- 
ences caused by high speed flow. This latter effect was interpreted as 
being due to compressibility rather than density difference. 


^.2.1.2 The plane jet. The plane jet exhausting into a moving free- 
stream was examined experimentally by Bradbury (ref, 42), It was shown that 
departures from self-preservation for small velocity ratios of .07 and .16 
were of no significance. Bradbury and Riley (ref. 43) considered a. wider 
range of velocity ratios and concluded that the flow changed from a pure jet 
flow near the nozzle exit to a self-preserving wake flow far downstream. 
Weinstein (ref. 44) also examined the plane jet and made mean axial velocity 
measurements for velocity ratios between .5 and 2. 


4. 2. 1.3 The round jet in a wind tunnel. Curtet and Ricou (ref. 45) 

and Antonia and Bilger (ref. 46) made measurements of the structure of a 

round jet in a wind tunnel. The former investigation was principally made 

to study the effect of the ratio of wind tunnel radius to jet radius on the 

free stream and jet development. Clearly measurements of this type do not 
simulate the 'ideal wind tunnel' with an infinite moving stream and the 
results cannot be used for prediction of forward flight effects. However, 
the results were used for an experimental comparison by Hill (ref. 47 ) whose 
analytical work enables the influence of the tunnel walls to be estimated. 
Antonia and Bilger's (ref. 46) experiments were for low velocity (28 and 
42 m/s) round jets at velocity ratios of ,22 and .33. The investiga- 
tion concentrated on whether assumptions of turbulence similarity could be 
made far downstream of the jet exit where the flow approached a small 
deficit wake. They concluded that the flow far downstream depends strongly 
on the complete past history of the flow and that no turbulence similarity 
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assumption was possible. Forstail and Shapiro (ref. 2) also made measure** 
ments of a round jet in a wind tunnel. Their results will be of consider- 
able interest in this current investigation. Measurements of mean velocity 
and mass concentration (helium was used in a tracer in the primary jet) 
were made for low speeds with velocity ratios in the range .2 to .75. On 
the basis of these measurements the following empirical formulae were 
proposed. 


x^. = 8 + 24 X, 
^a “ X _ 

~T^~T 

f.5 = (S/Xc)^’^, 


(4-3) 

('»-'*) 

(4-5) 


Vl- ^ I 

^a - A ■ 2 

Landis and Shapiro (ref. 48) and Pabst (ref. 49) considered the development 
of a heated air jet exhausting into a moving stream. The measurements by 
Landis and Shapiro (ref. 48) were of the mean axial velocity and were made 
for velocity ratios of .333 to .852. Pabst's (ref. 49) measurements are 
reported in detail by Szablewski (ref. 50). The primary jet velocity In 
these measurements was 400 m/s and the jet temperature was 300°C . 

Velocity ratios in the range .045 to ,47 were examined. The measurements 
of mean axia] velocity and temperature showed good collapse radially when 
plotted as (V^ “A)/(Vq“X) and T/f^ against f/x, downstream of the potential 
core. Von Glahn et al (ref. 51) made pitot-static surveys for various 
nozzle configurations, including a convergent circular nozzle, in a wind 
tunnel. Primary jet Mach numbers between .45 and 1.02 were used and Mach 
number ratios of between .1 and .75 were examined. It was concluded that 
the centerline velocjty could be correlated for the various Mach number 
ratios by plotting (V3-A)/(1 - A) as a function of [x/(2Cp/ 1 + Mj) where 
Cp, was the effective nozzle coefficient and b was an empirically developed 
exponent given by, 

b = 1 + y [X2 - 1]■^ (4-7) 


1 + COS 


Trr 


2F 


.5 


(4-6) 


4. 2. 1.4 Coaxial Jets, Alpinieri (ref. 52) made measurements of 
axial mean velocity and concentration for coaxial jets. Hydrogen and 
carbon dioxide were used as the inner jet gases and the outer jet was of 
air at 198 m/s. Velocity ratios in the range .8 to 2.13 were considered. 
The hydrogen temperature was 270 ±5*5^C and the carbon dioxide temperature 
was 254 ±5*5°C Coaxial jet flows have also been studied by Peters et al 
(ref. 53) using hydrogen and air, Fricke and Schorr (ref. 54) using helium 
and nitrogen and Zawacki and Weinstein (ref. 55) using freon and air. 
Velocity ratios in these experiments ranged from .12 to 39*5. in each case 
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axial mean velocity and concentration were measured. Zawacki and Weinstein 
(ref. 55) also measured the axial turbulence intensity and the mean density. 
Eggers and Torrence (ref. 56) made pitot static traverses in a coaxial jet 
with the outer Jet exit velocity always greater than the primary jet 
velocity. Measurements were made for supersonic outer jet velocities and 
slightly subsonic inner jet velocities. 

^. 2 . 1.5 Summary. Several points emerge from this review of previous 
experimental studies- Because of the relative ease of measurement, the 
mean axial velocity is the most often measured quantity. Measurements of 
turbulence intensity have been made in relatively few experiments and these 
have been limited, by available i nstrumentat ion , to low velocities. The 
main interest in round jets in moving streams appears to have been concen- 
trated in two distinct areas. Firstly, the rate of the decay of the jet 

centerline axial mean velocity, has been studied. Secondly, the assumptions 
of similarity and self-preservation have been examined. This has, for the 
post part, involved measurements downstream of the end of the potential core. 
Clearly, the following areas are in need of attention: 

(1) the effect of the secondary stream on the initial development of a 
round jet, 

( 2 ) the variation of turbulence intensity with velocity ratio, and 

(3) the variation of turbulence intensity with jet Mach number. 

The progress made toward understanding in these areas will be described 
in the subsequent sections. 

4.2.2 Description of Experimental Facilities 

In this section the facilities used for the measurement of the effect 
of velocity ratio and Mach number on the flow properties of an axisymmetric 
jet are described. 

The wind tunnel used in the experiments, shown in Figure 4-1, was a 
closed-circuit, s i ng le- return , low-speed facility located in the Research 
Laboratory. The test section was approximately .76 m x 1 . 09 m with a length 
of 1.22 meters. The flow was driven by a 1.83 m diameter, 10-bladed, single- 
stage axial flow fan. The fan was driven by a 3x10^ Watts, 1200 rpm, 
synchronous speed induction motor. The fan speed could be varied over the 
range from 0 to II 50 rpm by means of an eddy-current, coup 1 i ng-type variable 
speed unit. Steady state speed regulation was 1/8 of 1 percent. A standard 
friction brake rated at 1017 Joule braking torque could bring the fan to 
rest from full speed in less than 15 seconds. The tunnel provided an empty 
flow velocity capability of 0 to 94.5 m/s. 

A nozzle installation is illustrated in the photograph in 
Figure 4 2. An air duct 20.3 cms, in diameter was installed to supply air 
through a 5 .O 8 cm. automatically controlled value to the 5 .O 8 cm. diameter 
nozzles which exhausted axially down the center of the wind tunnel working 
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Figure h.2 Pylon, plenum and nozzle installation. 
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section. As a result of the proximity of the tunnel screens (shown cross 
hatched in the schematic) to the metal contraction section just upstream of 
the working section, access for the air duct and pylon mount system was 
extremely limited. A rectangular duct was designed with a cross-sectional 
area equivalent to a 15^2 cm. round duct. This rectangular duct, made from 
.95 cm. steel plate, formed the basis of the pylon mount which was fitted 
with wooden leading and trailing edge fairings and contoured with aluminum 
fill. The airfoil section was a 21Z thick zero lift design. A 20.3 cm. 
diameter nacelle made from standard pipe was mounted on top of the pylon 
with a wooden bullet upstream fairing as shown. A 5-08 cm. thick section 
of .32 cm. cell honeycomb was installed in the 20-3 cm. plenum as a flow 
straightener. A 20.3 cm. to 10.2 cm. reduction was attached to this plenum 
and smoothly f a i red to the end of a short length of 10.2 cm. pipe to which 
the model Jet nozzles were attached. This air supply ducting was aligned 
axially in the tunnel by use of a low power laser mounted at the downstream 
end of the working section. The nozzle exit plane position was governed by 
the extent of the optical glass wall through which the laser beams were 
passed. Because of the limited extent of this glass panel a 22.9 cm. long, 
10.2 cm. I/D spacer pipe was made to be inserted Just upstream ot the 
nozzle so that a total of 17 diameters of the jet flow could be explored by 
removal or insertion of the extension. The laser velocimeter, shown 
schematically in Figure 4.3 and described in detail below, was mounted on a 
hydraulic table outside the wind tunnel working section providing vertical 
positioning. Axial positioning was achieved with a lead-screw and way 
system, similar to that of a center lathe, aligned parallel to the tunnel 
centerline. Lateral traversing capability was provided in a similar 
manner. 

Slots in the tunnel wall Just downstream of the working section were 
provided to vent the working section to ambient. Since a considerable amount 
of extra air was to be Injected into the tunnel via the Jet nozzles, some 
means of [)reventing a pressure build-up was needed. This was achieved by 
providing a “bleeder’' door Just upstream of the fan as shown in Figure 4.1 . 
The opening or closing of this door was achieved by means of a cable and 
winch system controlled from the test stand. 

The measurement of the mean properties of a turbulent flow has a practi- 
cal importance, but it is only by measuring the fluctuating components of the 
flow that it is possible to obtain a better understand i ng of the mechanism of 
the turbulent flow and provide Justification for any mathematical hypothesis 
used in representing the motion. Measurements of the time averaged pro- 
perties of the flow may be achieved by total and static pressure surveys in 
high velocity and temperature flows. However, measurements of turbulent 
fluctuations have normally been limited to fairly low speed flows since in 
these velocity regions the conventional measuring devices, the hot-wire or 
film anemometer give uncontaminated velocity measurements and are structurally 
sound. It is only with the development of optical measurement techniques 
that accurate measurements have been viable in hot or high speed turbulent 
flows. There are two main advantages in using optical techniques. Firstly, 
there is no need to insert a probe in the flow which might disturb the flow 
regime or be physically incapable of withstanding stresses exerted by the 
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Figure 4.3 Plan view of LV optics platform. 
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flow. Secondly, an unambiguous measurement of velocity is obtained and 
measurements can be carried out in regions where the density and tempera- 
ture may fluctuate. 

There are a number of optical techniques for remotely detecting the 
velocity of particles and these methods are well documented (ref. 67)- 
approach to be described here, which has been developed by Lockheed-Georg i a 
Company, is based on an interference pattern of light formed in the measure- 
ment volume by the intersection of two coherent monochromatic light beams. 

As a microscopic particle passes through this fringe pattern, light is 
scattered and detected by a photo-sensor. The detector output signal burst 
has a frequency depending on the spacing of the interference pattern 
and the velocity component of the particle normal to the fringes. Since the 
fringe spacing is set by the geometry of the optics, the normal particle 
velocity is readily derived from the detector signal frequency. 

A detailed account of much of the optics, electronics and data process- 
ing used in the current laser velocimeter system is given in references 58 
and 59- However, a brief description of the configuration used in these 
experiments and some modifications that have been made will be described 
here. 

A plan view of the optical configuration is sketched in Figure 4.3 and 
a photograph of the complete unit is shown in Figure 4.4 . A 4-watt argon 
laser is mounted beneath the optics platform. The beam splitter/color filter 
assembly separates the two predominant lines of the argon laser into four 
beams; two blue and two green. These beams are turned through a dove prism, 
which allows rotation of the fringe patterns, and are then simultaneously 
focussed and caused to cross by the transmitting lens. This creates two 
coincident ellipsoidal measurement volumes. Light scattered by the particles 
passing through the measurement volume is collected by the receiving optics, 
mounted beside the transmitting optics, and is fed to a photomultiplier 
assembly. This assembly filters the collected back-scattered light into blue 
and green components which are sensed by two photomultiplier tubes. The 
resulting signals are subsequently processed to provide the particle compo- 
nent veloci ties. 

In this configuration the laser velocimeter is able to provide measure- 
ments of two orthogonal components of turbulent velocity as long as the 
deviation of the instantaneous velocity vector from the axial direction is 
sufficiently small. Since the angular criterion is not met by the expected 
velocities in a turbulent jet, these fringe patterns were rotated 45° for the 
current measurements and only the set of fringes normal to the jet axis were 
used. This meant that measurements of axial velocity only were made. The 
angular restriction can be alleviated by movement of the fringes in the 
measurement volume. This can be achieved by use of a Bragg cell in the beam 
splitter optics which frequency shifts one color light beam with respect to 
the other. This optical device has been incorporated in a current laser 
velocimeter system by Lockheed-Georg i a . 
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Figure 4.4 LV optical configuration. 
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A most Important aspect of the electronic processor is the data valida- 
tion circuitry. The original design, described in reference 58, was based on 
a digital comparison of the time for the first four fringe crossings with the 
time for the second four as well as an odd ratio comparison, three to five, 
of the fringe-crossing time count. Because of the digital nature of the com- 
parison, the acceptable time difference varied up to 3% of the period and at 
higher frequencies (velocities), where the fringe crossing period is very 
small, the accepted error was an increasing percentage of the period. This 
error window was found to be unacceptable and a new analog validation system 
was constructed. This error window is constant at ±1^ over a variable 

velocity range of thirty to one. This system was used in the current 

measurement program. 

In order to obtain satisfactory particles, the flow must be seeded arti- 
fically. It was found to be convenient to use different particles in the jet 
flow and the wind tunnel air. The jet air was seeded by introducing aluminum 
oxide particles, coated with CAB-O-SIL to reduce agglomeration, into the jet 
plenum. The particle size probability distribution peaked at 1 micron with 

an upper limit of 3 microns. The wind tunnel flow was seeded using an 

aerosol of 1010 hydraulic oil generated by a Laskin nozzle (ref. 60) . 


A. 2. 3 Experimental Program and Summary of Measurements 

The test conditions were designed to provide the maximum amount of use- 
ful data within the limits of the facility and instrumentation. The effects 
of two major changes were of interest. Firstly, how does the secondary 
stream velocity affect the jet development at fixed jet exit velocity, and 
secondly, how does the jet exit Mach number affect the jet development? In 
order to answer these two questions the experiments were divided into two 
sections. In the first the jet exit velocity was held constant and the wind 
tunnel velocity varied. In the second the jet exit Mach number was changed 
and the velocity ratio was kept constant. The test program is shown in 
Figure 4.5* In order to achieve a large range of velocity ratios and limit 
the spacing between the jet operating Mach numbers, the first series of 
measurements were carried out at Mj = .47 snd nominal velocity ratios of .1, 
.2, .3, .4 and . 5 . The second series of measurements were performed at a 
fixed velocity ratio of .1 and jet exit Mach numbers of .47> 1-37 and 

1 . 67 * The jet air was unheated in all the experiments. 

It was anticipated that the general effect of increasing both velocity 
ratio and Mach number would be to stretch the jet flow in the axial direc- 
tion. In order to obtain a representative set of turbulence measurements in 
series 1 (Mj = .47 for various velocity ratios), the following procedure was 
adopted. First a jet centerline traverse was carried out to give an indica- 
tion of the potential core length, Xq. Radial traverses were then performed 
at axial locations x^, x^/2 and x^/4. Radial traverses were also made at 
x = 8l.3 cm. and x=0.$l cm. These two locations were close to the limits 
of the measurement range set by the size of the LV optics and the dimensions 
of the working section window described in section 4.2.2. 

In order to reach the correct operating velocities in as short a time 
as possible, the following procedure was adopted. The correct values of 
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Figure h.S Jet turbulence test program. 
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static pressure ratio and tunnel dynamic pressure were calculated in the 
manner described in Appendix 2B . The jet and tunnel were simultaneously run 
up to their correct operating conditions on the assumption of a working see- 
t i on pressure equal to the ambient pressure. A door in the return section 
of the wind tunnel pressure, which was left open at the beginning of the 
run-up procedure to prevent pressure build-up in the tunnel as a result of 
excess mass flow injected by the jet, was then progressively closed until 
the static pressure in the working section was equal to ambient pressure. 
Minor adjustments to the tunnel speed and door position then enabled the 
correct test conditions to be quickly reached. 

At the time that a measurement was taken, a magnetic tape was written 
which recorded each individual velocity and its time of occurrence. Before 
analyzing the data, it was necessary to process this recorded data. The 
processing was required to correct for the natural biassing of the velocity 
data. Since the amount of fluid passing through the measurement volume was 
proportional to the local fluid velocity then, even if particles were uni- 
formly distributed in the fluid, a bias towards higher velocities occurred. 

It can be shown (ref. 6l) that if there is no dependence between the 
particle velocity and the particle number density then the true mean veloc- 
ity is found by taking the reciprocal of the arithmetic mean of the particle 
periods (time for the particle to cross eight fringes) multiplied by a 
constant of proportionality. It can also be shown that this biassing effect 
leads to a weighting of the velocity probability distribution function which 
is inversely proportional to the velocity. An example of the correction of 
a sample set of data is given in Appendix AD. A correction for this weight- 
ing of the probability function was made in the data processing, and the 
moments of the corrected distribution function were calculated giving the 
mean, standard deviation, skewness and kurtosis. At the same time any 
obvious noise points which had been validated by the LV processor were 
removed. A typical example of this procedure is also given in Appendix AD. 
The effects of particle size distributions are also discussed in Appendix AD. 

There are two characteristic parameters of the flow which were calcu- 
lated numerically from the mean velocity profile: the vorticity thickness 

and the dividing streamline radius. These two parameters as well as the jet 
half-width, the dividing streamline velocity, and the potential core radius 
were calculated at the same time. The calculation of the vorticity thickness 
involved finding the radial derivative of the mean velocity. In order to 
avoid unrealistic values of this derivative, the velocity profiles were 
smoothes prior to the calculation. 

Since no measurement of temperature in the flow was made, an approxi- 
mate relationship had to be used between the local mean density and velocity 
in order to calculate the dividing streamline radius. The dividing stream- 
line radius is defined as the radius at which 

^ds 

p Vj rdr = Pj "Y '■j^- 

/ 

o 
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assuming that the velocity and temperature are constant across the jet exit 
Using Crocco s relationship based on a Prandtl number of unity, it can be 
shown that the local mean density is given by 




(1 -Vj) 

ll -X) 


(y - 1) 

2 




(Vi -x)(i 



1 


Assuming constant total temperature, 
g i ven by 


the jet/tunnel temperature ratio is 


= 1 


+ — ' ^ ^ M 


/ 0-A2). 


With the density defined in this manner, the dividing streamline radius was 
read I 1y calculable. 

The velocity measurements will now be described. The analysis of these 
measurements using various scales will be carried out in section 4.2.4. 
However, before describing these measurements the influence of the tunnel 
walls, if any, on the development of the jet will be examined. 

The wind tunnel experiment was aimed at reproducing the flow conditions 
in and around a jet in an infinite stream of constant velocity. However 
since the tunnel working section had only finite cross-sectional area the 
tunnel walls were expected to influence the jet development. When a Jet is 
confined in a duct, there will be a region in the flow, before the jet 

boundary reaches the duct walls, where the entrainment of air from the duct 

flow into the jet is rapid enough to decelerate the freestream and establish 
a Posuive axial pressure gradient. The magnitude of this pressure gradient 
epends on both the ratio of the jet velocity to the duct velocity and the 

ratio of the jet diameter to the duct diameter. Hill (ref. 47) showed 

analytically that the behavior of jets confined in constant area ducts can 
^ function of only two independent variables, x/rj and 
m/ Mp)i, where rj is the radius of the duct, m is the total mass flow per 

iinii- average sum of momentum and pressure forces per 

axi;,ll»r by integrating the momentum equation across the duct and 

.k flow In the jet and duct in the exit plane of the 

nozzle, the value of m/(Mp)i in the jet exit plane is given, in the notation 
Or this report, by 


m 


(Mp) ■ 


^ + (1 -X) (rj/rj)^ 

(x2 + 2(1 - X2) (rj/r-r)2ji 


In the present experiments the equivalent duct radius to jet radius 
ratio was 20 to 1. The range of values of m/(Mp)i was between .84 for the 
lowest velocity ratio and .995 for the case of A =. 5 . Hill's (ref. ^7) cal- 
culations based on an integral approach show that in the first duct diameter 
tnere are some small changes in the duct velocity near the wall for values of 
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(m/Mp)^' of this order. The larger the value of tn/{Mp)i, the smaller will be 
the pressure gradient. 

Static pressure measurements were made in the tunnel wall for each 
velocity ratio and jet Mach number (measurements are not available for 

= A = .l due to a computer malfunction). These wall pressure measure- 
ments are shown in Figure ^.6 as differences from ambient pressure. Their 
relative location is a function of the accuracy with which the working 
section pressure was adjusted to ambient using the method described 
previously in this section. The maximum static pressure gradients are shown 
in Figure 4.7. It can be seen that there is little effect on the pressure 
differences when the jet plenum extension, mentioned in section 4.2.2, is 
removed. The trend of increasing axial pressure gradient with tunnel veloc- 
ity is contrary to Hill*s predictions. However, the pressure differences 
are small and Hill dealt with the problem of a circular duct with no wall 
boundary layers. Also, in the wind tunnel used in these experiments, slots, 
open to ambient conditions, are found in the tunnel walls at the end of the 
working section. This would also be expected to vary the pressure gradient 
from that predicted by Hill. The change in wall pressure was equivalent to 
a 1 % variation in free-stream velocity at the highest tunnel speed and a 
variation at the lowest tunnel speed. Thus, it can be concluded that for 
the range of experimental conditions covered, the wind tunnel free stream was 
a close approximation to an infinite uniform secondary flow. 

The variations of the jet centerline axial mean velocity and turbulent 
intensity for fixed jet Mach number and various velocity ratios are shown in 
Figures 4.8 and 4.9, respectively. Some basic trends in the data are 
immediately seen. As the velocity ratio X increases, so the potential core 
length of the jet increases, and the rate of decay of the centerline velocity 
with axial distance decreases. The development of the jet has been slowed 
with respect to time and since the jet exit velocity is constant the jet 
development occurs over a region of increased axial extent. The maximum 
turbulence intensity along the jet centerline also decreases with increasing 
velocity ratio and the location of the peak moves downstream. All the in- 
tensity traverses show a secondary, smaller peak, closer to the jet exit 
which may have been due to a regular vortex shedding from the jet lip. 

The variations of mean velocity and turbulence intensity on the jet 
centerline for fixed velocity ratio and various jet Mach numbers are shown 
in Figures 4.10 and 4.11, respectively. Again some basic trends are immedi- 
ately seen. Increasing the jet Mach number increases the potential core 
length and decreases the rate of decay of the jet centerline velocity. There 
also appears to be a clear division between the centerline behavior of the 
subsonic and supersonic jets. This difference is also seen in the axial 
behavior of the turbulence intensity. For the subsonic jets, there is very 
little change in the peak intensity level and only a slight downstream move- 
ment of the peak. The maximum turbulence intensity levels for the supersonic 
jet exit velocities is lower than the subsonic peak values and the location 
of the peak is further downstream. The secondary peak in the Mj=.47 
intensity data has disappeared for the higher jet exit velocities. This 
suggests another possibility, i.e., that the peak was caused by separation at 
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AXIAL DISTANCE x 


O Mj = -^7, X = .1 [7 Mj = .47, X = .5 

□ Mj = .47, X = .2 Q Mj = 1.37, X = .1 

^ ^ j ~ '^7, X = .3 Om. = 1.37, X = .l(Extens?on Off) 

OMj = .47, X = .4 OMj = 1.37, X = .1 

Figure 4.6 Working section wall pressure differences. 


117 



118 


i ng section pressure gradient as a function of tunnel Mach number. 
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Figure 4.8 Axial mean velocity on jet centerline, 
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Figure ^.9 Axial turbulence intensity on jet centerline, 






TANCE 



the main air supply valve, where the pressure drop for M| = .47 was from 
2x10^ N/m^ to 1.6x10^ N/m^ , causing regular pulsations in the air supply. 

It was indicated above that radial traverses were made at various 
multiples of the estimated potential core length. For the series I tests, 
(Mj = .^7i variable A) the axial locations at which the radial traverses were 
made are shown in Table I. 


Mj A (NOMINAL) AXIAL LOCATIONS OF RADIAL TRAVERSES (RADII) 


.^7 

.1 

3.25 

.47 

.2 

3.50 

.47 

.3 

3.75 

.47 

.4 

4.50 

.47 

.5 

5.25 


6.5 

15 

32 

7.0 

14 

32 

7.5 

15 

32 

9.0 

18 

32 

10.5 

21 

32 


TABLE I 


The mean velocity profiles and turbulence intensity profiles in the 
annular mixing region of the jet for Mj = .A7 and five velocity ratios are 
shown in Figure 4.12 through 4,21. The corresponding radial mean velocity 
and turbulence intensity profiles at x = 32 are shown in Figures 4.22 and 
4.23* Initial profiles of mean velocity and turbulence intensity at x=.2 
for Mj = .47 and A = .384 and .480 are shown in Figures 4.24 and 4.25, 
respectively. These measurements of mean velocity and turbulence intensity 
for a fixed jet Mach number and various velocity ratios will be examined in 
some detail in section 4.2.4. However, several points of explanation and 
comment will be made here: 

(1) The measurements for a velocity ratio of .1 are expected to be the 
least accurate for several reasons. Firstly, the tunnel speed was hard to 
keep constant at that low velocity. Secondly, the LV processor, in particu- 
lar the error circuitry, had to be readjusted for this case to enable 
velocity measurements to be made at the low speed. From the velocity proba- 
bility distribution functions, it could be seen that instantaneous veloci- 
ties well below the tunnel speed were encountered. It is possible that some 

truncation of the lowest velocities may have occurred. 

(2) The radial traverse, for Mj = .4? and A = l, was made at x=15 rather 
than x = 13 due to a positioning error. 

(3) The initial velocity profiles were made at x=.2 rather than x = 0 
for two reasons. Firstly, the two laser beams, which intersect at the 
measurement point to create a fringe pattern, approach this point at a small 
but finite angle. The closest location of the jet exit at which a traverse 
could be made was set by this angle. This is sketched below. 
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RADIAL DISTANCE r 


Figure A. 12 Axial mean velocity profiles, Mj = .^*7, A=.1 . 
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AXIAL TURBULENCE INTENSITY 



Figure 4.13 Axial turbulence intensity profiles, Mj ■ .47, A ■ . 1 
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Figure 4.16 Axial mean veloci ty prof i les, Mj = .47, X = .3 . 
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AXIAL TURBULENCE INTENSITY 
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Figure A. 19 Axial turbulence i ntens i ty prof i les , Mj = .47, A = .4 . 
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Figure 4.21 Axial turbulence intensity profiles, Mj = .47, 
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Figure 4.23 Axial turbulence intensity profiles, Mj = .47, x = 32. 
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Figure 4.2^ Initial velocity prof i les , Mj = .4? , X = .384, x = . 2 . 
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Figure A. 25 Initial velocity profile, Mj = .A7, A = .A80, x=.2. 
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Secondly, since particles only reach a location in the flow from inside the 
jet or in the wind tunnel, seeding is very uneven in the stagnation region 
at the lip. 

Since the lengths of the potential cores for the series II measurements 
(fixed velocity ratio, variable Mach number) were expected, in some cases, to 
stretch beyond the length of the working section, radial traverses in this 
measurement series were made at fixed distances from the jet exit. The 
locations were at x = ^, 8, 16, and 32. The mean velocity and turbulence in- 
tensity profiles for X =.1 and three jet exit Mach numbers are shown in 
Figures 4.26 through 4,31* The velocity measurements are tabulated in 
Appendix 4A. 

In the next section the data is reduced on the basis of several scaling 
parameters . 

4.2.4 Characteristic Dimensions and Scaling Parameters 

In this section the measured mean velocity and turbulence Intensity 
distributions, described in section 4.2. 3» will be used to calculate the 
variation of several characteristic properties of the flow. These parameters 
will be subsequently used as scaling parameters for the velocity distribu- 
tions. 

It was noted in the previous section that the effect of increasing the 
free-stream velocity was to stretch the mixing region in the axial direction. 
In Appendix 4B the decay of the centerline velocity is shown to follow the 
relationship 




1 - X 



(4-1) 


where = .Ok (1 - .32 X) x - .35. (4-2) 

The data reduced in this manner is shown in Figure 4.32. This leads to an 
expression for the potential core length, X(-, of 
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Figure 4.28 Axial mean velocity profiles, Mj = 1.37» X = .1 . 
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= 8.77/0 - -92 X) . 


(4-3) 


The empirical form for the variation of with A given by Forstal I and 
Shapiro (ref. 2) was 


= 8 + 24 A. 

This expression was obtained by fitting an expression of the form 


(4-4) 



(4-5) 


to their measurements. The same method has been applied to the present 
measurements and this is shown in Figure 4.33- It can be seen that even for 
the limited axial extent of the present measurements an expression of the 
form (4-5) fits the data adequately once a transition from the annular mix- 
ing region has been passed. Comparison of the expressions (4-1) and (4-5) 
shows that the former exhibits no discontinuity in the derivative of the 
centerline velocity at the end of the potential core whereas the latter does. 
Whenx=x^, (V^ ’ A) / ( 1 - A) = 1 , the exponential term inequation (4-1) can be 
written as exp{-A(A)/(x - x^) }, where A(a) is readily obtained from equations 
(4-2) and (4-3)* For a fixed velocity ratio the exponential may be expanded 
for large values of (x.-x^) giving, 


V 


a 


- A 


1 - A 


oc 



(4-6) 


which is significantly different from expression (4-5)- The variation of 
potential core length as a function of velocity ratio, found by these vari- 
ous methods is shown in Figure 4.34. There is a considerable scatter in 
Forstall and Shapiro's measurements. However, if equation (4-1) is compared 
with their measurements of the centerline velocity, there is fair agreement. 
This is shown in Figure 4.35* Since the empirical constants In equation 
(4-1) are dependent on the Mach number of the jet, some differences are to 
be expected. 

The variation of the potential core length with Mach number is shown in 
Figure 4.36. Witze (ref. 3) made use of a two-region turbulence model to 
predict the development of subsonic and supersonic free Jets. The empirical 
relationships he derived from a large number of experiments are shown in 
Figure 4.36. Clearly, there is some problem with the predictions for 
slightly supersonic Jet exit velocities. However, the general trends of the 
current measurements for A = .1 and these predictions are in good agreement. 


The rate of spreading of the jet is also defined at the axial variation 
of the Jet width. In the present investigation the width of the jet is de- 
fined as the vorticity thickness, 6^, which is given by 




(4-7) 
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Brown and Roshko (ref. /*1) point out that this definition of thickness is 
appropriate since the growth of the turbulent jet is governed by the un- 
stable motion induced by the vorticity. It has been found convenient for 
the purposes of presentation to consider the jet in two separate regions: 
the annular mixing region up to the end of the potential core and the flow 
region downstream of the potential core. 


It is to be expected that for sufficiently small ratios of vorticity 
thickness to potential core radius the initial mixing region will spread in 
the same manner as the two-dimensional mixing region. In fact, the spread 
rate of the axisymmetric jet agrees closely with that of the two-dimensional 
mixing region up to the end of the potential core. The variation of 6 /x 
for the two-dimensional mixing region is compared with the present measure- 
ments in Figure 4.37. It can be seen that the variation of mixing region 
thickness with axial distance agrees well for the two flows and can be 
simply related by 


6 


OJ 


const. 


(1 - X) 

(l + A) 


X 


(4-8) 


where Brown and Roshko (ref. 4l) give the constant in equation (4-8) a value 
of .181. Equation (4-8) is also shown in Figure 4.37. 

A second parameter which describes the rate of spread of the jet is the 
dividing streamline radius. This radius was defined in section 4.2.3. The 
location of the dividing streamline is shown in Figure 4.38 for the annular 
mixing region for the series I, Mj = .47 measurements. The value of the 
virtual origin of mixing, Xq , was found to vary between -.5 and 1.5. The 
calculated location of the dividing streamline as given by Korst and Chow 
(ref. 67) is also shown. It can be seen that these calculations for the 
two-dimensional mixing layer indicate a much slower divergence of the divid- 
ing streamline. The dividing streamline velocity is shown in Figure 4.39 
anjl is compared with Korst and Chow's calculations. In all cases, except 
A-.l, there is some agreement between the measured values of and the 
calculations for small values of x - x^. In the axisymmetric jet^case, the 
dividing streamline velocity decreases with axial distance. It should be 
remembered that, in the two-dimensional case, the mean velocity is similar 
whereas, in the axisymmetric jet case the mean velocity profile is not 
similar. 


In order to examine the effect of velocity ratio on the rate of growth 
of the jet downstream of the potential core, it is convenient to examine 
first the shape of ;he mean velocity profiles in this region. The local 
velocity difference normalized by the difference between the centerline 
velocity and the tunnel velocity is plotted in Figure 4.40 as a function of the 
radial distance divided by the half-velocity radius. It will be shown that 
there is a linear relationship between the half-velocity radius and the 
vorticity thickness so either can be used as the nond imens iona 1 i zi ng quan- 
tity. However, it is convenient to use the half velocity radius in the 
momentum equation below. It can be seen that there is good collapse of the 
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Figure 4.37 Variation of vorticity thickness with velocity ratio 





Figure 4.39 Variation of dividing streamline velocity with velocity ratio. 
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data except for the lowest velocity ratio, A = .096. Two analytical functions 
have been shown for comparison. The first is the cosine profile 


V^-A 




2 


1 + COS 



r < 2r 


* 5 


(A-9) 


This is the representation used by Squire and Trouncer (ref. 68). The second 
function is the error curve 

Vl-A 

= exp 

Va-A 

It Cdn be seen that both representations are a good fit to the data over most 
of the jet. The cosine curve, equation (^“9) > will be used in the subse- 
quent analysis. 

An expression for the variation of the jet half width as a function of 
axial distance can now be obtained from a simplified form of the momentum 
equation. If it is assumed that the tunnel wall radius is sufficiently large 
in comparison to the Jet radius so that there is no induced pressure gradient 
and the tunnel velocity is constant, then the momentum equation can be inte- 
grated to obtai n 

p rdr = const = y {pj + X^(r^^- l)} . (^"11) 

o 


^69315 


. 5 




(4-10) 


Assuming the mean velocity takes the form, 

Vi = X + (V3 - X) f (r/r^^) 

and letting 

then for the incompressible case, equation (4-11) can be written, 
r ^ 5 ^ {2(Va - A)A3i + (Vg - — 

where '^T 


Si = 


J 


f(n) ndn 


and 


6? - 


f2(n) ndn 


(^- 12 ) 

(4-13) 


(4-14) 


(4-15) 


i *;6 


o 



For a mean velocity profile of the form of equation (A-9) , these integrals 
are given by, Pi = .595 and 32® ‘345. The corresponding values for the 
profile (i|-10) are 3i».721 , 32 “. 361. Hill (ref. kj) used measured free 
jet velocity profiles to evaluate the integrals and obtained the values: 

3i -.7oo and 32 “-369. (The exact values depend on the ratio of the half 
velocity radius to the radius at the "edge" of the jet). Use will be made 
in this analysis of the first result corresponding to the cosine profile, 
equation (4-9). The axial dependence of (V, - A) is given by equation 
(4-1) so that equation (4-l4) becomes, 




2 


(1 + x) 


2X3ig(x) + (1 -X)329 ^(x) 


(A-16) 


where 


g(x) = 1 - exp {-1/25„}, 


(^-17) 


and is defined in equation (4-2). The axial variation of r 5, given by 

(4-6) tha^^^J^-^i/n earlier in equation 

® ® 9 (x) ^ A(x)/(x - Xc) for large values of ( i? - Xc) 

where A(A) was defined above equation (4-6) • For large values of x - x 
equation (4-6) shows that g(x) « (x-j^)-!. jhis leads to an interesting 
result for the asymptotic form of (4-1 6). For X equal to zero the half- 
veloci ty radius is linearly proportional to (x - x^) • For A not equal to zero 
and sufficiently large values of (x - x^) , r.5 is proportional to the square 
rwt of (x-Xj.). Close examination of the asymptotic form of (4-l6) shows 
that the function deviates from a linear growth with (X->tc) and approaches 
the square root dependence depending on the relative magnitudes of 
2X3i(x - Xj-)/(l -X)32A(X) and unity. The approach to these two limits is shown 
in Figure 4.4l and the present measurements are also shown to agree with the 
prediction. Szablewski (ref. 50) used a similar method to that used here. 
However, the asymptotic dependence is somewhat different, being r 5 <* x^/3 . 
This result was found by solving a first order differential equation in x 
for the centerline velocity which was obtained from an integral form of the 
mechanical energy equation with the mean velocity described by a shape 
function. In fact, Szablewski 's asymptotic dependence is more reasonable 
since it does corj-esj)ond to the axisymmetric wake flow value. However, for 
small values of (x - x^) , particularly close to the end of the potential core, 
equations (4-1 6) and (4-17) describe the Jet growth more accurately since 
the expression for g(x), the centerline velocity, is more accurate than 
that used by Szablewski. Clearly the accuracy with which the jet width is 
predicted depends on the accuracy of the centerline velocity prediction. 


Since the mean velocity profiles in the developed region of the jet 
have been shown to have a similar form, the vorticity thickness is to be 
expected to be linearly proportional to the half velocity radius. The 
vorticity thickness is plotted against the half velocity radius in Figure 
4-42. Also_shown are the relationships given by equations (4-9) and (4-10) 
which give 6^ = 4r.5/7r and 6;^ - 1 . 400r . 5 , respectively. The good agreement 
between the measurements and the ratio predicted by the error curve indi- 
cates that the slope of this function is closer to the measured slopes. 
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Figure k.k2 Relationship between vorticity thickness and half-velocity radius. 
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In the annular mixing region the mean velocity profiles can also be 
collapsed with the appropriate choice of coordinates. This does not apply 
to the mean profiles which occur in a region of flow establishment, where 
the wake velocity defect downstream of the jet lip sti 1 1 _exi sts . The mean 
velocity data in the annular mixing region in the form (V| ^\)/{yQ~X) as a 
function of (f -r^ 5 )/ 6 (^ are shown in Figure Though Vg is equal to one 

1 n most cases, some data just downstream of the lip of the potential core 
has been included. The collapse is seen to be very good. A good fit to 
the data is given by the expression, 

= y (1 - sin 2 n); - x 1 n $ X (A-l 8 ) 

Va-X ^ 


where 


^ “ ^5 



The simple expression (4-l8) has a derivative of -1 at n=0 which is a con- 
dition enforced on the data by the choice of the vorticity thickness as the 
local characteristic length. 

If the mean velocity profile is described in terms of a shape function 
and there is no pressure gradient, then two characteristic dimensions of 
the flow are required to specify the profile locally. These two parameters 
can be related through the momentum integral equation. Downstream of the 
potential core the centerline velocity and the half-velocity radius or 
vorticity thickness are the obvious choices. From equation (^-19) it can 
be seen that in the annular mixing region, the half-velocity radius and the 
vorticity thickness have been used. Either of these parameters could have 
been replaced by the potential core radius. The variation of the half- 
velocity radius and potential core radius with velocity ratio and axial 
distance is shown in Figure ^.4^. The potential core radius was defined as 
the radius at which (Vi"X)/(l - X) = .98. An analysis of the axial varia- 
tion of F ,5 and ?c» the core radius, can be easily carried out in the same 
manner as for the developed jet. It is interesting to note that the half- 
velocity radius appears to decrease with axial distance at some axial 
locations. This is not Inconsistent with the prediction of the half- 
velocity radius at the end of the potential core given by equation (4-l6) 
as, 


^ 1 + X 

2 2XBi + 0 -x)e 2 


X * Xf 


(^-19) 


The solution to this equation is also shown in Figure k.kk. The measure- 
ments at the end of the potential core also indicate that for the higher 
velocity ratios the half-velocity radius does lie within the jet lip line. 
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Figure 4.43 Radial distribution of axial velocity in the mixing region, 

Mj = .47. 
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The effect of change in jet Mach number for a fixed velocity ratio Is 
now examined. Equation k.S gives the effect of velocity ratio on the growth 
rate in the annular mixing region for uniform density. In Figure 4.45 the 
vorticity thickness is plotted as a function of axial distance for X 
nominally equal to .1 and various jet Mach numbers. The slope of the lines 
drawn in Figure 4.45 is called the spreading rate and, following Brown and 
Roshko (ref. 4l) , will be denoted by 6^*. For the two-^dimensional shear 
layer the spreading rate is related to Gortler's (ref. 69) parameter a by 
the relation, 

= ttVo. (4^20) 

Brown and Roshko (ref. 41) used an order of magnitude estimation to arrive 
at the dependence of the spreading angle on the Mach number. This depen- 
dence was of the form, 

(1 - X)^, (/*-2l) 

for high Mach numbers. In Figure 4^45 the variation of spreading rate 6^' 
has been plotted as a function of jet Mach number. Clearly an inverse 
dependence of the spreading rate on Mach number is evident for both 
measurements of A=0 and X = .l. The data collected by Birch and Eggers 
(ref. 70) indicates that for subsonic Mach numbers the spreading rate is 
constant. Thus, it is concluded that the expressions for the spreading rate 
of the annular mixing region of the jet are: 


^0)' ^ 0 - ^)/(l + A) Mj < 1 


(4-22) 


The variation of vorticity thickness at x = 32 for various Mach numbers 
is shown in Figure 4-4/. Both these results and the measurements in the 
annular mixing region suggest that for A the Mach number affects the 
spreading rate like Mj ^ for values of Mj greater than .7* However, the data 
points are too few to suggest any empirical relationship for this breakpoint 
and in the developed jet the vorticity thickness will be assumed to have the 
same Mach number dependence as in the annular mixing region, namely 

6^ ~ Mj > 1- (4-23) 

The dividing streamline radius, defined in section 4.2.3, is shown in 
Figure 4.48. The dividing streamline radius is seen to decrease with in- 
creasing axial distance in the supersonic cases in the annular mixing 
region. The half-velocity radius and the potential core radius are shown 
in Figure 4.49. The potential core lengths have been calculated using the 
constants given in Table III and equation (4B-22) of Appendix 4B . Clearly, 
the initial conditions are varying as a function of Mach number. No 
corrections for changes in virtual origin have been made. 
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Figure ti.kS Variation of spreading rate with Mach number. 
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Figure 4. A? Variation of vorticity thickness with Mach number 
at X = 32 , X = . 1 . 
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Figure k,kS Variation of half-velocity radius and potential core radius 
with axial distance, X=.l . 
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In Figure 4.50 the mean velocity profiles for X = .1 and various Mach 
numbers have been plotted. The profiles for Mj = .47 and Mj = .9 agree closely 
and follow a form similar to the profile obtained for various velocity ratios 
and Mj=. 47 in Figure 4.42. The corresponding function fit given by equation 
(4-l8) is also shown. However, the data corresponding to the supersonic jet 
exit conditions, though they agree closely with each other, do not follow the 
subsonic profiles. This is also the case downstream of the potential core, 
as is shown in Figure 4.51. The centerline Mach number for Mj = .47* *90, and 
1,37 is locally subsonic whereas that for Mj = 1.67 is locally supersonic. 
Anderson and Johns (ref. 75) have proposed a turbulence mixing model where 
the eddy viscosity in the supersonic region was different from that in the 
subsonic regime. Brown and Roshko (ref. 4l) also noted that at high Mach 
numbers a different model for the eddy viscosity could be used. The value of 
the eddy viscosity in a model used by Witze (ref. 3) was also found to take 
different values in subsonic and supersonic regions. This would lead to 
different profiles depending on the local conditions. 

The change in the turbulence level with velocity ratio and Mach number 
will now be examined. Two possible scaling velocities will be considered. 

The first is the usual velocity difference and the second is the difference 
between the Jet centerline velocity and the dividing streamline velocity. In 
Figure 4.52 the maximum axial turbulence intensity is plotted as a function 
of both these velocity differences. A greater amount of scatter is to be 
expected in the dividing streamline velocity difference due to the extra com- 
putations. A good fit to the data plotted as a function of velocity 
difference (V 3 - X) is given by 

V, = .053 + .118 (V^ -A). (4-24) 

^max 3 

The effect of Mach number on the maximum axial turbulence intensity is 

shown in Figure 4.52. It can be seen that the jet Mach number has little 

effect on the value of Vi . This confirms Brown and Roshko's supposition 

'*'max 

that the value of Vi is a function of the velocity difference only. 

^max 


For the case of an axisymmetric jet the variation of with velocity 

difference for small velocity differences can be inferred from measurements 

of axisymmetric jets into still air far downstream of the jet exit where the 

jet centerline velocity is very small. The measurements by Wygnanski and 

Fiedler (ref. 4) show that by 40 diameters downstream of the jet exit the 

value of vi is a constant fraction of the centerline velocity, namely 
^max 

= *3 (Va - A) ; (V 3 -A) < .163. (4-25) 


In Figure 4.54, Wygnanski and Fiedler's and the present measurements of 
Vj^ax been plotted as a function of (V 3 - X) . It can be seen that the 

curve can be divided into two regions. The first given by Equation (4-25) 
and the second given by, 


V 


Imax 


= .175 (Va " 


(Va - A) > . 163 . 


(4-26) 
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Figure 4.50 Mean velocity profiles in annular mixing region, X=.1 . 
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The slope of the linear portion of the curve, equation {k-2S) , and hence the 
crossover between the two regions will be a function of the shape of the wake 
generating body. Various asymptotic values of vi/(Va~A) have been noted by 
Antonia and Bilger (ref. A6) . 

No measurements have been made of the radial and azimuthal velocity 
components in the present study. However, it will be assumed that the 
radial velocity fluctuation is proportional to the rate of spread of the jet 
and the jet velocity. This assumption is bas ed on obtaining an estimate of 
the order of magnitude of the covariance VjV 2 from the momentum and continu- 
ity equations using, as we have shown above, the relationship between Vj and 
the velocity difference. It will be further assumed that the azimuthal 
velocity fluctuation is related to the other two components by [Townsend, 
ref. (76)] 

= (vi^ + V2^)/2. (4-27) 


Thus, the total mean square turbulent velocity is given by, 

q2 = 3 (vj2 + V2^)/2. 


(4-28) 


Using the expressions (4-22) for the spreading rate, it can be shown that 

.6 


q2 = 1 vj2 


1 + 


(1 -A) 
(1+A)2 


Mj < 1, 


(4-29) 


and 


, 2 . 2 . 


• 6 

1 .11^- 

M^d +x)2 




Mj > 1. 


(4-30) 


In the next section a numerical method for describing the structure of a jet 
in a moving stream will be discussed. 


4.3 NUMERICAL STUDY OF THE EFFECT OF A SECONDARY 
STREAM ON THE STRUCTURE OF A TURBULENT JET 


In this section a well-established numerical technique is applied 
to the problem of the jet flow in a moving stream. The various analyti- 
cal and numerical solutions that have been previously used will first be 
briefly reviewed. 

4.3*1 Summary of Previous Work 

There are four main classes of approach which have been used to analyze 
the jet boundary- layer problem. They are: 
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(i) similarity analyses, 

(ii) integral methods, 

(iii) eddy viscosity models, and 

(iv) turbulent kinetic energy models. 

In most cases these techniques, particularly the more simple analytical 
methods, have been used to treat special cases of free-shear flow such as the 
two-dimensional mixing layer or the free jet. However, some of the basic 
assumptions of each method can be considered to be relevant to the problem of 
a round Jet in a moving stream. 

The first method, the similarity analysis, makes the assumption of jet 
similarity. This does not imply that all components of the flow conform to 
the rigid definition of self-preservation but simply that the mean flow be 
similar. These methods have been widely used and are described in the works 
of Abramovich (ref. 77), Pai (ref. 78), and Schlichting (ref. 79)- The mean 
flow properties are described in terms of a single similarity coordinate 
which grows with the developing flow. This technique does imply that the 
local flow has no dependence on the initial conditions and as such, is only 
applicable in fully-developed regions of the flow where, in the case of jet 
flow, for example, the conditions at the jet exit no longer affect the jet 
development. A virtual origin of the flow may also be identified and 
measured experimentally which accounts for changes in the initial flow 
cond i t i ons . 

The second approach, the integral method, has been applied to the case 
of jets in a co-flowir.g stream. This general class of analysis uses a form 
of the boundary layer equation integrated across the jet. This method is 
used to calculate the downstream development of certain jet scaling parame- 
ters, such as the jet width, defined in a number of ways, or the jet center- 
line velocity. The equations of motion are usually written in terms of a 
similarity variable, and the shape of the mean velocity profile must be 
provided. As such, this method has limited usefulness since it is only 
applicable in regions where the mean flow is similar and it does not 
calculate the mean velocity profile but rather requires it as an input. 
Squire and Trounce r (ref. 68) integrated the momentum equation between the 
jet centerline and a radial location within the flow. The shear stress at 
this location was evaluated using the mixing length theory. Hill (ref. 47) 
used a imiment -of -momentum integral equation which was obtained, for the 
axisymmetric jet by multiplying the momentum equation throughout by the 
square of the similarity coordinate and integrating across the flow. The 
resulting integrals were evaluated using measured mean velocity profiles of 
a jet exhausting into still air. Szablewski (ref. 50) also used an integral 
approach to describe the motion of a heated jet into still air. He assumed 
that the mean velocity and temperature profiles could be represented by the 
function, (1 where the definition of n was different for the two 

properties. Integral methods have also been used by Patel and Newman (ref. 
80) , Gartshore (ref. 8l), Vogel (ref. 82), Bradbury (ref. 42), Bradbury 
and Riley (ref. 43), and Antonia and Bilger (ref. 46). Except in those 
works where the shear stress is measured, an assumption is made as to the 
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relationship between the shear stress and the local mean shear. This is 
achieved by the adoption of a hypothesis for the eddy viscosity or the turbu- 
lent Reynolds number which will be described below. Peters and Phares (ref 
b3) used a hybrid integral approach where the integrated turbulent kinetic 
energy equation was solved simultaneously with the integral equations of the 
mean flow. The turbulent shear stress was linearly related to the turbulent 
kinetic energy. It can be seen that though integral methods provide a simple 

development of the shear flow they cannot describe the 
whole flow region and do not permit calculation of radial flow profiles An 
integral approach has been used in section 4.2.4 to describe the relationship 
etween the jet width and the jet centerline velocity downstream of the 
potential core. 


The third class of analytical approach concerns the use of an eddy 
viscosity model. In the time-averaged boundary layer equations, the terms 
which present the greatest problems are the time-averaged products of turbu- 
lent fluctuations such as the Reynolds shear stresses and the product of 
velocity and temperature fluctuations which occur in the thermal energy equa- 
tion. Unless these products are evaluated experimentally, some hypothesis is 
required which relates them to the mean properties of the flow. Such an 
hypothesis is that original I y proposed by Boussinesq (ref. 84) which is now 
commonly referred to as the eddy viscosity hypothesis. The eddy viscosity 
was used to relate the turbulent shear stresses to mean shear in the same 
manner as the shearing stress in a laminar flow is related to the local rate 
of strain through the dynamic viscosity of the fluid. However, it is still 
necessary to relate the eddy viscosity coefficient to properties of the flow 
An attempt to do this was made by Prandtl (ref. 85) in his mixing length 
theory. This relates the kinematic eddy viscosity to the square of a length 
whose value must be established for the particular flow problem, multiplied 
by the modulus of the local mean shear. This meant that the shear stress 
would be zero at points in the flow where the local mean shear was zero A 

second hypothesis was proposed by Prandtl (ref. 86 ) which stated that the 
kinematic eddy viscosity was proportional to the width of the mixing region 
and the maximum difference in the mean velocities of the flow. 


Gartshore (ref. 81 ) used Townsend's (ref. 76 ) large eddy hypothesis to 
compute local values of the eddy or turbulent Reynolds number. This method 
related the change in turbulent Reynolds number from a self-preserving jet 
to a wake flow through the ratio of the strain rate of the mean flow at some 
typical station in the outer region of the jet. 


There have been many proposed eddy viscosity models of this type for 
both compressible and incompressible flows. These include those of 
Kleinstein (ref. 87 ), Warren (ref. 88 ), Alpinieri (ref. 89 ), Ferri et al 
(ref. 90 ), Witze (ref. 3), and Zelazny et al (ref. 91). The solutions 
obtained using this method agree very well, in each specific case, with the 
experimental results for the equivalent problem. However, no satisfactory 
universal constant has been found which binds together all the models. This 
method is simple to use and does give good agreement with experiments if the 
appropriate eddy viscosity model can be found for each flow region. 
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The fourth class of analytical approach is that involving the use of the 
turbulent energy equation. The use of this equation is once again for the 
purpose of describing the turbulent shear stresses. It is hypothesized that 
the turbulent kinetic energy controls the magnitude of the eddy viscosity. 

This suggestion was originally also due to Prandtl and has since been con- 
sidered by Glushko (ref. 92), Harlow and Nakayama (ref. 93_) Beckwith and 
Bushnell (ref. 9^) , Mel lor and Herring (ref. 95), and Spalding (ref. 96). In 
these models the eddy viscosity becomes proportional to the square root of 
the turbulent kinetic energy. Townsend (ref. 76) made an alternative sug 
gestion in which the shear stress was assumed to be proportional to the 
turbulent kinetic energy. Bradshaw et aZ (ref. 97) used this hypothesis 
successfully in the two-dimensional boundary layer. This hypothesis was also 
the one used in the integral approach of Peters and Phares (ref. 83 ) described 
above. The introduction of the turbulent kinetic energy equation does re- 
quire the specification of a number of relationships which describe the 
diffusion and dissipation of Reynolds stresses (this is described in detail 
by Nash and Patel (ref. 98 ). Harsha and Lee (ref. 99) and Lee and Harsha 
(ref. 100 ) have shown that the constant relating shear stress to turbulent 
kinetic energy applies over a very wide range of jet flows. However, the use 
of this equation does give the local shear stress a dependence on the up- 
stream history of the flow. The solution of the equations derived can be 
achieved using a finite difference approach and downstream marching since the 
equations are parabolic and form an initial value problem. The most compre- 
hensive study of this approach for compressible jets has been carried out by 
Heck and Ferguson (ref. 101) and Heck and Merkle (ref. 102) which also allows 
for the presence of shocks in the flow. 

The numerical method on which Heck and Merkle (ref. 102) based the tur- 
bulent mixing solution was that originated by Patankar and Spalding (ref. 
103 ). The GENMIX programs have been used extensively to solve many turbulent 
mixing problems. The method is based on the simultaneous solution of a para- 
bolic system of equations and the number of equations used is not limited. 
Various multiple equation models have been proposed and equations for all the 
Reynolds stresses ujuj (ref. 10^4), the decay rate of turbulence energy and 
the turbulence energy (ref. 105, IO 6 ). 

In the next section a two equation model of the turbulence which uses a 
transport equation for the dissipation rate is described. 

4 . 3.2 Description of the Numerical Method 

In a turbulence model which makes use of a conservation equation for 
the turbulent kinetic energy, it is necessary to provide several empirical 
relationships. Among these is the definition of the dissipation length 
which is required to enable the dissipation term in the T.K.E. equation to 
be evaluated. This length is usually linearly related to the local width of 
the shear flow and the constant of proportionality varies depending on 
whether the flow is plane or axl symmetr i c. The measurements of Antonia and 
Bilger (ref. 46) indicated that for the axisymmetric jet-wake, there were 
large changes in the ratio of the dissipation length to local thickness. 

This suggests that a turbulence model based simply on the kinetic energy 
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equation is not adequate for the problem of a turbulent jet into a moving 
stream. For this reason a system of equations is solved which includes an 
equation for both the turbulent kinetic energy and the turbulence energy 
decay rate. 

The analysis and the numerical technique used have been described in 
great detail by Pantakar and Spalding {ref. 103) and Launder et- aZ (ref. 

107) and no reiteration will be given here. However, it has been found that 
the method of defining the initial conditions Influences the solution. The 
method for obtaining these initial conditions from measurements will be dis“ 
cussed here. 

In many cases all the Initial values of the dependent variables will 
not have been measured. The parameters to be described are the axial mean 
velocity, the stagnation enthalpy, the turbulence kinetic energy and the 
turbulence dissipation rate. Measurements of axial mean velocity are 
usually available and the stagnation enthalpy can be calculated from these 
measurements and temperature measurements or by use of a Crocco's relation to 
obtain the temperature. If measurements of the covariance <VjV 2 ^ are not 
available, they can be estimated from the relation 

<V|V 2 > = Pt ^Vj/ar (^-31) 

where is an assumed eddy viscosity. The turbulence kinetic energy can 
then be estimated from the relation 

k = |<ViV2>|/.3. (4-32) 

The energy dissipation rate is used in the definition of the eddy viscosity 
and this relationship can be inverted to estimate the dissipation rate as 



where Cy and the remaining empirical constants required for the solution 
are given by Lauder et al (ref. 107). 

Using these definitions of the initial values for the dependent vari- 
ables, the flow field can be obtained by numerical integration of the con- 
servation equations. In the next section some solutions are described and 
it will be seen that some expertise is required in the choice of initial 
cond i t i ons . 


4 . 3.3 Numerical solutions and data comparisons 

In this section two numerical tests will be described. The first is a 
test case [case 9, (ref. IO 8 ) ] where the definition of the initial cond i - 
t ions has been given by the authors of the numerical technique. The second 
case is a use of the solution with the starting conditions taken directly 
from the current measurements. 
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In the first test case the mean velocity profile was measured. Helium 
was added as a tracer to the jet flow. The measurements were made by 
Forstall (ref. 109). The initial turbulent kinetic energy was calculated 
from the mean velocity data using equations (^“3l) and (4-32). The dissipa- 
tion rate was then found from equation (4-33)- The value of y ^ used in the 
calculation of the initial profile was that used by Launder et al (ref. 107) 
and was given by 

Pj = 3-2 X 10“^ prj AVi, (4-34) 

where AV^ was the velocity difference between the jet exit velocity and the 
minimum velocity and the freestream velocity and the minimum velocity in the 
inner and outer initial boundary layers, respectively. The initial grid 
values of V], k and c are shown in Figure Using these starting condi- 

tions the distribution of mean parameters in the jet was calculated. The 
centerline velocity decay of the mean velocity and the mass concentration of 
helium is shown in Figure ^.56. The agreement between the predictions and 
the measurements of Forstall (ref. 109) was good. 

The second test case used starting conditions taken from the current 
measurements for Mj = . 47 and X = .382. in order to examine the effect of 
initial grid spacing on the solutions, two sets of starting conditions were 
used. The turbulence kinetic energy and energy dissipation rate were calcu- 
lated from the initial mean velocity profile in the same manner described 
above for the first case. The initial mean velocity profile for a 20-point 
grid and a 27“point grid are shown in Figure 4.57* The extra 7 points have 
given greater definition to the region of rapid velocity variation. The 
predicted variation of jet centerline velocity is shown in Figure 4.58. It 
can be seen that the decay is over-predicted by the finer initial grid 
whereas the agreement between the 20-point grid calculation and the measure- 
ments is good. Similar trends are shown In Figures 4.59 and 4.60 for radial 
mean velocity profiles at x = l8 and 32, respectively. The agreement between 
the measurements and the more widely spaced finite difference grid is quite 
good in all cases. The variation of the turbulence kinetic energy along the 
jet centerline is shown in figure 4.61. The prediction has been compared 
with the measured values of 3 vi^/2, the turbulence kinetic energy for 
isotropic turbulence. Agreement between the absolute magnitudes cannot be 
expected because of the known anisotropy of the flow; however, the prediction 
of the axial location of the peak turbulence kinetic energy is seen to agree 
well. 


In this section it has been shown that the prediction of the time- 
averaged properties of a round jet in a moving stream may be numerically 
predicted using the method of Spalding and Patankar. However, it is clear 
that considerable expertise is required in the choice of initial boundary 
conditions if a good prediction is to be obtained. 
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Figure ^.SS Centerline decay of axial mean velocity and mass concentration of helium 
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Figure 4.57 Initial mean velocity profiles, A = .382, 
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Figure 4.59 Radial velocity profiles at R = l8. 
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Figure 4.60 Radial velocity profiles at x — 32. 
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Figure 4.6l 


Jet centerline variation of turbulent kinetic energy. 
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4.4 ACOUSTIC SCALING LAWS 


In this section the relationships between the flow parameters and the 
velocity ratio derived in section 4.2.4 will be used to estimate the varia- 
tion of the radiated noise due to alteration in the noise source magnitude 
and distribution. In Appendix 4C another model of this kind has been 
developed based on the properties of the dividing streamline. In this 
section attention will be paid to the more usually used scale of the velocity 
difference. One major objection to attempts to scale the measurements with 
the difference between the jet centerline velocity and the freestream veloc- 
ity has been the apparent unrealistic prediction of zero turbulence level 
and spreading for zero velocity difference. The existence of a mixing 
process for equal freestream and jet velocities in a realistic situation is 
due to the boundary layers on the inner and outer surfaces of the jet nozzle 
and the finite thickness of the jet lip. In the absence of these real 
phenomena, there would be no mixing for zero velocity difference. Thus, in 
theory, the turbulence level in the jet and the spreading rate could be 
related to the velocity difference between idealized freestream and jet 
velocities. In section 4.2.4 it is shown that for zero velocity difference 
a finite turbulence level and a finite rate of spread are predicted from 
data for lower velocity ratios. However, when dealing with a mean velocity 
profile where the minimum mean velocity is less than the freestream velocity 
such as that which exists in the vicinity of the jet exit, the real velocity 
difference, which governs the mean shear and hence the vorticity in the 
inner and outer flow regions, is not Vj " Vj but the difference between the 
maximum and minimum velocities. From the Initial velocity profiles shown in 
Figures 4.24 and 4.25, values of of .144 and .154 are found for maxi- 

mum velocity differences of .76 3nd .83, respectively. Using equation 
(4-30) values of of .144 and .154 exactly are fortuitously predicted. 

If the prediction had been made on the basis of the velocity difference be- 
tween the jet and the freestream values of of .125 and .111 would_^have 

been obtained. Clearly the use of the maximum velocity difference Vj-Vmin 
is more appropriate. 

The method to be used for noise estimation is based on the original 
dimensional analysis of Lighthill (ref. 110) and follows the same lines as 
that used by Cocking and Bryce (ref. 23) though the scaling laws for the 
flow variables will be those derived in section 4.2.4. I t wi 1 1 be assumed 
that the effect of freestream velocity on the radiated noise will be typi- 
fied by its effect on the initial development of the jet not too far down- 
stream of the potential cone. This limitation, which may still be repre 
sentative of conditions further downstream, should include the major noise 
source region. 


The noise per unit volume of the jet is given by 

T2 03^^ 

noise/unit volume —p-* , 

Po ' 


(4-35) 
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where a» Is a typical sound frequency, and V Is the volume of an assumed eddy 
The stress term T can be related to the totll turbulence level by 


T - Pn q 


o '^max 


(4-36) 


Assuming that the velocity difference is not too small so that vu^.. Is given 
by equation (4-26), equation (4-29) can be used to give 


^max 


, (1 - A) 
(1+X)2 


.6 


(1 -X) 


1.4 


(4-37) 


The volume of a typical eddy will be assumed to be related to the local 
thickness of the shear layer so that 


Ve ' 


(4-38) 


It will also be assumed that the characteristic frequency of an eddy noise 
source is given by. 


^^max (1 - X)*^ 


Vj- 


(4-39) 


Substitution of these relationships into equation (4-35) leads to 


noise/unit volume - 


.6 


1 + ii---") 
(1 +^)2 


(1 -x)^ ‘^ 
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ao' 


5‘ 


(4-40) 


The cross-sectional area of the mixing region will be approximated by 
so that the noise per unit length of the jet is given by 


.6 


noise/unit length ~ p_r i |l + i-Llii 

° M (1 +X)2 


(1 -X) 


5.6 


a 5 
®o 






If the total noise producing length of the jet is characterized by the 
potential core length, then the total noise radiated is given by 


1 + - 0 

0 +X)2 


total noise ~ 

Equation (4-3) can be used to obtain 


total noise 


Po^J' 


.6 


1+ illil 

(1 + X )2 


^ (1 -a)5-6 


V.8 

Ti - .92 X) ''J ■ 


(4-42) 


(4-43) 


Since the noise measurements at 90® to an isothermal jet are practically free 
fr^ any convective or refractive effects, the predicted variation of total 
noise with velocity ratio can be compared with the measurements at this loca- 
tion. The change of the noise level relative to the level when X»0 is 
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plotted as a function of velocity difference in Figure A. 62. The prediction 
given by equation (4-^3) tends to slightly over predict the effect of velocity 
difference on the radiated noise at small velocity differences. The exponent 
of the velocity difference from the noise measurements was 5*5* The exponent 
given by equation (4-43) is readily calculated from equation (4-43) and is 
shown in Figure 4-63. The exponent rapidly approaches the measured value of 
5-5 from its value at X=0 of 7*28. The difference between the prediction and 
the measurements never exceeds 1.5 dB and, since the exponent is less than 
5 . 5 , the measured value for values of A greater than one half, 
the difference becomes smaller at higher velocity ratios. However, 
since the relative velocity effect is only of the order of a few dB, the 
prediction could be improved. 


The most probable cause of the inaccurate prediction is the form chosen 
for the variation of the turbulence level with velocity ratio. More speci- 
fically, the estimated variation of the radial and azimuthal turbulence 
intensities with velocity ratio appears to be too great at small values of A 
These variations are the only ones for which adequate experimental informa- 
tion was not available. It is of interest to note that if variation of 
spreading rate, and hence radial velocity fluctuation, of the form suggested 
by Yule (ref. 38 ) is used, namely. 


6 


I 


U) 


(1 -A) 

(1 + A)i 


(4-44) 


a better fit to the data is obtained. Use of equation (4-44) results in a 
noise prediction of the form. 


total noise 



1 + 


(1 - X)- ^ ^ (1 - 

(1 + a) (1 - .92 X) 


V/. 




The result for this prediction is also shown in Figure 4.62. The variation 
of the corresponding velocity difference exponent with velocity difference is 
shown In Figure 4.63. The small change in the expression for the variation 
of the spreading rate clearly has a significant effect on the resulting noise 
prediction. The variation of spreading rate and radial and azimuthal turbu- 
lence levels with velocity ratio for small values of velocity ratio is thus 
seen to be very important. 


4.5 SUMMARY AND CONCLUSIONS 


An experimental investigation of the effect of a secondary stream on the 
turbulent structure of a round jet, and the noise radiated by the jet, has 
been performed. The major highlights and conclusions are given below. 

( 1 ) A two- inch diameter jet was mounted on an aerodynamical ly faired 
plenum and support in a 0.76 m x 1 .09 m low-speed wind tunnel . Tunnel 
velocities of up to 76.2 m/s could be achieved. Static pressure 
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Figure 4.62 Predicted variation of radiated noise with 
velocity difference. 
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tappings were located along one wall of the working section. A window in the 
opposite working section wall allowed observation of the jet flow up to 0.86 m 
from the Jet exit. The static pressure in the tunnel was controlled by a door 
in the tunnel return section. 

( 2 ) Wall static pressure measurements indicated a maximum variation of 
two per cent in the working section tunnel velocity. 

( 3 ) A laser velocimeter system was used to perform the velocity 
measurements. The velocimeter optics were mounted on two lathe guides and 
a hydraulic table which allowed traversing in three directions. The system 
operated in single-channel back-scatter mode enabling measurements of instan- 
taneous axial velocity to be made. The Jet supply air was seeded with one 
micron aluminum oxide particles and the wind tunnel air was seeded with a 
hydraulic oil aerosol. 

(4) Two series of test points were considered. The first were for 
fixed Jet exit Mach number of .4? and velocity ratios of . 1 , . 2 , . 3 , .4, and 
. 5 . The second series were for a fixed velocity ratio of .1 and Jet exit 
Mach numbers of .47, .9, 1.37, and 1.67. All experiments used unheated air. 

( 5 ) for the fixed Jet exit Mach number measurement series, radial 
velocity traverses were performed at axial locations of x^/k, X(./2, x^ and 

“ 32 , the potential core length being determined by a centerline axial 
traverse. In the second measurement series for variable Mach number, radial 
traverses were performed atx/rj=4, 8, 16, and 32. 

(6) The Jet centerline velocity was found to decay as. 


(Va - X) 


exp 


.08 (1 - .92X)x - .7 ’ 


downstream of the end of the potential core. 

(7) The potential core length for M=.47 is given as a function of X 

as, 

x^ = 8. 77/(1 - .92X). 


(8) For fixed velocity ratio, .1, and variable Jet Mach number, the 
potential core length followed the same variation as a Jet exhausting into 
stationary air, being almost constant for Mj < 1 . 

(9) In the annular mixing region of the Jet, the rate of spread of the 
mixing layer agreed closely with that of a two-dimensional shear layer. 

(10) For supersonic Jet exit velocities, the rate of spread was in- 
versely proportional to the Mach number. 
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(11) The axial mean velocity radial profiles in the annular mixing 
region could be described by a shape function whose dimensions were fixed by 
the half-velocity radius and the vorticity thickness. This only applied 
after the region of flow establishment where the effects of the initial 
boundary layers and the finite jet lip thickness were no longer present. 

(12) The radial distribution of axial mean velocity downstream of the 
potential core could also be described by a shape function and the local 
vorticity thickness and centerline velocity. 

( 13 ) The growth of the jet downstream of the potential core region was 

related to the centerline velocity decay through the momentum integral equa- 
tion. The resulting relationship showed that the width of the jet was 
linearly proportional to axial distance for X=0 and large values of x, but 
was proportional to x 2 for A not equal to zero. It was noted that in 
reality the asymptotic limit would give an relationship. 

(14) it was observed that the mixing charac ter i s t i cs of locally super- 
sonic flow were different from that for locally subsonic flow. This was 
evident from axial and radial mean velocity profiles. 

( 15 ) The peak axial turbulence intensity at any axial location was 
found to be given by, 

- (Va - A) 

^max ' ^ ^ * 

except at small values of (Vg - A) where the peak intensity was linearly 
proportional to the relative velocity. 

( 16 ) The axial turbulence intensity was found to be independent of the 
jet exit Mach number. 

( 17 ) Using a simplified model of the noise generation mechanism, the 
radiated noise at 90 ^ to the jet axis was found to be given by. 

Radiated noise at 90° ~ 1 + li” j (i - ' 

This agreed closely with the acoustic measurements obtained in the anechoic 
f ree-jet faci 1 i ty . 

(18) The relative velocity exponent at 90^ predicted by the expression 
above rapidly approached the measured value of 5-5 from a value of 6.3 at 

X =0. 

(19) The turbulence measurements indicated that the measured relative 
velocity reduction in noise at 90° to the jet axis, which is observed in 

f ree-jet and wind tunnel measurements, may be attributed to changes in the 
source level due to alteration of the turbulence structure by the secondary 
stream. 
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5. SYNOPSIS OF TECHNICAL WORK 




5,1 INFLIGHT SIMULATION EXPERIMENTS ON JET NOISE 


5 . 1.1 Anechoic Free-Jet Facility 


A new aneohoio free- jet facility y funded by Lockheedy has been designed and 
constructed for investigating forward motion effects on jet noise. The 
facility was subjected to rigorous aerodynamic and acoustic performance 
evaluation tests prior to the inflight simulation experiments on jet noise^ 


The facility was powered by a jet ejector and was capable of providing 
tunnel velocities up to 75 ni/s with a test section of .76 m by 1.07 Prior 
to the construction of the full-scale facility, an exact one-fifth scale model 
was built to confirm various aerodynamic performance concepts and to aid in 
the design of the free-jet working section as well as the shape of the collect 
tor. A planview schematic of the complete facility is shown in Figure 5*1- 
In order to provide low background noise levels, the tunnel ducting between 
the anechoic room and the jet ejector incorporated several sound absorbent 
sections. For minimum blockage (and therefore minimum flow disturbance) In 
the working section, the air supply ducting for the primary jet was installed 
axially in the intake/contraction section. The ducting was carefully designed 
to avoid any flow separation within the accelerating free-jet flow in the con- 
traction section. The upstream internal noise levels were minimized by the 
two mufflers connected in series. The microphone arc was at 5^ nozzle diame- 
ters, and was placed outside the free-jet flow. Noise data were normally taken 
in the range 3O®<0m<9O^, but forward arc measurements could be accomplished 
by moving the primary jet downstream. 

In order to confirm the design criteria and to ensure the accuracy of the 
subsequent jet noise measurements, the facility was subjected to rigorous 
performance evaluation tests, and the major findings were as follows: 

Flow visualization tests using a smoke generator established that the 
free jet was stable throughout its length, and the air-flow circulation veloc- 
ities in the anechoic room were negligible. The mean flow properties of the 
free-jet working section were examined quantitatively by using a specially 
designed pressure probe rake. The mean velocity profiles were mapped in 
detail, and the tunnel calibration was derived from these tests. 

The acoustic performance tests were designed to examine the anechoic 
quality of the facility, the background noise levels, and the rig Internal 
noise levels. It was established that (l) the facility was anechoic down to 
200 Hz, ( 2 ) the facility background noise in the measurement arena was low 
(see Figure 5.2), and (3) the jet noise results would not be contaminated by 
internal noise at least down to Vj/aQ = 0.32. (All symbols are defined in 
Appendix 5.) 

Limited calibration experiments using a point source established that in 
the frequency range from 1 KHz to 5 KHz, no detectable internal reflections in 
the free-jet test section were observed for tunnel velocities up to 60 m/s. 

It was also confirmed that at these frequencies, the effect of turbulence 
scattering in the facility was not significant. 
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5.1.2 Jet Noise Experiments 


The effects of forward motion on the characteristics (both directivity and 
sipectraZ) of turbuZent mixing noise from jet exhausts were examined by con^ 
ducting inflight simulation experiments in the anechoic free- jet facility. 
The resulting scaling laws, after applying facility corrections y were found 
to be in close agreement with the scaling laws derived from theoretical and 
semi- empirical considerations • 


In the jet noise experimental program, the free-Jet to nozzle area ratio 
was 400. The nozzles (d = 5«08 cm.) employed were a convergent nozzle, and 
M = l.^ and 1.7 convergent“di vergent nozzles. Acoustic measurements were con- 
ducted at forty different test conditions (combinations of Vj/a^ and Vj/ao) ; 
eight values of jet exit velocity were chosen in the range 0.4 < Vj/a^ < 1 . 345 , 
ana the tunnel velocity was varied up to Vj/aQ = 0.2. All acoustic data were 
carefully scrutinized prior to detailed analysis, and all data likely to be 
contaminated by extraneous noise sources (background/! nstrumentat ion/shock 
noise) were not utilized. 

The uncorrected results (i.e. data to which no facility corrections were 
applied) were found to be in good agreement with the results from published 
free-jet experiments. 

The measured results were subjected to a systematic data correction pro- 
cedure, which was derived in the present program, and which converted results 
from a free-jet facility to the corresponding results that would be obtained 
in a large-scale wind tunnel simulation. In essence, the correction proce- 
dure took proper account of source distribution effects in a jet flow, the 
downstream convection of sound waves by the tunnel flow, and the refraction 
of sound caused by the free-jet shear layer. Using the measured results at 
fixed measurement angles, 9^, it finally yielded results corrected to constant 
emission angles, 0^, for an observer moving with the nozzle. In general, 
the magnitudes of the facility corrections were small at 90° to the jet axis, 
but produced a significant effect at lower angles in the rearward arc. 

At 0g = 90 ° , the effect of tunnel velocity on the corrected spectra was 
virtually independent of frequency, and the spectra at various tunnel veloc- 
ities were nearly parallel. This implied constant reduction in equivalent 
source strength at all frequencies. A typical example for Vj/ao = 0.9 is 
presented In Figure 5.3* 

The corrected overall SPL results indicate that the magnitudes of the 
inflight noise reductions increase as the observer moves from 0^ = 90° toward 
the downstream jet axis (Figure 5*4). The reductions at all emission angles 
were scaled on the relative velocity basis. At 0^ = 90°, the OASPL is found 
to be proportional to 5.5 powers of the relative velocity, a result which 
agrees closely with the theoretical and semi-empirical scaling laws, derived 
purely from source alteration considerations. 
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Figure 5-^ Effect of tunnel velocity on the directivity of 
corrected overall SPL at Vj/aQ = 0.9. 


200 



The results from the present free- jet experiments are in very good agreement 
with published wind tunnel results. However ^ there are large differences 
between the inflight effects as observed from model simulation experiments 
and full-scale flight tests. These discrepancies are tentatively attributed 
to the doubtful validity and accuracy of the flight data. 


The relative velocity exponents, obtained from the present corrected 
results, are compared with the exponent values obtained in the NGTE wind 
tunnel experiments (where the test section to nozzle area ratio was 5000) in 
Figure 5*5. The agreement is good at all emission angles. From this compari- 
son it is concluded that (l) the effect of scattering caused by turbulence 
in the free-jet mixing layer does not appear to be significant; (2) the free 
jet to nozzle area ratio of kOO used in the present experiments was adequate 
for accurate evaluation of inflight effects on Jet noise; and (3) a free-Jet 
facility is capable of simulating many of the inflight effects on Jet noise 
providing adequate and accurate facility corrections are applied to the 
measured resu 1 ts . 

The comparison of results from inflight simulation experiments with 
full-scale flight results (Figure 5.5) shows that although the agreement is 
reasonable at low angles to the Jet exhaust, there are significant discrep- 
ancies at larger angles. At there is little or no change in the 

flyover noise levels, whereas the simulation experiments indicate significant 
noise reductions. All attempts to explain this discrepancy in terms of source 
motion (relative to a fixed observer) effects lead to the conclusion that the 
flyover results at 90° will not be affected. It is therefore concluded that 
the flight results considered here may not represent pure turbulent mixing 
noise. It is recommended that the "acoustic cleanliness" and the measurement 
accuracy of all flight results be examined thoroughly, both for existing 
flight data as well as any future flight test results. 

Finally, a theoretical formula for the prediction of inflight noise 
reductions at emission angles outside the so-called zone of silence has been 
derived, and the correlation between calculated and measured OASPL reductions 
is obtained to an accuracy of ±1/2 dB, as shown in Figure 5»6. 
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5.2 ACOUSTIC PROPAGATION THEORY 


5 . 2.1 An Acoustic Model of the Free~Jet Flight Simulation Facility 


A geometric acoustics model of sound propagation through the free- jet shear 
layer has been evaluated. The model yields analytic formulae for angle and 
amplitude corrections that convert the measured free-jet data to estimated 
flight data. The recommended correction procedure takes into account the 
axial distribution of the jet mixing noise source as a function of Strouhal 
number, 


In the Lockheed flight simulation facility, sound waves radiated from 
the primary jet are convected and refracted by the mean flow in the potential 
core and in the shear layer of the free^jet. The change in angle and ampli“ 
tude of the direct radiation across the shear layer must be measured or cal- 
culated in order to convert the free-jet measured data to estimated flight 
data. A theoretical approach has been chosen, based upon Geometric Acoustics 
(GA) or "ray-tracing." The GA model has been partially justified with an 
order of magnitude study of the terms in the homogeneous Li 1 ley equation 
which are neglected in the geometric approximation. 

Results from the GA model Indicate that axial variations and axial 
gradients of the mean flow velocity in the free-jet shear layer have a negli- 
gible effect on sound refraction and that the change in angle can be calcu- 
lated with the analytic, stratified flow equation to a very good approximation 
(better than 1° accuracy). Examples of ray paths on the x-r plane (cylindri- 
cal coordinates, origin at free-jet nozzle) are shown in Figure 5.7t before 
and after propagation through the free-jet shear layer; 0j is the wavenormal 
angle within the uniform flow and corresponds to the flight emission angle 
0^. The angle 6^ is the measured polar angle based on a microphone arc of 
radius 5^ d centered on the primary jet nozzle at x = 6d (d is the primary 
nozzle diameter). The ray origin or source location (X^Xg) is determined 
from an empirical model in which the effective center of the turbulent mixing 
noise source in the primary jet flow can be calculated as a function of the 
modified Strouhal number, S^. In Figure 5*7 the source location corresponds 
^m ” P^ths inside the shear layer, although calculated, have been 

replaced by an extrapolation to the "lip-line" of the entry and exit ray 
paths. In general extrapolated ray paths, as in Figure 5-7, intersect the 
"lip-line" at almost exactly the same point and therefore, to a good approxi- 
mation, ray paths can be drawn as if the shear layer were replaced by a 
vortex sheet at the "lip-line" as indicated in Figure 5*8. 

However, this vortex sheet model is not used to calculate amplitude 
changes across the shear layer, only ray paths and changes in angle. The 
sound amplitude, according to GA, varies smoothly through the real shear 
layer and the total change in amplitude is calculated from an analytic ex- 
pression based upon (i) the law of energy conservation in a ray tube and (ii) 
the change in ray-tube cross section across the shear layer which is defined 
by the ray paths as shown in Figure 5-8. 
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Free-jet shear layer can be replaced by a vortex sheet 
at the "lip-line," R = 1, for ray tracing purposes. Ray 
tube cross-section changes from dAj^ to UAq across 
shear layer; wavenormal angle from ej to 6 q. 


5 . 2,2 The Influence of Forward Motion on Flow-Acoustic Interactions 
as Described by the Li 1 ley Equation 


The influence of forward motion on flow-acoustic interactions has been ex- 
amined through a re-interpretation of the ^^static^* numerical solutions to the 
Lilley equation* Outside the cone of silence y the Lilley equation radiation 
levels change with forward motion in a way that can be calculated, to a good 
approximation y with the geometric acoustics analytic result* 


The influence of forward motion (simulated or actual) on the theoretical 
directivity of isothermal jet mixing noise can be summarized as follows. The 
expression for the proportional bandwidth mean square pressure in the static 
case is a product of two factors 

(Flow-acoustic interactions: V3) x (Eddy convective 

ampl i f i cation: V^) 

and, in the flight case, three factors 

(Flow-acoustic interactions: Vg - x (Eddy convective 
amplification: - Vy^) x (Dynamic effect: Vy^) , 

where Vy^ is the (uniform) aircraft flight speed. The eddy convection velocity, 
Vc, controls the static convective amplification (e.g. |l * V^cosOe/ao | “^) ; in 
the flight case it is the eddy convection velocity relative to the observer, 
Vc”Va> that determines the amplification (e.g. 1 1 “ (V^ - V^) cos0g/aol “^) . 

Both this effect and the dynamic effect are well established results. 

Similarly, the static flow-acoustic interaction effects are described by 
Lilley equation solutions based upon a mean velocity profile with a center- 
line velocity Vg and that is reduced to Vg = Vg - for the flight case 
although it is still a "static*' profile. Hence, forward motion effects can 
be examined through a re- i nterpretat ion of the "static" numerical solutions 
to the Lilley equation. Some typical numerical results outside the cone of 
silence for the variation of the flow-acoustic interactions factor, F^g], 
with Va^rel ^re shown in Figure 5 * 9 ; the Strouhal number, takes two 
values that are appropriate for the four diameter axial station profile used 
in this example. This emission angle 0g varies with Vg pg| in these results 
since a source emission angle, 0g, is held constant; 0g lies between 50^ and 
60*^ over the range 0.75 < Vg^pg]/aQ < 1 .25. 

The geometric acoustics analytic result is also shown for comparison in 
Figure 5 . 9 * While there are significant absolute differences between it and 
the Lilley equation numerical results, the slopes are very similar. There- 
fore, outside the cone of silence, the Lilley equation flow factor variation 
with forward motion can be calculated, to a good approximation, with the 
geometric acoustics analytic result. 
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5.3 EFFECT OF FORWARD VtLUClIY ON THE 
STRUCTURE OF A TURBULENT JET 

5 . 3.1 Turbulence Measurements 


Ve loaity measurements are made in a jet in a wind tunnel, with a laser velo- 
o^meter, to describe the effect of velocity ratio and jet eant Mach number on 
the development of a jet in a moving stream. 


The aim of this work has been to describe the structure of a turbulent 
jet, in the presence of a secondary stream, in relation to the noise- 
producing properties of the jet. 

A 5.08 cm. diameter jet was mounted on an aerodynami ca 1 1 y faired plenum 
and support in a .76 m x 1.09 m low- speed wind tunnel. The experiments were 
conducted for a range of wind tunnel and jet exit velocities shown in Figure 
5^0. One series of measurements, at fixed jet exit Mach number, examined 
the effect of velocity ratio on the jet development and the second series, 
at fixed velocity ratio examined the effect of jet exit Mach number on the 
jet structure. 

Measurements of wall static pressure in the tunnel working section 
showed the axial pressure gradient to be less than an equivalent freestream 
velocity change of 2% at the highest tunnel velocity. Measurements of the 
Instantaneous velocities in the jet were made with a laser velocimeter system 
mounted outside the tunnel working section. 

The variation of jet width with axial distance in the annular mixing 
region of the jet was found to agree with the variation in a two-dimensional 
mixing layer. The variation of potential core length with velocity ratio was 
found to be given by, Xj-=8. 77/(1 -.92 X) . This variation is shown in Figure 
5 . 11 , where it is compared to the linear variation of potential core length 
with velocity ratio proposed by Forstall and Shapiro (ref. 2 ). The method 
used by Forstall and Shapiro to determine the potential core length variation 
was applied to the present measurements and was found to overestimate the 
potential core length. The variation of potential core length with jet exit 
Mach number for a velocity ratio of 0.1 was found to be In agreement with the 
two region model prediction of Witze (ref. 3 ) for free jets. The radial 
mean velocity profiles in the annular mixing region of the jet and downstream 
of the end of the potential core were found to be represented well by shape 
functions using the half velocity radius and vorticity thickness and the jet 
centerline velocity as the characteristic properties of the flow. The 
collapse of the radial mean velocity profiles in the annular mixing region is 
shown in Figure 5 . 12 . 

For the most velocity ratios of practical interest the axial turbulence 
intensity was found to follow a relationship of the form, axial turbulence ~ 
(relative velocity)^-^ This is shown in Figure 5.13. For small velocity 
differences the turbulence level is expected to be proportional to the 
relative velocity. The axial turbulence intensity was found to be independent 
of the jet exit Mach number. 
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Figure 5*13 Variation of maximum 
turbulence in tens i ty 
with velocity 
d i f ference. 




5.3*2 Jet Noise Scaling 


The variation in radiated noise ^ at 90^ to the jet axisy with velocity ratio 
is predicted using turbulence measurements and a simplified Lighthill noise 
radiation model. 


The measurements of jet turbulence made in the subsonic wind tunnel 
experiments may be used to estimate changes in the radiated noise with 
velocity ratio. 

The Reynolds stress sources in a subsonic jet are assumed to be propor- 
tional to the ambient density and the total turbulence intensity. In order 
to calculate the total turbulence level, three components of velocity fluc- 
tuations are required. The axial turbulence intensity is obtained directly 
from the measurements. The radial velocity fluctuations are assumed to be 
proportional to the rate of spread of the jet and the azimuthal velocity 
fluctuation is assumed to be the average of the other two components. The 
source frequency is related to the axial velocity fluctuation and the local 
shear layer thickness. The eddy volume is taken as proportional to the cube 
of the local thickness. The variation in the noise producing volume is 
considered to be described by variations in the mixing layer thickness and 
the potential core length. 

Two prediction formulae are obtained depending on the description of the 
variation of local thickness with velocity ratio. Best agreement is obtained 
if the spreading rate of the mixing layer is assumed to be proportional to 
(l -A)/1+X)i. In this case the radiated noise at 90^ to the jet axis, where 
convection and refraction effects are at a minimum and changes reflect source 
alterations, is given by 


Radiated noise at 90° - 


This prediction is compared with the free-jet noise measurements in Figure 
5*1A. The predicted variation of relative velocity exponent is shown in 
Figure 5.15. It can be seen that, except at small values of velocity ratio, 
the agreement between the predicted exponent and the measured value of 5-5 
i s good . 
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Figure 5.1^ Predicted variation of radiated 
noise with velocity difference. 



Figure 5-15 Variation of velocity difference 

exponent with velocity difference. 
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6. RECOMMENDATIONS FOR FUTURE WORK 


The work conducted during the present program has examined many impor- 
tant aspects of forward velocity effects on turbulent jet mixing noise from 
a fundamental viewpoint. Several important recommendations for future work, 
necessary to achieve the ultimate goal of this research program (described 
in Section 1), have emerged during the course of this study. These recom- 
mendations are now discussed in detail. 

The recommendations, essentially fall into two categories. The first 
category is centered on the problem that significant differences exist be- 
tween the inflight effects as observed from the full-scale flight tests and 
those obtained from model-scale inflight simulation experiments conducted in 
static facilities like a wind tunnel or a free-jet facility. These discrep- 
ancies need to be resolved in a systematic manner, since the basic concept of 
predicting inflight noise levels from flight simulation tests becomes mean- 
ingful only when there is a good correlation between the flyover results and 
the corresponding forward velocity simulation results. The second category 
or set of recommendations form the continuation of the fundamental study 
that has been conducted in the present phase. These fundamental studies on 
the effects of forward velocity on jet noise are considered to be vital in 
achieving the long-term goal of being able to predict the sound field of a 
jet exhaust in flight without actually flying the jet configuration. 

It is strongly recommended that in order to reach the ultimate 
objective of this research effort in an efficient manner, the two categories 
of recommendations described above should be implemented with 6C[u/xZ emphasis. 
That is, both sets of recommendations should be followed concurrently if 
possible, and it could be detrimental if one set of recommendations is ex- 
ecuted at the expense of the other. 


6.1 FIRST CATEGORY RECOMMENDATIONS 

The recommendations necessary to reconcile the differences between 
flight results and model simulation results are discussed first. From the 
present study, it can be concluded that the observed discrepancies may have 
arisen due to one of the following two reasons: (l) a relative source 
motion effect, which is present in the flyover tests, but which is absent in 
the simulation tests where the observer is effectively moving with the jet 
nozzle; (2) the doubtful validity of the cleanliness of flight data, which 
may give rise to inflight effects that may differ considerably from those 
that are associated with pure turbulent mixing noise. The basic philosophy 
behind future work on this aspect lies in the assessment of the relative 
roles played by these two possible reasons in producing the currently ob- 
served differences. This can be achieved by conducting a carefully 
controlled experiment where the two phenomena can be isolated. To do this, 
it is recommended that an experiment involving the measurement of the noise 
field, by stationary microphones, from a moving jet should be conducted. 
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The jet conf igurat ion to be used in chis experiment should be free from in- 
ternal noise, and the experiment should be conducted at model scale in a 
carefully controlled environment where the measurements and corrections can 
be obtained accurately. The same Jet configuration should then be tested in 
a flight simulation facility, for example, in an anechoic free-jet facility, 
and the results from the two sets of experiments should be compared in 
deta i 1 . 

The outcome of this comparison will establish the exact nature of subse- 
quent research as follows: (l) If the results from the moving jet experiment 

agree with the correspond i ng results from the flight simulation experiment, 
then it can be concluded once and for all that the currently observed dis- 
crepancies arise due to some problems associated with a full-scale engine 
(e.g. internal noise, measurement problem, propagation effects), and are not 
related to any relative source motion phenomenon. (2) On the other hand, if 
the results from the simple moving jet experiment do not agree with the simu- 
lation data, but are found to be in reasonable agreement with the trends 
observed from existing full-scale flight data, then it will become clear that 
the currently observed differences are associated with some relative motion 
effects which have not been accounted for correctly in translating static 
simulation data to the flight results. Subsequent research can therefore be 
directed towards a re-examination of various phenomena associated with source 
motion relative to a fixed observer. In this manner, it will be possible to 
define the most fruitful avenues of approach for reconciling the flight data 
with model simulation results in the long-term research efforts. 


6.2 SECOND CATEGORY RECOMMENDATIONS 

Regardless of the outcome of the fundamental investigation, described 
above, of the noise of a pure jet in flight and in the flight simulation 
facility, it is necessary to continue the basic work of defining the effects 
of forward velocity on jet mixing and shock-associ ated noise, especially for 
heated jets. Preliminary aspects of this work have been completed during 
the current contract; however, considerably more detail must be obtained in 
al 1 aspects of the work recommended below. 

The specific technical objectives of the proposed fundamental study of 
flight effects on jet noise should be: 

(1) to define the limitations of the anechoic free-jet facility tech- 
nique for measuring flight effects on jet noise, using experimental informa- 
tion from current facilities in conjunction with appropriate theory; 

(2) to further develop a methodology for using such facilities, in 
particular, to determine theoretically the proper transformat ion for convert 
i ng facility data to the corresponding in-flight case; 

(3) to measure the appropriate jet flow quantities for the in-flight 
condition, such as mean and turbulence velocity distributions, convection 
speeds and turbulence spectra, for input to a jet noise source alteration 
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model, and to define, by utilizing the data and available theoretical tech- 
niques, jet noise source alterations resulting from aircraft forward motion; 

(4) to extend the acoustic inflight simulation experiments on turbulent 
mixing noise and shock-associated noise, in particular, to examine and quanti- 
fy the effects of forward velocity on heated jets; and 

( 5 ) to calculate the noise field of a jet In flight and compare with 
experimental results. 

The following recommendations are designed to meet the specific 
technical objectives. 

6.2.1 Phase I - Methodology and Corrections for Transformation of 
Free- Jet Data to the Equivalent In-flight Condition 

This phase of the program should be directed toward further development 
of the appropriate corrections and calibrations to convert noise measurements 
from^a free-jet anechoic facility to equivalent linear measurements (parallel 
to direction of flight) for a jet in flight. The effort should use a combi- 
nation of theoretical and experimental approaches as outlined below. 

6. 2. 1.1 Theoretical program. Effects of the free-jet on the sound 
field of the model jet should be studied with the aid of appropriate flow- 
acoustic interaction models. Charts and calculation procedures should be 
produced which relate free-jet facility results to the corresponding in- 
flight noise levels. 

(i) The effect of transmission through the mean flow field of the 
combined primary and secondary jets s.hou1d be calculated by numerically 
solving the Li 1 ley equation. This will provide a standard of comparison for 
the approximations below (ii, iii). 

(ii) The effect of transmission through the outer mean shear layer 
alone should also be studied using the Li 1 ley equation. In this way the 
effect of neglecting the primary jet flow field can be assessed. 

(iii) . The applicability of the geometric acoustics (GA) approximation 
for transmission through the combined (primary + free-jet) jet flow field 
should be assessed by comparing GA results with the full Lilley equation 
results, for realistic mean flow profiles. 

(iv) The above calculations could be performed using a mean flow model 
which neglects spreading of the inner and outer shear layers. Growth of the 
outer free-jet shear layer in the downstream direction may significantly 
affect the interpretation of free-jet facility data. This possibility should 
be assessed using a geometric acoustics model of the acoustic-mean flow 
interaction process. 
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(v) The items listed above refer to interaction of sound with the mean 
flow field in a free-jet facility. The effect of turbulent scattering in 
the outer free-jet shear layer should also be examined in order to determine 
its importance. 

(vi) Based on the results of the above studies, charts and calculation 
procedures should be developed and presented for converting jet noise measure- 
ments obtained in an anechoic free-jet facility to the corresponding in-flight 
s i tuat i on . 

6. 2. 1.2 Experimental program. Experiments should be performed in an 
anechoic free-jet facility to investigate the following aspects of its 
acoustic performance. 

(i) Prior to any other testing, azimuthal symmetry of the jet noise 
field, after transmission through the rectangular free-jet boundary, should 
be experimentally investigated. If the jet noise field proves to be azimuth- 
ally asymmetric, the free-jet exit shape should be modified so as to give 
azimuthal symmetry for the jet noise field. 

(ii) Refraction of sound by the mean flow field of the free-jet should 
be studied, using a single-frequency source placed in the jet flow. Phase 
measurements should be used to determine the orientation of the refracted 
waves. 

(iii) Scattering of sound by turbulence in the free-jet shear layer 
should also be studied in more detail using a single-frequency source placed 
at various locations in the jet. Information in this case can be obtained 
by measuring the frequency spreading effect at different angles outside the 
free-jet. 

(iv) The validity of in-flight jet noise data obtained from the free- 
jet facility should be checked empirically by means of a scaling test, in 
which results are compared from different primary nozzle sizes under the 
same test conditions. This will reveal the limitations imposed by having a 
finite area ratio (secondary/primary jet exit area). 

6.2.2 Phase II - Mean Flow and Source Alteration Effects 

This effort should be primarily experimental but should utilize previ- 
ously developed theoretical results in conjunction with experimental results 
obtained under the current program to develop a source alteration model. 

6.2.2. 1 Experiment. Using a four-channel, two-point laser velocimeter 
system (or other acceptable instrumentation), extensive measurements should 
be made in the jet mixing and fully-developed regions for a cold jet mounted 
in a low-speed wind tunnel. The range of test parameters should be 
0.45<Mj<2 and 0.05 < Vj/Vj < 0.5. In the supersonic experiments both conical 
and convergent-divergent nozzles should be tested so that the effect of shock 
containing flows can be determined. The following quantities should be 
determined. 
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(I) Mean velocity and turbulence intensity distributions for both the 
axial and radial components should be measured throughout the jet. 

(ii) Spectra, cross-spectra and phase speeds (comparable to convection 
speed) should be determined from the two-point measurements for u', v* and 
u* -V* combinations. 

(iii) Fourth-order cross-correlations and cross-spectra should be 
determined for the equivalent source function (vj vj). This will provide all 
the information necessary to substantiate and/or empirically modify the 
source alteration model. 

6. 2. 2. 2 Theory. Theoretical programs are required to provide predic- 
tion of (a) the mean velocity and turbulence intensity distribution, (b) the 
turbulence jet noise source function distribution, and (c) the large-scale 
coherent turbulence structure source function (or more appropriately, the 
noise radiated from the large-scale turbulence structure in the presence of a 
relative velocity). 

6.2.3 Phase til - Determination of Relative Velocity Effects on Jet Noise 

In this phase, the effort should be primarily experimental, but should 
utilize previously developed theory to obtain a method of calculating the 
noise field of a jet in the forward flight condition. 

6.2.3* 1 Experimental program. An extensive set of noise measurements 
for both hot and cold jets should be made as follows: 

(i) In the current exploratory program, cold jet noise measurements 
were taken at forty test conditions which consisted of combinations of five 
free-jet velocities and eight jet exit velocities. The range of this test 
program should be extended, including forward arc measurements if feasible 
in the available facility. 

(ii) Shock associated noise measurements should be made over a range of 
free-jet and jet-exit velocity conditions. 

(iii) The free-jet anechoic facility model jet air supply should be 
capable of delivering heated air up to SOO'^C. In this series of tests, 
approximately 200 combinations of jet velocity, jet temperature and free-jet 
velocity should be tested. The data should be examined in conjunction with 
the cold jet data to determine the appropriate scaling parameters. 

(iv) Limited shock-associated noise measurements should be made for 
the heated jet case to validate the results on temperature effects from 
static jet tests. 

6.2. 3*2 Calculation of forward motion effects on jet noise from source 
alteration data. The information on source alteration effects associated 
with jet mixing in flight obtained from the Phase II program could be used 
to predict forward motion effects on far-field jet noise. For this purpose 
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the source alteration data should be combined with the jet noise radiation 
model, based on 111 ley's equation, which has already been developed for 
static jets under Air Force/DOT contract (Contract F3361 5“C-2032) . 

6.2.3- 3 Comparison of theory and experiment. The experimental jet 
noise data acquired under 6.2. 3-1 should be corrected by the techniques 
developed under 6.2.1 and should be compared with the calculated result of 
6. 2. 3. 2 and with available full-scale in-flight data, as appropriate. 
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APPENDIX 2A 


FLYOVER DATA ACQUISITION CONSIDERATIONS 




The measurement of noise received on the ground from flying aircraft 
is very complicated and expensive. In order to ensure the accuracy of the 
measured data, and the subsequent comparison with static results, It is 
necessary to examine various aspects of data acquisition and correction 
procedures. In particular, detailed consideration should be given to the 
following five aspects of flyover noise tests: 

(1) acquisition and analysis of non-stat lonary data, 

(2) aircraft position and attitude, 

(3) ground reflection effects, 

{h) atmospheric attenuation corrections, and 

(5) real atmosphere effects. 

Furthermore, a study of each of the above aspects should consider (a) the 
procedure to be adopted, (b) the instrumentation requirements, and finally 
(c) the statistical accuracy of the measured data. 

The problems associated with the acquisition and analysis of non~ 
stationary flyover noise data are discussed In this Appendix, 


ACQUISITION AND ANALYSIS OF NON-STAT I ONARY DATA 


The acoustic signal received by a fixed fai — field microphone from a 
static jet exhaust is random, but with stationary properties in the statis- 
tical sense. The accuracy of the measured data can therefore be increased 
to any desired level, simply by increasing the length (or time) of the data 
record. In contrast, in the case of a flying aircraft, the sound emission 
angle seen by a fixed ground-based observer is constantly changing with 
time. The sound pressure signal recorded by the microphone is thus non- 
stationary in nature, and the statistical accuracy of the measurements is 
I imi ted. 

To illustrate this feature, consider the sketch shown in Figure 2A.1, 
which shows the relationship between the sound reception time, t, and the 
reception angle ip (i.e. the source-to-observer angle relative to the down- 
stream jet axis at the time of reception). If the reception time is 
referenced to the time when the aircraft Is directly above the microphone 
(t=0 when <// = 90°) , then this relationship Is 


tan ip = 


H 

VAt » 


(2A-1) 


where H is the normal height or distance of the flight path from the micro- 
phone, and is the aircraft velocity. A sketch showing the variation of 
reception angle (j; with reception time t is also included in Figure 2A.1. 
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Figure 2A.1 Relationship between reception time (t) 
and reception angle (ij;) • 
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From data acquisition considerations, the quantity of interest is the 
rate of change of angle i|». Differentiation of equation (2A-1) with respect 
to time gives 



(H/Va) 

(H/Va) 2 + t2 * 


(2A-2) 


Clearly, for fixed values of H and Va, the rate of change of angle Is 
maximum when t*0. That is, the worst case is presented when the aircraft 
Is directly above the microphone location. 


In the study of inflight effects on jet noise, however, the com- 
parison between static and flight results must be conducted at constant 
emission angles. It is therefore necessary to establish the relationship 
between the emission angle, 6g, and the reception time, t. This can be 
obtained by deriving an expression which gives reception angle ip as a 
function of emission angle 0g. With reference to the geometry between 
these two angles defined In Figure 2A.2, the required result is obtained as 


sine- 

tan Ip * - — ^ . 

Ma + cos0g 


(2A-3) 


Combining equations (2A-1) and (2A-3) , for a pressure-time history recorded 
by a microphone, the time scale is related to the sound emission angle 0g 
by 


t 


PI 

~ (Ma cosec 0g + cot 0-) . 
’'A 


(2A-4) 


In order to obtain the inflight result at a specified emission angle 
6^ from the non-s tat ionary pressure-time history recorded by the microphone 
during a flyover test, the result has to be extracted from the continuous 
record. This is done by dividing the complete record Into several smaller 
records, and the time-dependent pressure fluctuations within each of these 
records are assumed to possess locally stationary statistical properties. 
The length (or time At) of each of these records is clearly governed by the 
angular resolution (tolerable range of angular variation about a fixed 
value of 0g under consideration) specified; as the resolution is increased, 
the available record length reduces. 

Let us consider some typical examples. For an aircraft flying over- 
head at a height of M = 152.4 m and at a Mach number of 0.2 (Va“68.6 m/s) 
simple arithmetic, using the above equations, shows that In the worst 
cases, the record lengths (At) for various angular resolutions (A0g) are 

A6g: 4° (±2°) 6° (±3®) 10° (±5°) 

At: 0.155 sec. 0.25 sec. 0.4 sec. 
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It is clear that for tolerable angular resolution, the typical record 
length is rather small, usually of the order of 1/A second. The accuracy of 
the frequency analysis of such small records is now examined. 

The statistical quality of the data can be evaluated, at first in a 
preliminary manner, by examining the mean square error. For frequency 
analysis in constant percentage bands (for example, one-third octave 
analysis), the mean square error where 


2 * mean square error 
^ true mean square value * 

can be expressed by the well-known result 


e 


2 


BT • 


(2A-5) 


In the above equation, B is the bandwidth of the ideal filter under consid- 
eration, and T is the record length in seconds. For one-third octave 
analysis, the result can be written as 


0.23 T fc ’ 


(2A-6) 


where f^ Is the center frequency of the one- third octave band. It is clear 
that for a fixed record length of time T, the mean square error increases as 
the frequency decreases. For a typical flyover test (T = 0.25 second), the 
mean square percentage errors are 

fc = 50 Hz, _ 35 ^ 

fc = 10 KHz, ~ 0.2%. 

The statistical accuracy can further be expressed In terms of the con- 
fidence limits as determined from standard Chi-square distributions. The 
results can be obtained from any standard textbook on statistical analysis 
of stationary data. 

For samples with n equivalent number of statistical degrees of freedom, 
where n*2BT, the measured mean square value and the true mean square 
value 0 ^ are related by 



(2A-7) 


where the value of (Chi-square) for various confidence intervals may be 
obtained from standard statistical tables. 


For the typical flyover test considered previously, the number of 
degrees of freedom at the lowest and highest one-third octave bands of 
interest are 
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■ 50 Hz, 
fc - 10 KHz, 


B « 11 .5 Hz. 


n -6 


for T = 0.25 second 


B = 2300 Hz, n ^ 1150 ] 

The value of n at low frequencies Is therefore rather small, and examina- 
tion of statistical tables for chi-square distributions shows that the 
confidence level in the measured values of mean-square sound pressure at low 
frequencies is not high. 

The statistical accuracy of the analysis at low frequencies can be 
improved by increasing the number of samples, and computing the arithmetic 
mean after analyzing all samples (record lengths). A description of the 
improvement in the statistical accuracy of the results derived from N 
samples can be obtained from any standard textbook on statistics, and the . 
final result is that the standard deviation is inversely proportional to 
Thus, for large N, the measured (or sample) mean value approaches the true 
mean value. 

In a flyover test, the required number of samples, for some specified 
statistical accuracy, can be obtained during a single aircraft flyover by 
utilizing a line of microphones set in the direction of the flight path. If 
the microphones are sufficiently far apart from each other to record uncor- 
related (and hence necessarily independent) signals, and if each is analyzed 
and compared with the others after suitable time delays, a significant 
improvement in measurement accuracy can be achieved. 
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APPENDIX 2B 

CALCULATIONS OF JET OPERATING CONDITIONS 




In the experiments performed In the free jet facility and the low-speed wind 
tunnel, two factors alter the operating conditions from those used In a 
normal static jet test facility. Firstly, due to aerodynamic considerations, 
the plenum-to-jet area ratio is not large so that there is a finite velocity 
in the plenum or reservoir. Secondly, the jet is exhausting into a non- 
stationary stream so that the total pressure in the surrounding air does not 
equal the static pressure. The following analysis has been used to calculate 
the jet operating conditions. It is assumed that the static pressure in the 
plane of the primary jet exit is constant and, since the jet is unheated, the 
total temperature is also assumed to be constant. A sketch of the typical 
test setup is shown in Figure 2B.1, From the energy equation it is found 
that , 


so that, 


CpTR 



CpTj 







+ 


(y-1) 

2 



(a^-1). 


(2B-1) 


(2B-2) 


where A = Vy/Vj. 

From the assumption of constant total temperature, 


Tj ^ 1 + (y-1) Vj^XV2aT^ 

Tt 1 + (y-1) Mj2/2 

Combining equations 2.B.2 and 2.B.3> we get 

VjVa 2 

2 - d I 

1 + (y-1) Vj2 (x2 - i)/2aj^ 


or 




aj^ 1 + (y-1) Mj2 (1-x2)/2 


(2B-3) 


(2B-4a) 


(2B-4b) 


In order to calculate the velocity in the reservoir, the continuity of mass 
flow condition is used to give 




Letting the area ratio Aj/A|^ = a (a<1) then 


227 




(2B-5) 



Since the total temperature is constant, the jet and reservoir temperatures 
are related by. 


Tj _ 2 + (y-1) Mr^ 

Tr 2 + (y-1) 

so that from (2B-6) and (2B-5) we obtain 

(t^)" 

For a given value of Mj (or Vj/aj from equation (2B-^)), the value of (Tj/Tj^) 
can be calculated from equation (2B-7) hy an iterative technique. 

The static pressure ratio is given by 

Pr/Pt “ (2B-8) 

and the total pressure ratio is given by. 


(2B-6) 


(2B-7) 


Pj 



2 + (y-I)Mr^ 

2 + (yl)x2Vj2/aj2 


Y/(Y-1) 


(2B-9) 


Consider the following example. We need to find the static and total pres~ 
sure ratios for the following conditions, 

X = .25, ot = .25 and Mj * .9. 

From equation (2B-4b) (Vj/aj) = .838157- 

Iteration of equation (2B~7) using Newton-Raphson two times gives 
“ •8633^, so that = .11*518 from equation (2B-5) 

Thus, finally: Pr^Pt “ 1*667921 

and P(^/Pj = 1.641591 * 


Clearly the required jet operating conditions can be set using a measurement 
of either total or static pressure, whichever is experimentally more 
convenient. 


229 




APPENDIX 2C 

TEST CONDITIONS AND MEASURED 
(UNCORRECTED) TURBULENT MIXING NOISE DATA 




Run 

Number 


7 

8 

9 

10 

11 

12 

13 

1i» 

15 

16 

17 

18 

19 

20 
21 

22 

23 

Zk 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


TEST CONDITIONS 
M = 1 . 0 , CONVERGENT NOZZLE 




\l= —] 


hl-Li] 

PR 

ao 

3 o 



Pt 

0.400 

0.019 

0.047 

0.407 

0.969 

1.1143 

0.400 

0.051 

0.127 

0.406 

0.970 

1.1139 

0.400 

0.100 

0.251 

0.406 

0.970 

1.1139 

0.400 

0.150 

0.376 

0.405 

0.970 

1.1135 

0.400 

0.200 

0.500 

0.404 

0.970 

1.1128 

0.500 

0.025 

0.049 

0.513 

0.952 

1 .1870 

0.500 

0.051 

0.101 

0.513 

0.952 

1.1867 

0.500 

0.100 

0.200 

0.513 

0.952 

1 .1865 

0.500 

0.150 

0.300 

0.512 

0.952 

1.1860 

0.500 

0.200 

0.399 

0.511 

0.952 

1.1855 

0.600 

0.030 

0.050 

0.623 

0.930 

1 .2858 

0.600 

0.051 

0.085 

0.623 

0.931 

1.2851 

0.600 

0.100 

0.167 

0.622 

0.931 

1 .2847 

0.600 

0.150 

0.250 

0.621 

0.931 

1 .2837 

0.600 

0.200 

0.333 

0.620 

0.931 

1.2823 

0.700 

0.036 

0.052 

0.738 

0.905 

1 .4176 

0.700 

0.051 

0.072 

0.737 

0.905 

1.4173 

0.700 

0.100 

0.143 

0.737 

0.905 

1.4163 

0.700 

0.150 

0.215 

0.736 

0.905 

1 .4151 

0.700 

0.200 

0.285 

0.734 

0.905 

1 .4134 

0.801 

0.042 

0.052 

0.858 

0.875 

1.5951 

0.800 

0.051 

0.063 

0.857 

0.875 

1.5934 

0.800 

0.100 

0.125 

0.856 

0.875 

1.5921 

0.800 

0.150 

0.188 

0.855 

0.875 

1.5903 

0.800 

0.200 

0.250 

0.853 

0.876 

1.5873 

0.901 

0.048 

0.054 

0.985 

0.840 

1.8344 

0.900 

0.051 

0.056 

0.984 

0.841 

1 .8313 

0.900 

0.100 

0.111 

0.983 

0.841 

1 .8291 

0.900 

0.150 

0.167 

0.981 

0.841 

1.8259 

0.900 

0.200 

0.222 

0.980 

0.842 

1 .8216 
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M = l.i*, CONVERGENT-DIVERGENT NOZZLE 


Run 

Number 

^0 


II 

«j 

(■?) 

LlLLl) 

pr 

Pt 

37 

1.17^ 

0.063 

0.05^ 

1 

.380 

0.725 

3.0610 

38 

1.175 

0.100 

0.085 

1 , 

.381 

0.725 

3.0640 

39 

1.176 

0.150 

0.128 

1 , 

.380 

0.725 

3.0611 

i»0 

1.178 

0.200 

0.170 

1 

.380 

0.725 

3.0596 




M = 1.7, 

CONVERGENT-I 

DIVERGENT NOZZLE 



Run 

Number 


Vt 

^0 


) " 4 =^) 

LlLLl) 

Pr 

Pt 

42 

1.343 

0.077 

0.057 

1 .683 

0.638 

4.7725 

43 

1 .344 

0.100 

0.075 

1 .682 

0.638 

4.7709 

44 

1.346 

0.150 

0.112 

1 .683 

0.638 

4.7764 

45 

1 .348 

0.200 

0.148 

1.683 

0.638 

4.7760 
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APPENDIX 3A 

GEOMETRICAL ACOUSTICS FOR PARALLEL, SHEARED FLOW 




The notes in this Appendix describe the derivation of analytic, geomet- 
ric acoustics results for the radiation from various types of acoustic source 
distribution surrounded by a parallel sheared flow. 

The notation in this Appendix and Appendix 3B differs in some respects 
from that used elsewhere in the report. The main difference is that here the 
external flow speed or flight speed is 'Vq' as opposed to 'Vj' or 'V/\' used 
elsewhere; in general, external flow conditions are referred to here by sub- 
script 0 . This means that 9 q here is equivalent to 6j in, for example. 
Appendix 3C and not 6q. Notation in general is explained with the aid of 
Figure 3A.1. 


3A.1 SOURCE IN SHEAR FLOW - GA LIMIT 
(FIXED FRAME ANALYSIS) 


The stationary-fluid result, in terms of the (k, oj) cross power spectral 
density (cpsd) of the acoustic source distribution q(x, t) , would be (with 

Vs = 0) 

2 2 

l^r Is ~ 7 Ps ^ (3A-1) 

where P(u>) is the psd of the fai field pressure in the source-region fluid, 

0 ) is the fixed-frame frequency, and the functional form of the source (k, w) 
cpsd is Q(k, cj) =p| <!> (k, 0 )) . 

Provided q contains no 3/3t operators, the Doppler transformation gives 
the corresponding moving-fluid result “-in terms of fixed-frame frequencies^ 
as 

Rp P(<^)ls = J Ps (“s“s'^^s» (3A-2) 

where 

= Ds = {1 + gig • Vg/ag} ^ 

= {1 + Vgcos0g/as}"^ (3A-3) 

and WsOfg/ag is the radiation wavenumber in the source- reg ion fluid. In 
equation (3A-2) and elsewhere, denotes the distance travelled by the 
wavefronts in the wavenormal direction. 

To go from the source region to far outside the jet, we assume that the 
external radiation pattern is axi symmetric. It then follows from energy 
conservation that 
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WAVENORMAL 
UNIT VECTOR a 



RECEPTION 
V Rp/a POSITION 



SOURCE POSITION (FIXED) 


Figure 3A. 1 Definition sketch. 


(3A-if) 


^ PoC 
{Rr PM}s PsDs** 

where 


®o - {1 + 2o • = {1 + VoCoseo/ao}"^ = oi^/oi (3A-5) 

is the Doppler factor too/u) relating frequencies relative to the external 
flow to frequencies in the fixed frame. 

Combining equations (3A-2) and (3A-4) gives the radiation psd in the 
external flow as 

{Rp P(u) }q “ Po Ps Y ®(<‘>s3s/^s* ®o ®s^ ■ (3A~6) 


3A.1.1 Compact Volume Acceleration Source at Rest 


2 

If we put Q(k, 0)) * Pg ^(k, w) and assume $ compact, we get 

{R? PU)>o Po Ps |-®(0. 0.) Ds-2 


compared with 


{Rr Po(o))}o = Po^ |-$(0, 0)) 


(3A-7) 

(3A-8) 


if the entire flow field were replaced by fluid (pq, a^,) at rest. 

Thus for the same emission direction and distance, the effect of the 
j^Zow fieZd on the radiation from a specified compact volume acceleration 
distribution at rest is given by 


PM 



(stationary point volume- 
acceleration source). 


(3A-9) 


This result is not to be confused with any effects which arise from having 
a convected eouroe pattern; only the fZuid is moving here. 


3A.1.2 Compact Volume-Acceleration Quadrupole at Rest 
If q = pg (a^S j j/3x| 3Xj) , so that 
Q(k. 0)) - pg^k. kj kj^ k^ 


(3A-10) 
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is the functional form of the source (k, w ) cpsd, then (3A-6) is replaced by 

Rf P{«j) 1 q “ ~2 “si “sj “sji “sm ^i j tm(<*>s2s^®s* 

^s 

X u>5‘* Dq** Ds"^> 

and if $ is compact, we get (using u>g = wOg) 

2 P S 7T i \ 

Rp P(“)lo ■ iT 2 ’ 

^s 

X 0)*^ Do** Ds2 . (3A-11) 

For a given emission direction ag and distance Rj., the effect of the flow 
field on the radiation from a specified Sjj(x, 1) compact distribution is 


P((u) ^ 

Po\ao/ 


(stationary point volume- 
acceleration quadrupole)* 


(3A-12) 


3A.1.3 Compact Vol ume-Di splacement Quadrupole at Rest 

If q = PgCa^/Bt^) (9^Bj j/3xj 8Xj) relative to the flow, so that 

Q(k, 0 )) “ Ps^ k| kj k^ kj^ (k, to) 

is the functional form of the source (k, to) cpsd, then the compact-source 
result is (in the GA limit) 

P P 

P(“)io "" ■■^i*“si “sj “s£ “sm «’ij£m 
^s 

X (U® Oq** Dg® . 

Thus, the effect of the flow field, for a given Bjj(x, t) distribution, is 
given by 


y = — ( — ) (stationary point volume 

Po \^o/ displacement quadrupole), 


if 25 and are held constant. 


278 



3A.2 NON-COMPACTNESS EFFECTS AND 
SOURCE CONVECTION 


The simplest estimate of non-compactness effects In aerodynamic noise 
IS obtained by assuming Isotropy of $ (or with respect to the direc- 

tion of the vector argument (kj Vei, kz Ve 2 , kg Vga, - kg V_i). (For 
convenience a convection velocity (Vci, 0 , 0 ) Is assumed, i.e. parallel to 
the mean flow velocity). Thus, we put 

Q(k, m) = ps 2 $(k, 0 )') (o)' = u-kiVj-i 

= con vec ted f name 
frequency) ; 

♦ (k, 0)') = $(0, {(kiVei)2 + (k2Vg2)^ + (kgVeg)^ + o)'2}^) 

(isotropic function assumption) . (3A-13) 

Putting in (3A“13) gives the required value of $ as 

4>((0sgis/as, u - 0)5051 V^-j/ag) = $(0, o) J , say, (3A-14) 

which defines the modified frequency o)^ and the modified Doppler factor 
Note that (for kj matching). 

Thus, the stationary compact-source result (3A-11) is modified to 


Rr 


P(o)) 


- 7T 


(0 


2 ^sj o^sm ^\jim ^ 




(non-compact convected volume-acceleration quadrupole) 


Also, 


(3A-15a) 


P(w) lo 


PoPs IT , 

"TIT 2 ^“sl « 

^s 


(non-compact convected 

Note that in (3A-15), Is a 

density factors. 


sj asm “m)“m®> ~ 


8, Do** D56 


D 8 
m 

vol ume-dlspiacement quadrupole) . (3A-15b) 

kinematic quantity, i.e. it contains no 
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3A.2.1 Recovery of Lighthill Radiation Model 

In the Lighthill radiation model, the sources radiate directly into the 
surrounding fluid (density p^, sound speed, Bq, velocity yj . The foregoing 
results may be converted to the Lighthill radiation model by simply putting 


Ps Po 

as -> 3o 

Vs ^ Vq (3A-16) 

«s ^ “o 

and Ds *^o 

Further details are given below, in sections 3A.2.2 and 3A.2.3- 

3A.2.2 Evaluation of D^, (GA Limit) 

From (3A-13) and (3A-14), 

» 1,2 „ 2 ^ + fc 2 

+ [u) - (wpapi/ao) Vcl^ ' 

where 

|<2 + kl “ (ws/a^)^ ~ ki 

= (u)/as)2 Ds - (wo“ol/3o)^ 

= (w/as)^ D 5 - (w/ao)^ Dq^ cos^0o 

(“01 ^ cos0q) . (3A-18) 

Thus if we make the simplifying assumption that Vg2 = Ve3 = Vgt> say, the 
modified Doppler factor follows from (3A-17) and (3A-18) as 



The following special cases of (3A-19) are of interest. 

(a) Vq ~ ® (source region at rest relative to external fluid): then 

Dq = 1 . Also, 


'e3 


(3A-17) 
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(3A-20) 


-1 

Dg = (1 + VgCosSj/ag) = (1 - VgCos0Q/ao) 
from consideration of phase speeds in the axial direction. 

(b) Lighthill radiation model: all s subscripts in (3A-19) are 

changed to 0. Thus, 




i. 


Note that 


^ COS0Q = 1 COS0O 


(3A-21) 


(3A-22) 


from the definition of Dq in equation (3A-5) . 

(c) General GA modified Doppler factor with V^j = If the convec- 

tion velocity is set equal to the local flow velocity in the shear layer, 
equation (3A-19) simplifies to 


^m ^o 


cos20, 




V 


ei vetl . .^2 / *et . , 

: J + Dg^rel I — + 1 


^o / \ as^ 

where the relative motion Doppler factor is defined by 


^s,rel 


= 1 


(^) 


COS0, 


(3A-23) 


(3A-24) 


V - V 

1 + COS0g 

ac s 


0 ) n 

o “o 


-1 


(3A-25) 

(3A-26) 


Equations (3A-25) and (3A-26) follow from the phase-speed relationship 
across the shear layer, Vg + as/cos0g = Vq + ao/cos0Q. 


3A.2.3 Check on the Ffowcs Williams (ref. 7) result obtained using 
Lighthill radiation model 

Changing all s subscripts to 0 in equation (3A-15) gives the radiation 
from a Lighthi 1 1 -type convected quadrupole distributon in a completely 
uniform moving medium (velocity Vq parallel to the convection direction). 
Thus, in proportional bandwidth form. 
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(3A-27) 


{Rr P(oj)a)}o 


Pq r ^(0/n 

{o^oi ^oj ^ol ^om j 


X 


D 


6 

o 


D ^ 
‘-'m 


To convert from source-region frequency u) to the frequency o)q relative to 
the uniform flow, no change in the proportional-bandwidth mean square 
pressure is involved, but simply a frequency shift by a factor ajo/o) = DQ, as 
given by equation (3A-5) • The frequency gJq corresponds to the frequency 
heard by an observer on the ground, as the source region is moved past with 
velocity -Vq (assuming the atmosphere is at rest). 

The factor Dq^/D^^ follows from (3A-21) and (3A-22) above as 


(1 + 




(3A-28) 


which is the same as the factor obtained by Ffowcs Williams (ref. 7) for a 
homogeneous source model with a Gaussian space-time correlation function. 
The Ffowcs Williams form of source distribution is a special case of the 
more general form assumed in (3A-13)» which was first proposed by Crighton 

(ref- 111). 
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APPENDIX 3B 


GENERAL RESULTS FOR PARALLEL SHEAR FLOWS 
(SOURCE DISTRIBUTIONS WITH AXIAL COHERENCE ONLY) 




These notes describe the framework In which the LI 1 ley equation flow 
factor solution can be used to estimate flow-acoustic alteration effects due 
to forward motion. 

The notation used here and in Appendix 3A differs In some respects from 
that used elsewhere in the report. The main differences are described in 
the introduction to Appendix 3A. 


3B.1 the modified DOPPLER FACTOR 


When only axial source coherence effects are retained, the modified 
frequency, (o^,, defined in section 2 of Appendix 3A, reduces to 


I .e. 


»^-(^cos6o V^,)" + (. 


0) ^ COS0„ 

®o ° 


= 0) |( 1 - Dr 


2 .,2 




= wDn,, say 


In (3B.1), 


Do - (i + ^ cose„)-‘ . ^ 


(3B-1) 


(3B-2) 


3B.2 EQUIVALENT STATIONARY POINT VOLUME-ACCELERATION 
OR VOLUME-DISPLACEMENT QUADRUPOLE 

The quadrupole cross-power spectral density in axial wavenumber and 
frequency, evaluated at the actual frequency wT and the actual wavenumber 
ki * (ui/aQ) D qCos 6 q required for radiation at 6 q to the xj axis, is assumed to 
be the same as the cross-power spectral density at the modified frequency 
Urn and at zero wavenumber, ki =0. Coherence of the quadrupole distribution 
in the transverse plane (x2» X3) is assumed to be negligible 

Thus, if P(<»)) is the fai field pressure psd’^ radiated by a stationary 
point volume acceleration/displacement quadrupole of the aame instantaneous 
total strength, the radiation at the same distance from the actual non- 
compact quadrupole distribution is given by 


^RelAitive to coordinates fixed in the source region^ i.e. fixed 
with respect to the nozzle, 

^Azimuthally averaged. 
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P( 0 )) = Hi^) (v = 1) 

or 

p(u) = P(o)m) Dm“8 (v = 3) (3B-3) 

The Dp,"** or D^“® factor arises from the u** or oj® frequency dependence of 
P(a>). If proportional frequency-bandwidth values are required, the factor 


becomes or thus, 



P(o))u = P(tOni)com Dm~^ 

(v = 1) 


or 



P(a))w = P((A)n,)wn, D„,“® 

It 

(3B-A) 


which is a more convenient form for scaling purposes. 


3B.3 SCALING OF RADIATED PRESSURE 
POWER SPECTRAL DENSITY 

In a uniform fluid at rest, the radiation fran the stationary point 
quadrupole described above would be given by 

Di \ ^ ^ m 2 v +3 (3B~5)^ 

Rr PU)u) = 2 7^ “oi “oj “o£ 


Here, for example, 


(Vj^ d3)2 

(0,cu) 0) oc 

I j£m 


(3B“6) 


since cross-power spectral density (in oi) of the total 

integrated quadrupole strength, /vj'vj'dx or /Cj'Cj'dx, and it Is assumed 
in this example that vj* scales on Vj or that 5;' scales on Vj(Vj/d) ^=d, 
the nozzle diameter. Also, 


so 


(A) ^ V j/ d ; 

Rp P (q)) q) 
d"Po" ao" 


(S) 


(3B-7) 


(3B-8) 


^Note <1>; j£^(0,o)) is the li) -spectral density of the quadrupole strength. 
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Thus if F is the factor by which the point volume-acceleration/displace- 
ment quadrupole radiation intensity is multiplied, due both to finite Mach 
numbers in the flow and to departures of (p , a) from (p„, a_) , the scaling 
law derived from (3B-1*) and (3B-8) is 


Rp P(a))(ji) 


d^p 




o o 


:(v) 


U) 


■(2v+3) 




(for given w^d/Vj) 


(3B-9) 


3B.4 EFFECT OF EXTERNAL STREAM/ OF SPEED V^/ ON 
FLOW RADIATION FACTOR F(u) 

The flow factor refers to the situation (a) sketched in Figure 

3B.1 in which a stationary point quadrupole is placed in a shear layer with 
(in general) a non-zero external velocity. 

If the flow factor in the situation (b) sketched in Figure 3B.1 is de- 
noted by F^gi' (w) , where the flow velocities are reduced everywhere by Vq to 
bring the external stream to rest, then f(^) ( u) is given by 

F^''^((d) = (o)Do) (for same olq^) (3B-10) 

since (a) is equivalent to having a convected point quadrupole of frequency 
u) moving upstream at Vq as in (c) . The fixed-frame frequency corresponding 
to^ (c) is since (c) is simply (a) with a velocity -Vq superimposed; 

this is the reason for having wDq as the argument of F^Nj) . 

The factor in (3B-10) represents the basic effect of a com- 

pletely uniform flow of velocity ^ on the point quadrupole radiation+; once 
this has been allowed for the ^***®o^ factor accounts for departures 

from this uniform flow as described above. 

Note that the combination D^"^2 v+ 3) „hich appears in the 

radiation scaling law when (3B-10) is substituted into (3B-9) , may be 
written as 


Do(VDo) 


-(2v+3) 


D 


-(2v+3) 

m,re1 


(3B-11) 


where is the modified Doppler factor as in equation (3B-1), but with 

'^ci •‘splaced by V^j “Vq. Thus, the scaling law may be written in the more 
useful form 


As follows from the Doppler trans formation » 
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QUADRUPOLE ^ 

N) 





Figure 3B.1 Flow factor —velocity profile situations 



P(ti))o) 
d2 Po2 ao** 


(S)' 




„2v+3 

®m,rel 


(3B-12) 
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APPENDIX 3C 


A GEOMETRIC ACOUSTICS INVESTIGATION INTO THE 
INFLUENCE OF FREE-JET MEAN VELOCITY AXIAL VARIATIONS 
(NON-STRATI FIED FLOW) ON SOUND REFRACTION BY THE 
SHEAR LAYER (REUTIVE TO THE STRATIFIED FLOW CASE) . 




The simple, analytic relationships used elsewhere in this report, be- 
tween ray and/or wavenormal angles inside and outside an axi symmetric 
cylindrical sheared flow field are results which are valid only for the 
special case of a parallel or unidirectional, stratified*** mean flow field. 
These terms are used occasionally to mean that such properties persist at all 
locations upstream and downstream of the local flow region under considera- 
tion. Here that unrealistic flow field model would be referred to as an 
infinite stratified, parallel flow field model and is thereby distinguished 
from the more realistic (locally) parallel, stratified flow field model. The 
latter is a fairly realistic model of the initial mixing region of an axi- 
symmetric jet flow, since the mean velocity vector varies in direction by 
only 2 or 3^ from the axial direction and its magnitude is only a weak 
function of the axial coordinate compared with its strong dependence on the 
transverse or radial coordinate 

Here the restriction that the flow be stratified is removed so that the 
magnitude of the axial mean velocity vector is allowed to vary with axial 
position in a realistic way. With that dependence specified, the refraction 
or bending of sound rays as they propagate through a parallel, but non- 
stratified, shear layer is examined and emerging ray angle results are 
compared with those calculated for the parallel stratified shear layer case; 
ray displacement effects are also studied. 

In a sense this non-strat i f i ed flow model may be more unrealistic than 
the stratified one; for example, if the flow is incompressible, then the 
continuity equation can be written as 


3V 


r 


3x 


1 3(r Vt) 
r 3r r 3<fr 


(3C-1) 


where are the mean velocity components in the radial (r) and azimuthal 

((p) directions. Thus, when Vx=Vj((x,r), the continuity equation requires, in 
general, a non-zero radial mean velocity V^, i.e. the flow is non-parallel 
which is in conflict with the flow model as specified. In principle then, 
when the stratified flow restriction is removed, non-ax iai mean velocity 
components should be included in the mean flow model, in particular (for 
non-swirling flows) the component = Vt;(x,r) with V^ = 0. This would be a 
useful extension of the present work. 


3C.1 A JUSTIFICATION FOR UTILIZATION OF THE GEOMETRIC 
ACOUSTICS EQUATIONS TO DESCRIBE SOUND PROPAGATION 
THROUGH THE FREE-JET SHEAR LAYER 

According to Morse and Ingard (ref. 112) the geometric approximation is 
"...appropriate when a and V vary slowly and 'smoothly' (i.e. when the 


'^Defined as the flow properties being a function of the transverse coordinate 
only. 
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wavelength times the Laplacian of a and V are negligible compared to their 
gradients) The geometric acoustics equations used below are obtained 

by assuming a solution for the acoustic pressure, p’, of the form 

p‘(x,t) = A(x,cj) exp[“jk^(x) + joit] (3C-2) 


to the convected wave equation 


a2 Dt2 


v2p‘ 


(3C-3) 


(D/Dt E d/dt + Vi 8/axj) . 

$ is the so-called eikonal function which defines the surfaces of equal 
phases and A(x) is the amplitude function. The eikonal equation 

|grad<J>| = (3C-A) 

a + y • n 

(n = gradO/|grad<l> I ) (3C“5) 

is derived by substituting the solution (3C-2) into equation (3C”3) and neg- 
lecting all but the highest powers of k (=u)/aQ) - However apart from the 
trivial case of a completely uniform flow, V(x) = constant, equation (3C“3) is 
I tsel f an approx i mat i on : i f V = (V^(r) ,0,0), a=a(r), p=p(r), then the exact 

equation is simply the left hand side of the Li 1 ley equation set equal to 
zero (since sources within the free-jet shear layer are not considered here, 
only sound propagation through it), that is, 

v2n'! +i^L 

Dt2 p dr Dt 


9£l 

3r 


+ 2 




3x3r 


(3C-6) 


In this particular kind of flow field, an infinite, parallel, strati- 
fied flow, the geometric approximation is valid, in the first instance, when 
the last two terms can be neglected in equation (3C-6) thereby reducing it 
to a form closely resembling that of equation (3C-3); clearly, the geometric 
acoustic solutions of equation (3C-3) also satisfy this reduced form of the 
Li 1 ley equation. 

By considering the solution of equation (3C-6) at a particular frequency 
and for a certain radiation angle, it Is not difficult to show that the mean 
velocity gradient term is of order l/S^j relative to the two leading terms, 
for a given mean velocity profile shape. The Strouhal number is based 
upon the vorticity thickness of the free-jet shear layer, 6^j, and the 
centerline or exit mean velocity, Vj. 

In the free-jet configuration <5^jwill be typically an order of magni- 
tude larger than that in the primary jet for a given "ray path" at 90° or 
in the rear arc. The ratio of the tunnel or free-jet exit velocity, Vj, to 
the primary jet centerline velocity varies over a wide range of values (for 
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a given ray path) but is less than unity so that 

S(5T > 10 Sgp 

where S^p is the Strouhal number based on primary jet quantities. It can be 
evaluated by linking it with the axial location of maximum radiation in the 
primary jet and is found to be (empirically) of order unity. Hence, in the 
present problem, the gradient term in the homogeneous Lilley equation is at 
least an order of magnitude smaller than the two leading terms that form the 
basis for the geometric approximation. 

The Lilley equation solutions presented in section 3-2 for sound radia- 
tion from the primary jet alone do exhibit significant deviations from the 
geometric acoustic solutions but that represents, in principle, a combina- 
tion of source-flow interactions and propagation effects. That is, sources 
are located within the shear layer and the sheared flow modifies the 
emission process as well as that of propagation. In any case the corres- 
ponding free-jet Strouhal numbers would be an order of magnitude larger than 
the values adopted there and the trends in those results indicates a rapid 
approach to geometric acoustics with increasing Strouhal number. 


3C.2 BASIC EQUATIONS FOR A GEOMETRIC ACOUSTICS STUDY 
OF REFRACTION IN A PARALLEL NON-STRATIFIED MEAN 
FLOW MODEL OF THE FREE-JET 


The equations governing the propagation and refraction of acoustic rays 
in an otherwise time independent mean flow field can be written as three 
second order, ordinary differential equations [Ugincius (ref. 113)]: 


d / dr \ dr 

^ ^ (3C-7) 


Later Ugincius (ref. 1l4) corrected these to 
d / \ / \ dr 

'(I?- ^7 - <3c-8) 

where r = r(s) is the position vector of the ray as a function of arc length 
s and 

V 

W = y j (3C-9) 


y 


a + V 


(3C-10) 


N = y Vp/a 


(3C-11) 
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* 2|J|i| „s6. (3C-12) 

y E ej a/3x| and V‘ = ej3/axj where 

xj s= dr/ds. (3C-13) 

The speed of sound a(r) and velocity vector V(r)are arbitrary, well-behaved 
functions of the position vector, r; the speed of sound a^ is a reference 
value. The relationship between V, n, a, and 9 is indicated in Figure 
3C.1. The vector n is the unit normal to the wavefront (a surface of 
constant phase); it is the direction of the sound ray only if there is no 
flow, y E 0. With flow the ray direction is given by the vector Vj- where 

Vr - an + y (3C-14) 

i.e. the direction of the ray, dr/ds, is given by 

dr/ds = Vr = Vp/Vr . 

The angle between the wavefront normal direction and that of the mean flow 
velocity is designated 0 and that between the ray direction and the mean flow 
direction, 

Equations (3C“7) and (3C-8) differ by the term 



which Ugincius (ref. 1l4) argues is negligible for low Mach number flows. The 
present investigation is based on equation (3C-7) > since typical Mach numbers 
of the jet flows of interest are less than 0.25. However, the influence of 
that extra term should be considered in future work [by solving equation 
(3C-8)]. 


Equation (3C“7) simplifies considerably when the mean velocity vector is 
defined to be every where parallel to the x axis and axisymmetric, i.e. 



y ={Vj^(x,r) , 0, 0} 

(3C-l6a) 

with 

a = a(x,r) 

(3C-I6b) 

so that 

W = {yV^/a, 0, 0} 

(3C-17) 


^0 

^ * a + Vj^cose 

(3C-18) 


= 1 + (V^/a)^ + 2(Vx/a) cos6 . 

(3C-19) 
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Vp = + COS0, a sin9, 0 


(3C-20) 


dr 

ds 


+ a cose 


a sine ^ 

Vr ’ 


(3C-21) 


With equations (3C-16) through (3C-21) the x component of equation (3C-7) 
reduces to 


d cos 6 I 

ds a + V^cose I 


1 

Vp(a + V^cose) 


ax 


+ cose 



(3C-22) 


When a stratified flow field is assumed [Vx = Vx(r), a = a(r)], equation 
(3C-22) can be solved immediately to give 


cos 6 

a + V Jocose 


constant. 


It follows that if (i) e=0j, Vx^'Vj, a=aj and (ii) 0=0o» ^^ = 0, a = at 
two points on the ray path then 

^o 

— 1 — + Vt = — - — 

COS0J * COS0O 

a relation that also can be obtained by '^matching axial phase velocities" on 
both sides of a stratified shear layer (a vortex sheet being a special case). 

Equation (3C-22) was adopted by Schubert (ref. 115) for his "Numerical 
Study of Sound Refraction by a Jet Flow," based on ray or geometric acoustics 

With the following relations 


d ^ , a ^ , 3 


\J^+ a COS0 a COS0 

COSiJ^ = rj ; slni|; = — 

''r *r 


Schubert (ref. II 5 ) reduces equation (3C-22) to 


d6 

ds 


1 I • n 
t:- sine 


9a 


- 3a 

- COS0 T— + COS0 




3x 3r 

and then by differentiating 

tamji = a sin0/(V;^ + a cosO) 


3Vx 

ne _ - cose 


1^)1 


(3C-24) 

(3C-25) 


(3C-26) 


(3C-27) 


to obtain 
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(3C-28) 


^ “ u 2 


(1 + Vx/a*cose) ~ - sine ^ (V^/a) 


Schubert (ref. 115) finally obtains the following equation for the ray 
curvature dij//ds 

* ■ I? 'ir 


where M = Vx/a. The ray speed can be expressed in terms of 

Vp/a “ Mcosi|( + (1 -M^sin2^)^ 


(3C-29) 


(3C-30) 


so that the right-hand side of equation (3C-29) Is a function of the mean 
flow properties and is only an algebraic function of the angle between the 
ray and mean velocity vectors, ip. 

In order to solve equation (3C-29) numerically, the following identities 
are substituted 


dijj 

ds 


simj^ 


- d^x/dr^ 

[l + (dx/dr)^]^/^ 

1 

[1 + (dx/dr)2]i 


COS^I) 


dx/dr 

[1 + (dx/dr)2]^ 


(cotip = dx/dr) 


(3C-31) 

(3C-32) 

(3C-33) 


SO that equation (3C-29) is a nonlinear, second-order, ordinary differential 
equation of the form 


d2x 

dr^ 



X, 




(3C-3^) 


Equation (3C-3^) is reduced to two first order equations, in the usual way, 
by denoting dx/dr by Yj: 


X, r) 

^ (3C-35) 
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The mean flow model described below is isothermal {a = constant = B q) and the 
mean velocity gradients are specified with respect to coordinates non- 
dimensional i zed by rj — an effective free-jet nozzle radius. In these 
coordinates the equations to be solved are 

dY, 3V 3V 

Fi (Yi, X. R) F 2 (X, R) (3C-36a) 

§ = Yi (3C-36b) 


where X = x/rj, R = r/rj, 



The trigonometric functions have been retained here, although they are 
eva 1 uated wi th 

sinijj = [1 + Yj2] ^ ( 3 C- 38 ) 

cosil^ = Yi • sim|;; (3C-39) 

the ray speed, Vp, is evaluated with equation (3C-30) . 


3C.3 MEAN FLOW FIELD MODEL 


The free-jet is assumed to be isothermal and ax i symmetric with an 
effective nozzle radius, rj. The axial mean velocity profile, Vj^, is given 
by 

Vx = 7 {1 - erf (x)} (3C-40) 


where 



(3C-i*l) 


The constant by =.297 is chosen to yield a value of 0.663 at r ■ rj; that 
is, at the "lip-line" the mean velocity is chosen to be 66.3% of the center- 
line value. A nominal value of the spreading parameter, a, is 13.5. Rewriting 
equation (3C-4l) in terms of non dimensional coordinates 
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X 


(3C-42) 


o(R-l) . 
X *"T 


the required radial and axial derivatives with respect to R and X can both 
be expressed in terms of di?j^/dx: 


3Vv 3Vv 




3R 

~ dx 



3V, 




3X 

dx 

from 

equation (3C-40) 





1 



dx 


equation (3C-42). 





0 



3R 

‘ X * 

3X 

mean 

velocity gradi 

ents 

can be i 



a 

“i “X 


_ 

3R 

X 

IT e 


_ 

(R- 

Jj_ d^x 


3X 

X 

dR 


1 -X 
e 


(X+ bT) d 


(3C-43) 

(3C-i»i») 


(3C-45) 


(3C-i*6) 


dR 


(3C-47) 

(3C-48) 


3C.4 ACOUSTIC MODEL AND SOLUTION METHOD 


The flow-acoustic model is sketched In Figure 3C.2. A point source is 
located on the free-jet centerline, at an axial position, Xg (X has its 
origin at the free-jet nozzle exit). This idealized point source represen- 
tation of the turbulence noise source distribution within the primary jet 
is discussed in Appendix 3D. An acoustic ray radiated from the point source 
intercepts the inner edge of the free-jet shear layer at a point (Xj, Rj). 
The angle between the ray and the downstream jet axis, is related to its 
wavenormal angle, 6j, by 


cot 6j 


cos 0j + Vj/Bq 
sin 0j 



(3C-1»9) 
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The ray path through the parallel but non-stratif led shear layer is deter- 
mined from a numerical solution of equations (3C-36) with the initial 
condi tions 


Yl(Rl) = = cotipj 

and 

X(Ri) = Xi 


(3C-50) 


(3C-51) 


Numerical integration of equations ( 3 C- 36 ) from R = Rj to R = R 2 (the outer 
edge of the shear layer) then gives the ray angle i|)q or 0 q as it enters the 
ambient medium 


YaCRa) cotij^Q = cot0Q 


(3C-52) 


and its axial position 


X(R2) = X 2 . (3C-53) 

The initial and final values R = Rj, R 2 are determined from values of X that 
define the inner and outer edges of the shear layer, i.e. 


(1 - Vj^) < AV 

when X < Xi 


Vx < AV 

when X > X 2 

(3C-54) 

where AV is a small number, e.g. 

1 

0 



3C.5 SUMMARY OF RESULTS: COMPARISON OF NUMERICAL 

SOLUTIONS TO THE GEOMETRIC ACOUSTIC EQUATIONS 
WITH STRATIFIED FLOW RESULTS 

Numerical results are presented for four axial source locations 
X 5 = 0 . 6 , 1 . 1 , 1.6 and 2.1 for six ray paths corresponding to the wavenormal 
angles (within the uniform flow) 67 = 30°, , 60°, 75°, 90° and 105°. The 

numerical solutions are compared with the analytic, stratified flow results. 

The isothermal, ax i symmetric mean flow model of the free- jet is 
characterized by the jet exit velocity, Vj, or the Mach number, Vj/aQ, and 
the spreading parameter, a. In view of the Mach number limitation in the 
basic equations [according to Ug i nc i us (ref . Ilif)], the value of Mj/ao in the 
following results is limited to values less than or equal to 0 . 2 , which is 
the highest value of immediate interest in the present investigation. Further 
work is required to determine the significance of this Mach number limita- 
tion, but it is not expected to have any serious bearing on the Mach 0.2 
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results given below. The spreading parameter is not varied here but is held 
constant and equal to a nominal value of 13-5. Preliminary results, not in- 
cluded here, indicate that variations in the value of this parameter (which 
encompass effective values for the free-jet) have little or no effect on the 
numerical solutions. 

The solution method was partially verified by solving equations (3C*36) 
with axial variations suppressed and the axial gradient term omitted; the 
numerical solutions for say, Gq, should agree with values calculated from the 
equation based on a stratified flow model, that is, from equation (3C-23) 

0o = cot"^ {cosGj [(aj/a^ + Vj/aQ cosGj)^ - cos^Gy] (3C-55) 

or since aj = a^. 

Go = cot ^ {cosGj [(1 +Mj cosGj)^ - cos^Gj] ^}. (3C-5&) 

A set of Gq results are shown for comparison in Table 3C.I; the numerical 
and analytic results agree to four significant figures. For reference 
purposes values are included in Table 3C.I. 


Gt° 

Gq° (numerical) 

Gq° (analytic) 


105 

105.84 

105.84 

93.48 

90 

90.00 

90.00 

78.69 

75 

75.75 

75.75 

64.59 

60 

62.96 

62.96 

51 .05 

45 

51.72 

51.72 

37.94 

30 

42.42 

42.42 

25.13 

Table 3C.I. 

Verification of Numeri 

cal Solution Method: 

Compar i son 


of stratified flow Gq values given by numerical solu- 
tion with values given by equation (3C-56); Vj/ao = 0.2, 


a = 13.5. 

When equations (3C-36) are solved wvth axial variations and the axial 
gradient of is included, the resulting G^ values for this non-stratif ied 
case (N.S.) differ by less than from the corresponding stratified flow 
(S.) values, as shown in Table 3C.M. 

With reference to equations (3C-36) and (3C-37) the superficial reasons 
for the almost negligible influence of axial variations in are that (i) 
the axial gradient term is of order V-f/ao smaller than the (leading) shear 
gradient term and (ii) it changes sign at the lip-line, making a positive 
contribution to d(cotiJ^)/dR in the first half of the interval Rj <R<1 and a 
negative contribution in the remaining Interval 1 <R<R 2* These observa- 
tions do not concern axial variations in Vx which are present in the 
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0T° 

00° (N. 

105 

105.76 

90 

89.99 

75 

75.82 

60 

63.12 

45 

52.00 

30 

42.88 


0 

0 

CO 

A0o° 

105.84 

-.08 

90.00 

-.01 

75.75 

.07 

62.96 

.16 

51.72 

.28 

42.42 

.46 


Table 3C.II. Comparison of Non-Strati fied and Stratified Flow 
0o Values: Vj/bq = 0.2, a = 13. 5. 


iTftl'll Tm non-stratified flow 

effects at Mach numbers lower than Vj/a^ = 0.2 can be also safely ignored. 

In the results so far, no mention has been made of the influence of Xc 
the point source or ray origin position on the jet centerline. Its value 
determines, with the ray angle the point at which the ray enters the 

shear layer (Xi, Rj) . However, it can be shown that this has no influence 

on the magnitude of the refraction (Oq - iPj) . as follows. 

Equation (3C-36) may also be expressed in terms of x (in place of R) : 


d(cotij^) . dVy 

= F(cot^. X) 


(3C-57a) 


dX _ (cotij;) X 

dx a 


(3C-57b) 


where 


F(coti|/, x) = Fj(coti|^,x) + FoCx) 


(X+bj) I 


(3C-58) 


Hence the equations are decoupled; that is, equation (3C-57a) could be inte- 
grated separately between x = Xi and X2, given the initial value cot^;, to 

nm J e^’^'-ging ray angle, 0^, at X = X 2 ; since the 

limits X-Xi, X 2 define the inner and outer edges of the shear layer at anv 
axial location, 6^ - ,/,-p is independent of X^. The ray displacment informa- 
tion IS contained in the solution to the second equation. 

r.„c demonstrates that the refraction or bending of acoustic 

nim H H ^ r ’T'" Of axial location in the present flL 

opon the initial ray angle, the profile shape, , 
expressed in terms of the similarity coordinate, x, and Vj/a^. 

It follows that results such as those given in Table 3C. I I are valid for 
any axial source position. 
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Finally, solutions for the axial displaoement of the refracted ray, 
X 2 -X 1 , (see Figure 3C.2) are considered in the following way. The ray paths 
wlthln'the uniform flow and ambient medium are extrapolated without change in 
direction to the "lip-line," R=l, as shown, for example in Figure 3C.2, the 
axial intercepts being denoted respectively by Xg,, X^. The numerical results 
for Xu, Xn indicate that, to a very good approximation, these points are 
coincident. A set of results is given in Table 3C.I1I for the four source 
positions. These numerical results are consistent with the analytical result 
indicated by equation (3C-57b) that the ratio X 2 /X 1 is independent of axial 
location and, in particular, of the source position X 5 . 

To summarize, numerical solutions to the geometric acoustic equations^ 
indicate that to a good approximation two simple rules govern sound refrac- 
tion by a parallel, non-s trat i f i ed free-jet mixing region. 


Rule (a) The sound ray refraction (or angle change) can be 

calculated to a good approximation with the analytic 
formulae for refraction by stratified flow, and 

Rule (b) the ray path (or axial displacement) can be calculated 
as if the stratified flow refraction all takes place 
abruptly at the "lip-line." 



Xs = O' 

.6 



Xs = 1.1 


0T 


Xu 

AX 

Xjl 

Xu 

AX 

105 

0.5^ 

0.54 

.00 

1 .04 

1 .04 

.00 

90 

0.80 

0.80 

.00 

1.30 

1.31 

.01 

75 

1.08 

1 .08 

.00 

1.58 

1.58 

.00 

60 

1 .1*1 

1 .42 

.01 

1.91 

1 .92 

.01 

^5 

1.88 

1 .90 

.02 

2.38 

2.41 

.03 

30 

2.73 

2.77 

.04 

3.23 

3.28 

.05 


Xs = 1 

.6 



CM 

II 

X 


0T 


Xu 

AX 


Xu 

AX 

105 

1 .5^ 

1.55 

.01 

2.04 

2.05 

.01 

90 

1.80 

1.81 

.01 

2.30 

2.31 

.01 

75 

2.80 

2.09 

.01 

2.58 

2.59 

.01 

60 

2. 41 

2.43 

.02 

2.91 

2.93 

.02 

‘*5 

2.88 

2.91 

.03 

3.38 

3.42 

.04 

30 

3.73 

3.79 

.06 

4.23 

4.30 

.07 


Table 3C. 1 I 1. Axial displacement of acoustic rays: "lip-line" intercepts 

Xji, Xu of extrapolated ray paths for axial source positions 
Xs = 0 . 6 , 1 . 1 , 1.6 and 2.1. (Vj/ao = 0.2, o= 13-5). 
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Rule (a) is not an unexpected result; it justifies the utilization of 
results from Appendices 3A and 3B which are based upon a parallel stratified 
flow for free-Jet model calculations, tt is also an encouraging indication 
that uaye acoue tic models based upon infinite, parallel stratified flow 
models, such as Lilley^s equation, may be good approximations, particularly 
as they are often utilized for conditions approaching the GA limit. Finally 
Ru e (a) justifies the use of the vortex sheet model for sound refraction 
calculation since there is no distinction between the model and one consist- 
ing of a finite thickness stratified shear layer in this immediate context. 


Rule (b) is consistent with and provides the justification for models 
developed and used elsewhere (ref. 1l6 and 117) in which the free- jet shear 
layer is replaced by a vortex sheet. However, in those models both the ray 
paths and shear layer transmission effects on the sound amplitude are calcu- 
lated on the basis that the shear layer is replaced by a vortex sheet. Rule 
(b; only confirms that ray paths can be calculated on that basis. Transmis- 
sion effects on the sound amplitude are analyzed in Appendix 3E. 


The application of these two rules to the calculation of the angle e™ in 
the present free-jet configuration (e„, is the polar angle based upon 
an origin at Xg = 0.6) yields results which are in very close 
agreement with those calculated with axial variations and gradients of \?x 
included in the basic GA equations. In Figures 3C.3 through 3C.6, the ray 
paths outside the free-jet shear layer are shown for the four source posi- 
tions; these have been derived from the full numerical solution to the axial 
variation/gradient equations. Inside the shear layer rather than show the 
exact ray paths, the rays have been extrapolated to the "lip-line" to illus- 
trate the magnitude of the error involved if Rule (b) is applied. Also, in 
Figures 3C.3 through 3C.6, values of 6^ are shown; values from a calculation 
based on rules (a) and (b) give results that agree with those to the nearest 
degree. Exact figures are given in Table 3C.IV. 
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Figure 3C.3 Ray paths on the X-R plane for wavenormal (= emission) angles 
0t = 3O° ( 15 °) 105 °, with resultant microphone angles, 0^. 

VT/ao = 0.2, Xs=0.6 (xs'/d = 0). Code; as in Figure 3C.2. 







Figure 3C.5 Ray paths on the X-R plane for wavenormal {= emission) angles 
At = ^0° (1S°) 105°, with resultant microphone angles, 0^,. 

\/I/a. = 0 2. X. = 1.6 (x '/d = 10). Code: as in Figure 3C.2 




Figure 3^.6 Ray paths on the X“R plane for wavenormal (= emission) angles 
Of = 30^ ( 15 °) 105°» with resultant microphone angles, 9^. 
Vj/ao - 0.2, Xs = 2.1 (xs'/d = 15). Code as in Figure 3C.2. 
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o O 

®mab 


105 

90 

75 

60 

45 
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103.47 
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87.88 

73.47 
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60.17 

60.15 

47.63 

47.60 

34.93 
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98.32 

98.44 

82.48 

82.55 

68.28 

68.32 

55.38 

55.40 

43.37 

43.42 

31.20 

31.35 


0T 



®mab 



105 

93 

90 

77 

75 

63 

60 

50 

45 

39 

30 

27 


11 

93.26 

07 

77.16 

03 

63.09 

54 

50.60 

07 

39.18 

43 

27.68 


87.79 

87.97 

71.53 

71.65 

57.66 

57.76 

45.60 

45.71 

34.69 

34.87 

23.59 

23.94 


Table 3C.IV, Comparison of values of 0^ calculated from numerical 
solutions to axial variation/gradient GA equations 
with values 0mab calculated by Rules (a) and (b) : 
Vj/ao = 0.2, a = 13.5. 
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APPENDIX 3D 

A MODEL FOR THE EFFECTIVE POINT SOURCE - AXIAL 
LOCATIONS OF PRIMARY JET TURBULENCE MIXING NOISE 
AS A FUNCTION OF STROUHAL NUMBER 




11 o^ recent measurements by Grosche (ref. 118) and Laufer et al (ref 
1131 provide quantitative information on apparent or effective axial distri- 
butions of the turbulence mixing noise source in subsonic and supersonic iets 
as a function of the Strouhal number, fd/Vj. 

The utilization of this information is necessary in the present probiem 
because measurements in Fiight Simulation Facilities are usually taken rela- 
tively close to the jet, as close as 20 jet nozzle diameters (see Figure 
30.1), while flight data is usually acquired at much larger distances of the 
order of a hundred jet nozzle diameters or more. Suppose it is assumed that 
radiation at a oarticular Strouhal number originates mainly from a region 10 
diameters down: tream of the jet nozzle and the radiation level at a particu- 
lar polar angl is under investigation. Then with the origin at the jet 
nozz e, lig^'^ <ata (i.e. Rjn>100d)at angle 6^ wi 1 1 be a good measure of the 
radiation from the source region at angle 0 ^,. If the distance is reduced to 
Rm = o0 diameters, then in the example shown in Figure 30.1, the measured data 
would correspond to radiation at an angle 0^ + 6®, while at a distance of 
Rm-20 diameters it corresponds to Qm + 30°. Furthermore, the ray path length, 
the distance between the measurement point and the effective point source 
location, differs, in general, from the "measured radius" R„,. Clearly a 
correction procedure is required if "near-field" flight simulation data is to 
be extrapolated and compared with true far-field flight data. 

Morfey (ref. 120) has proposed a simple model in which the real axial 
source distribution for each Strouhal number is replaced by a point source 
located at the axial position Xs‘+ where the apparent distribution reaches 
a maximum. An empirical law for the variation of the point source axial 
position, Xj'/d, with Strouhal number is derived from the Grosche and Laufer 
(ref. 118, 119) data on the following basis. The radiation at a particular 
frequency is associated with a shear layer thickness and a centerline 
velocity. In the initial mixing region the centerline velocity is a constant 
and equal to the jet exit velocity, Vj, whereas the shear layer thickness is 
proportional to the distance from the jet exit. Consequently, the empirical 
law for the initial mixing region is of the form 

Xs'/d a j. 

In the fully developed region the reciprocal of the centerline velocity is 
also proportional to the distance from the jet exit Ut large distances) and 
hence the law takes the form 


Xs'/d a -y , 

The Strouhal number, S, is now replaced by the modified Strouhal number (to 
allow for the convected nature of the source), S„,, which Is defined below. 


prime is used to distinguish, this x coordinate, which has its origin at 
the primary jet exit ptane^ from that used elsewhere within Chapter 3 which 
has its origin at the free-jet exit plane* 
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Morfey (ref. 120) obtains the constants of proportionality from the Grosche 
and Laufer data (ref. 1l8, 119) and finds that, to a good approximation, 


Xg'/d 



for the initial mixing region and 


Xg'/d 


for the fully developed region. In Figure 3D. 2 those empirical expressions 
are compared with some of Laufer's (ref. 119) data. The agreement is parti- 
cularly good in the important low Strouhal number range, .15sSni<.6, 
(leading to the largest angle corrections) whereas the discrepancies at high 
Strouhal number are of little significance for the present application. The 
"cross-over" point between the two empirical laws should be, according to 
the physical basis for these laws (ref. 120), in the region of Xg'/d- i*-5. 
Here it is chosen to be the point Sm = 1.5, i.e. x^/d » 3.3 to achieve a 
smooth "transition." 

The modified Strouhal number is defined In this exercise as 



(^m ~ S D^) . 

^m “ Mj.cos0j)2 + 

(Mj. = 0.67 Vj/aQ, a = 0.3) 

Therefore, given d, f and Vj in consistent units, and 0j, the axial source 
location, xs'/d, is calculated via the following steps: 

S = fd/Vj 

“ {(1 -McCOS0j)2 + 

(Me = 0.67 Vj/a^, o = 0.3) 

= S Dm 

Sn, > 1.5: XgVd = 5/Sm 
Sm < 1.5 Xs'/d = VSmi 

In the Lockheed free-jet facility the measurement radius, R_, is 5l*d, 
where d is the primary nozzle diameter (d = 2") or 5.k rj where rj is the 
effective free-jet nozzle radius (to a good approximation rj - 20") . The 
primary jet nozzle protrusion beyond the free-jet exit plane, x„, is 6 d or 
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0.6 rj. The effect of the extended ax!al source distribution on the values 
of the measurement polar angle, 0^, which correspond to specified polar 
radiation angles, 0j, is illustrated in Figures 30.3 through 3D. 6 for the 
case of zero free-jet velocity, Vj. That is, Xg has its origin at the free- 
jet exit plane so that the first value of Xg corresponds to a source 
location at the primary jet nozzle, succeeding values correspond to 5d 
increments. The first value (Figure 3D. 3) represents the ideal case where 
the measured and actual polar radiation angles 0jjj, 0j and distances R^, Rp 
coincide. A source location of 5d downstream (Figure 3D.^») gives rise to 
angle differences of up to 5°, of 10 d to 11® and of 15 d to 16°. The 
effect of forward motion on the co-flowing stream (the uniform potential core 
of the free-jet) is not included In this source location model. A prelimi- 
oary study has indicated that this effect is unlikely to be important in the 
present correction procedure. 
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Figure 3D. 3 Ray paths on the X-R plane for wavenormal (= emission) angles 
0x = 30° (15°) 105°. with resultant microphone angles, 6^. 
Vx/a^ = 0.0, Xe “ 0.6 (xc'/d = 0). Code: Ray path. 







ce: 
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Figure 3D. 5 Ray paths on the X-R plane for wavenormal (= emission) angles 
0J = 30° ( 15 °) 105°, with resultant microphone angles, 0j^. 
VT/a„ = 0.0, Xc = 1.6 (xo'/d'IO). Code: -»■ Ray path. 




APPENDIX 3E 

FREE-JET-TO-IWT TRANSMISSION AMPLITUDE CORRECTIONS 
FROM A STRATIFIED FLOW - GEOMETRICAL ACOUSTICS MODEL 




The main result of the investigation described in Appendix 3C is that 
(in the geometric acoustics limit) the direction and position of rays emerg- 
ing from a parallel non-strati fied model of the free-jet shear layer can be 
accurately calculated with the stratified flow analytic expressions together 
with the condition that the rays are refracted abruptly at the lip-line. 

The result is illustrated by Figure for ray path calculation purposes 

the realistic situation depicted in (a) can be replaced by the oonceptuat 
model shown in (b) . 

Clearly the flow model is now a cylindrical region of uniform flow 
bounded by a vortex sheet. However, emphasis on this aspect of the result 
has been deliberately avoided since in one other important respect the 
vortex sheet feature is not utilized. In the analysis given below, the 
finite width stratified shear layer is retained for the purposes of 
deriving a relation between the amplitude of the direct radiation from 
the primary jet inside the free-jet uniform flow region and the radiation 
amplitude outside the free-jet in the ambient medium. 

The basis for this relation is the law of energy conservation in a ray 
tube. This law applies to the geometric acoustics limit; in the same limit 
the energy flux associated with a particular ray in the uniform flow region 
is completely determined by the amplitude of the direct radiation, since, by 
definition there are no rejections, single or multiple, fvom the shear 
layer. In a vortex sheet model [for example, see Amiet (ref. Il6)], there 
are reflections and the corresponding relation between the transmitted and 
incident pressure amplitudes at the vortex sheet, in effect, replaces the 
energy conservation law that will be used here. While the results from each 
method may not differ drastically, the arguments presented in Appendix 3C in 
support of a GA model apply equally well here and therefore this approach 
will be used in preference to that associated with a vortex sheet flow model. 

The principles employed and the end result are identical to those 
described by Schubert (ref. 115) see his Appendix B: “Ray acoustics for 

nonspreading Jet“ ~ although with the results of Appendix 3C it will be 
argued that the result is also valid for the real, spreading, free jet. 

In the geometric acoustics limit the law of energy conservation in a 
ray tube holds, however complex the mean flow field within the free-jet 
shear layer. On the other hand the analytic relations between angles and 
distances used below rest upon the stratified flow assumption. Conceivably, 
a more realistic flow model than that investigated in Appendix 3C might lead 
to results that are significantly different from those given by the strati- 
fied flow expressions. For the present, this possibility is set aside, 
pending further work in that direction. 

With reference to sketch (b) in Figure 3E.1, the law of energy 
conservation along a ray tube, given by Blokhintsev (ref. 121) is 


I dA = constant 


(3E-1) 



SOURCE 


bd = 

be - Rra 
fg = RrT 

LIP-LINE 
(R“l or r=rj) 



e' APPARENT SOURCE LOCATION 


Figure 3E.1 Free-jet shear layer can be replaced by a vortex sheet 
at the "lip-line," R = 1, for ray tracing purposes. Ray 
tube cross-section changes from dAj^ to dAo across 
shear layer; wavenormal angle from to 0o. 






where I is acoustic intensity and dA the cross section area element of the 
ray tube or conical sector (the direction of the intensity, I, is identical 
to the ray direction). The subscript “T" denotes evaluation at some point 
along the ray tube in the uniform free-jet flow and subscript "’I" denotes 
the conical sector (ipj, + di(^y) ; subscript “ 0 " denotes, as usual, evalua- 
tion in the ambient medium. It is preferable to rewrite the equation as 


•jt *T 6 “ *o ^*^0 


(3E-3) 


where Ijo is the intensity component normal to the wavefront; dAjQ is the 
elemental wavefront cross section area element formed by the intersection 
of the ray tube and the wavefront. 

Morfey (ref. 122) has expressed Blokhi ntsev ' s (ref. 121) acoustic 
intensity definition in terms of the wavenormal component, Ip,: 

In = + (Vn/pa2) ^ + (Vn/a)2 +p\l„Z^ (3E-4) 


where overbar denotes time average, p, a are the local density and speed 
of sound and u,^, Vp denote the components of the fluctuating and mean 
velocity normal to the wavefront. With the GA relation 


P* = P a Up 


the expression reduces to 




(3E-5) 


(3E-6) 


Thus, in the present problem the intensities l-p 0 , Iq are given by 
I- 


p't 


'Te 


Pt aj 


{1 + (Vj/aj) cos0j}2, 


(3E-7) 



The corresponding area elements (see figure 3E-1b) are 

dA -|-0 = (2irRpj s i n0 j) Rcj d0-j-, 
dAo = (27TRro sin0Q) d0^. 


(3E-8) 


(3E-9) 

(3E-10) 


where the distance from the apparent source location, is derived as 

f o 1 1 ows . 
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Writing Rpa Figure 3E.1b) 

R^_ = be + ce 

= (bd - cd) + ce 

= (Rpo " cosec0Q) + ce (3E-11) 

and noting that 

dX. d 

ce = - sin0o ^ = “ sinO^ (cotipj) (3E-12) 


where 
R PQ cen 

is the axial location of the point c, an analytic expression for 
be found if the stratified flow relation 


— { ( V j/a j) ( aQ/a j) + 3 ^cos0q}/kj 

(3E-13) 


(b 2 = 1 - (VT-/aT)2) 


is used 
derived 

to determine the derivative d(cot\{)j)/d0Q. Equation (3E“13) 
from the definition 

has been 


Vj + aj cosej 
sin0^ 

(3E-U) 

and the 

stratified flow relations (from "phase velocity matching") 



COS0J = (aj/a^) cos0q/D-p 

(3E-15) 

and 

sin0-|- = (a^/a^) kj/Dj 

(3E-16) 

where 

Dt = 1 - (Vj/ao) cos0„ = , ^ 

(3E-17) 


and 


Kj = {(aQ/a-]-)^ Dj 2 - cos^Oq}^ . 


Differentiating equation (3E~13) with respect to 0 q gives 


d cot •J'T 

del 


= - sin0Q {aJa^)^/K^' 


(3E-18) 


(3E-19) 
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so that from equation (3E-12) 


ce = sin^e^ (aQ/a-|-)^/Ky3 


(3E-20) 


and hence 


*^ra = Rro+ cosecS^ {(sin6Q/Kj)3 (ao/aj)2 - 1 }. (3E-21) 

Substituting the intensity expressions, given by equations (3E~7) and 
(3E-8), into the energy conservation iaw, equation (3E-3) , leads to an 
expression relating the pressure amplitudes pj and p^* : 


py2/pi2= (pt aT)/(Pp ap) dAp 

{1 + (Vj/aj) COS0 j}^ dAj 0 

From equations (3E~9) and (3E-10) the last factor can be written as 


(3E-22) 


*^^o _ *^ro ^ra d (cos9q) ^ Rj.^ # , , 2 

dAje „ 2 d(cosej) I 2 ’ 

’'rT ^rT 


(3E-23) 


equations (3E-15) and ( 3 E-I 7 ) have been used to evaluate d(cos0_)/d(cos0T) . 
Combining equations (3E-22) and (3E-23) the final result is 

Pt'^/Po'^ = (pj ^T^'/Po) (^ro Rra/R??)* (3E-2i») 


The isothermal forms (pj = Pq, aj = ap) of equations (3E-21) and (3E-24) are 
used in Appendix 3F to correct free- jet measured data to that which would be 
measured under IWT conditions. Ther e the distance RrJ is replaced by R , a 
specified emission radius" at which p-pi 2 is required. This may take any 
value, that is, an estimate of pj ' 2 can be made at any distance, since it 
simply varies as RTf-y under IWT conditions [as indicated by equation (3E-24)]. 
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APPENDIX 3F 


FREE-JET-TO-IWT DATA CORRECTION PROCEDURE 




The free-jet SPL measurements are taken at a number of evenly spaced 
polar angles S^ii* ' 2 ... on an arc of radius centered at the primary 

jet nozzle (see Figure 3F.1). The primary objective of the correction pro- 
cedure outlined below is to convert free-jet data to that which would be 
measured under IWT conditions for any value of the emission radius, Rf, and 
at specified emission angles 6 ej, i = 1 , 2 ... . 

I 

As a secondary objective, corrections are included to allow for the 
finite distance between the effective center of the turbulent mixing noise 
source distribution and the nozzle exit. This minimizes the errors in- 
volved in the extrapolation of small R„, flight simulation data to a large 
R,^ or R^ IWT conditions. 


PROCEDURE OUTLINE: FREE-JET TO IWT CONDITIONS 

Step 

1. Given 0ei=0T» '^t/^q. calculate 0 q: 

cosiJ;t + Vr/ao 

cotiiT = L— — 2. 

sinGj 

cot 0 Q = cos 0-|-/{(1 +Vy/aQ cosQj)^ - cos^Qj}^ . 

2. Given f*, d* , Vj‘, Mq and x, calculate Xg'/d with procedure given in 
Appendix 3D. 

3- Given R^/d, rj/d, Xp/d calculate Xg , and hence R^o* ®m' 

Xn = (xn/d)/(ry/d) 

Xs = X + (xs'/d)/ (rj/d) 

Xjj^ = Xg + COtl|jj 

Rro - ^^rrf ~ Jl^sin^G^}^ - £ cos0q 
cose^ = (Rn? + - R?o)/2)IR„ 


*Cons'Letent units. 
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Figure 3F.1 Definition sketch for f ree-jet* tO” IWT data correction procedure. 


where 




5. 


5 

"> (rj/d) 

i = (Xj^ - X^) - cot0Q. 


Interpolate for SPLI(0^) at frequency f; 

SPLl(0n,) = (1 -h) SPLM(0^j) + h SPLM(0^j^^) 


where 


®mi - ®m - ®mi+l 


h = 


^mi 

^mi+i ®mi 


Given R^/d calculate and apply amplitude correction given in 
Appendix 3E: 


*^ra “ f^ro cosec0Q { (s i nO^/x-j-) 3 - 1} 
SPLC(Rr, 0-p) = SPLl(0n,) + lOlogjo Cp Cr , 


where 

<T = Hi - V-p/a^ cos0q) 2 - cos^e^}^ 

Cp = Dt" 

" r+ (Vt/Io)cos 0^ = ’ ‘ ^V^o> ‘'os^o 
Cr = Rro Rra/Rr^ 

Rr = (Rr/d)/(ry/d) . 
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APPENDIX NA 

TABULATED VELOCITY MEASUREMENTS 
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APPENDIX AB 


MEAN VELOCITY DECAY ANALYSIS 




The analysis of this section is based on that by Kleinstein (ref. 123) 
who obtained solutions for mixing the turbulent axisymmetric free jets using 
the "modified Oseen method." In this method the conservation equations are 
linearized in the plane of the von MIses variables resulting in a system of 
equations of the heat-conduction type. 


The equations of continuity, momentum and energy are, neglecting the 
axial pressure gradient (though the experiments to be considered were per- 
formed in a wind tunnel, the axial pressure gradient was small since the 
ratio of tunnel area to jet exit area was large): 


9(pVir) a(pV2r) 

9x 3r 


(4B-1) 


av 3Vj 


i _L 

r 3r 


(pe) r 


ar 


and 


,, aH ^ an i a 


i£ll r 


ar 


(4B-2) 

(^tB-B) 


Equations (4B-2) and (^B-3) have an identical form, as would the equation 
for conservation of chemical species In a nonreacting turbulent fluid. This 
general form may be written. 


PV 


3Pl 

1 


ar 


i _L 

r ar 


A|^(pe)r 


3^ 

ar 


A dimensionless stream function is Introduced in the form 
3r ’ax 


(ifB-i|) 


('♦B-Sa) 


where ip = (f/2)^, (^B-5b) 

and overbars indicate dimensionless variables. 


The von Mises transformation with independent variables T and x is 
applied to equation ( 4 B-I 4 ) giving. 




(AB-6) 


The term [A|^(pe)p Vi?^] is written in a series expansion in about a 
point X on the axis yielding: 


[Ak(pe)p Vir^]- « [A|^(pe)]. + Oiw"*) . 

A j T X 9 U 


(AB-7) 
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Neglecting terms of fourth order and higher in this expression and 
substituting into equation (4B-6) yields 


fit * i -L 

* 'i' 3>F 

Finally, equation (4B-8) may be written 

fit - i JL Lflt 

f 3y ^ 3y 

where. 


3Pr 


[Akrai„,o’' 


(4B-8) 


(4B-9) 


X 

Cr “ 2 I [ArPeI^ Q dx . (4B-10) 

o 

Equation (4B“9) is the well-known heat conduction equation and solutions to 
this equation are available. We note that in the form of equation (4B-9) a 
constant shift of the dependent variables does not effect the solution in 
the ^ plane. Thus, a step function boundary condition in velocity 
equal to the velocity difference for r<1 and zero for r<1 is permissible. 
The boundary conditions on P|^ are, P|<^(H',o) =g|^(y) and Pl(®» 5|<) *0 with 
gRlY) by. 


9kW 


PkJ - Pko 0<y<yj 
0 ^ 


Subject to these conditions, the solution to equation (4B-9) is equivalent 
to the heat conduction equation solution for an instantaneous cylindrical 
source of heat at t=0, Carslaw and Jaeger (ref. 124), and is 



and along the axis 

PRa = 1 - exp (- 'J’/ACr) (4B-11b) 

where the subscript 'a' refers to jet centerline conditions. 
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Clearly, from equation (4B-11b) the variation of with R Is given by 
?k(x) - - ln[1 -Pjca (x)l, (AB-12) 

and experimentally obtained centerline distributions of the function P can 
be used to determine 

We will seek a form for ^ which will be regarded as a function of x 
and X only* To this end ^ is assumed to take the form 

^ r .5 [(p Vi)^ - (p Vj)^] (4B-13) 

in a manner analogous to the formulation due to Ferri, Libby and Zakkay 
(ref. 90 ) for Vj*0. r ,5 is the half-width of the jet, given by the point 

at which the momentum Is given by 

P Vi pu - J5 [(p Vi)a + (p Vi)^]. (4B-li») 

In order to obtain the form of ^ the behavior of the variables as 
Cv"*’® is examined. 

The solution, P(4', 5^) » given by equation (i»B-11a) is conveniently 
divided Into an axial variation and a radiation variation viz. 

P(y, 5k) * Pa['l'j/(2?k)^;0lP* T/(25,^)^]. 

Now as 5 |^ ->■<» 

Pa['‘'j/(25k)^.0] - (i»B-15a) 

and P* ['i'j/(25k)^I 'i'/(25|^)^] » exp (- (i|B-15b) 

Also the density at large axial distances becomes uniform so that from 
equation (4B-1A) 


Vi('l'.5. ?v) - ^ tv- + V-] 


or 


P. P*5 

a • ^ 


= Pa» 


thus P*['J'j/( 25 y)^, Y 5 /( 25 ^)^] - h 

The form of r .5 as is thus given by 


(4B-16) 
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'.5 


.5 


- 2 

r.5 = 


¥ df 1 


P-V, 


f df 


Po } 
o 


PaP*(Vj - Vo) +Vo 


For convenience in this case all quantities are nond i mens i ona 1 i zed with re- 
spect to jet exit conditions, so that Vj = 1 , P j = 1 , etc. and more importantly, 

= i. 

Denoting Vq = Vq/Vj by X the expression for r _5 becomes: 


_ 2 

^5 


Po 


«p d*!* 


exp(- ’f^AcJd -X) +X 


x?7 [(1-X) + 25wX] 


and finally using equations (4B-15b) and (4B-16) this becomes: 


_ 2 [(1-x) + HyX] 

■■•s “IT [(1-X) + 25^X] 


o - ■ ’V 

Using L'Hopital's rule, it can be shown that 

Lim (r ,- ) = — 


X->*o 


(1<B-17a) 


(4B-17b) 


which is the expression obtained by Kleinstein (ref. 123) for X = 0. Thus, 
the expression for equation (I 4 B- 1 3 ). becomes , 


■pe = kjPjj 



.(1-x) + 1»5yX 

.11^ 

— In 

^Po 

|.(1-X) + 2?^X 

2Sv 


as X -»-0, pe = kiP 


— l^iPo ^ (1 
i 


1 I (1-x) + 

2^ (1-X) + 2K^\ 


(4B-l8a) 

(1»B-l8b) 


The variation of Cy with x can be found from equation 12. The values 
of Cy 3 function of x are tabulated in Tables 4B-I and 4B-II and is 
plotted as a function of X in Figures liB-1 and 4B-2. The straight lines in 
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Table 4B . I Variation of with veloci ty ratio for Mj = . 47 . 
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Table 4B . I I Variation of 5y with Mach number for X * . 1 . 
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Figure kB.] Variation of Cy with axial distance Mj = .it7. 
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Figure kB-} , which are seen to fit the data weli, are given by the equation, 

Cv = fMx - .35, (4B-19) 

where f(x) is the slope of the straight lines for each value of A. The 
slope of f(A) Is plotted as a function of A In figure 4B-3, for Mj = .1*7. A 
good fit to the measurements is given by 


f(x) = .0399 (1 - .9223A). (l»B-20) 

It is Interesting to note that the value of f(A) for A = 1, that is equal 
jet exit and moving stream velocities, is non-zero. This implies that a 
mixing process is taking place for equal velocities and the turbulent shear 
stresses produced give rise to a finite eddy viscosity. This mixing is due 
to the boundary layers which exist on the inner and outer surfaces of the 
jet nozzle upstream of the j'et exit. If no such boundary layers existed 
and the jet lip was of zero thickness, then the form of equations (l|B-l8) 
and (AB-20) would suggest that f(A) be given by 


f(A)/f(0) = (l-A). (lfB-21) 

It is also of interest to note that In spite of the dependence of the term 
under the square root sign In equation (AB- 18 ) on both A, p^, and x the 
eddy viscosity ?5r is independent of x for the developed region of the jet 

fl^. However, the eddy viscosity is clearly a function of velocity ratio 
and density ratio. 


Since the exact value of the constants in the expression for f(A) will 
have a direct dependence on the jet and free stream initial profiles, a 
more general form is given by equation (4B-21) and this form will be used. 
Making use of this expression in equation (4B-16) leads to: 


(5a - i) I , 


(4B-22) 


where K is a function of Pq and A and, from equation (^B-19), X^ = .70. This 


to that given by Wltze (ref. 3) for A=0 except that no 

For the case of A = 0 
and was able to derive 
ij and PT. Since the current 
measurements only considered three cases at a fixed velocity ratio for 
different Mach numbej^, there is insufficient data to provide the variation 
of K with A, Mj and p^. The values of K calculated from the present 
measurements are listed in Table 4B.III, 


’ ^rer 

dependence has been given in equation (l»B-22) . 

Wltze correlated a large amount of experimental data 
expressions for the dependence of < on M , and 


The value of k for Mj = .li7 in Table 4B.III is taken from the analysis by 
Wltze (ref. 3) for A-0 and Mj - .i»7. Using this set of values for K the 
measured centerline velocities were plotted as a function of k(1-x)x In 
Figure kB-^. The general agreement is fairly good since the experiments 
were carried out in a wind tunnel with a small, but finite, pressure 
gradient and there existed boundary layers at the jet exit. 
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Figure 4B.3 Variation of dgv/dx with velocity ratio, Mj*.47. 
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Mj 
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Po 
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.096, .206, .305, .^13, .'•97 
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1.37 

GO 
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1.378 

.01»83 

1.67 

.098 
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.Oi»39 


Table i»B.III Variation of decay coefficient x with A and Mj. 
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(Y - L)/(Y - ®A) 
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Figure 4B.4 Jet centerline velocity decay. 



APPENDIX AC 

JET NOISE IN FORWARD FLIGHT - SOURCE 
ALTERATION EFFECTS 




introduction: lighthill's model for 

JET NOISE SOURCE SCALING 


The basic source of jet noise is taken to be the volume-acceleration 
quadrupole density V|Vj. In the absence of the flow-acoustic interaction 
U.e. for radiation in a uniform medium of density po and sound speed c^) , 
the acoustic intensity per unit volume of turbulent flow is predicted 
ILighthill (ref. 110) equation 13] to scale according to 

I ''e 4 

'voi “ ^ • (^c-1) 

e 


Here Vg is the correlation volume for v'jvj fluctuations; Xg is a character- 
istic time scale, and the covariance of v'jvj has been assumed proportional 
to the fourth power of the rms turbulent velocity q (= <v|^^i). Source 
non-compactness and convection effects are disregarded for the time being- 
they wi 1 1 be discussed later, when the source alteration effects implied bv 
equation (4C-l) have been studied. 


QUADRUPOLE CORRELATION VOLUME AND TIME SCALE 


quantities Vg and Xg are assumed to scale on the local shear layer 
thickness 6 and the local rms velocity fluctuation q. Thus 


X0 oc 6/q 

and hence from equation (1|C-1) 

•voi 


('♦C-2) 

(^C-3) 


[hZ-k] 


Alternatively, the intensity radiated per unit area of shear layer is 
pendent of 6 (and hence streamwise position), and proportional simply 
*Imax* 


i nde- 
to 


INTENSITY RADIATED BY THE ENTIRE JET 


The quantity of interest is not l^ol* the intensity radiated by the 
whole volume of turbulent flow. We introduce for this purpose the axial 
scale parameter a (as used to scale mixing layer profiles, where the trans- 
verse coordinate is written as ay/x). The parameter a is used to indicate 
the axial extent of the initial region of the jet; thus the total intensity 
contributed by this region (up to the end of the potential core) is expected 
to scale according to 
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(4C-5) 


I 


8 

^max 


a, 


where the maximum value of q across the shear layer. 


FREQUENCY SCALING 


In a truly two-dimensional mixing layer, the relation (^C-5) would 
apply not only to the overall intensity but to the proportional-band inten- 
sity at any frequency, since the radiation at any frequency f would be^ 
associated with a typical shear layer thickness q^^g^/f. However, the jet 
exit diameter d sets an upper limit to the value of 6 for which the shear 
layer can be modelled as two-dimensional. The corresponding frequency 
scale q^nax^^ sets a lower limit to f , if (^C.5) to remain valid. 

For this reason, the intensity scaling prediction (4C-5) should be 
applied to corresponding frequency bands on the non-dimensional frequency 
scale 


St 


fd 
' » 
^max 


(4C-6) 


with a lower limit on set by the end of the jet potential core. Sound 
radiated predominantly from the region downstvecun of the potential core is 
expected to scale differently (uut is not considered here) . 


FLOW PARAMETERS 

In order to predict forward flight effects on jet noise at various 
temperatures on the basis of (4C-5) > we require the variation of q^ax^ with 
the following parameters. 

(a) Shear layer velocity ratio, X 

(b) Shear layer density ratio, Pq/Pj 

(c) Mach number (in case compressibility effects are important). 

Subscripts (i,o) refer to the inner (or high-speed) and outer (or low-speed) 
sides of a two-dimensional shear layer. For the case of a round jet in 
forward flight, Vj is identified with the jet exit velocity Vj, and Vq with 
the forward flight velocity. 
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EMPIRICAL CORRELATIONS FOR (q^a) 


Re table measurements of q over a combined range of A and p /p. are not 
befole^ *■'' '"easurements may provide the missing information 

ll-i 126?"^ constant-density mixing, however, the data of Yule (refs, 

nfn sufficient to show that the simple scaling law q - (Vj - V ) is 

inadequate; an alternative relation is proposed below. ' ° 

the axial scale parameter a, a limited amount of data - not 
rely self-consistent - does exist on the combined effects of a and 
P^/Pj variations, but there is no generally accepted prediction method The 
situation IS summarized in a useful review by Birch and Eggers (ref. 70). 

empirical relations are offered as a provisional fit to 

as a dividing-streamline velocity Vh! 

as a correlating parameter (see following section). ^ 


(a) For the rms velocity fluctuation, the relation 


^max 


(Vi - Vj 3 > 


(i*C-7) 


IS proposed. The actual profile of q across the shear layer is assumed to 
(V. function of the mean-velocity profile variable (Vj - Vj / 

Dredic?! Parametev, a simple relation which roughly 

(velociL H . ® parameters listed above 

ivelocity ratio, density ratio and Mach number) is 


a . (-> y 
Wi -Vdsj 


A somewhat better fit to the published data is given by 


a a 


/_v^\ V3 / Pj V 1/3 

Wi-Vds^ ' 


(4C-8) 


(^C-9) 


puJposes!"^'" 'O'- P-'-'n-y prediction 

to a '"tensity given by (W-5) is much more sensitive 

%ax variations than to o variations. 


CALCULATION OF DIVIDING STREAMLINE 


In a self-preserving two-dimensional 
parallel streams, the dividing streamline 


free shear layer between two 
is defined by the property that 
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there is no mass transport across it, on average. Thus, the dividing 
Streamline separates the mass flows of the two streams; it also locates the 
position of maximum turbulent shear stress, and thus has a physical signi 
ficance which makes the choice of (V^s, Pds) correlation purposes appear 

somewhat less arbitrary. 


The position of the dividing streamline, in terms of the mean velocity 
ratio Vj /V- = ^ds ’ deduced from the mean profiles of density and 

veloci ty i n' the shear layer. Values of for various combinations (X, 

p^/pj), have been calculated by Korst and Chow (ref. 67) using the following 

assumpt i ons . 


(a) For incompressible mixing, the mean profiles of specific volume 
p-l(y) and velocity V(y) are the same shape. 

(b) For compressible gas mixing with y constant across the shear 

layer, the mean profiles of stagnation temperature T(;ot(y) velocity V(y) 

are the same shape. 


(c) The profile shapes in each case are represented by an error- 
function curve. 

Values of Ads obtained in this way were used in equations (4C"7) a^d 
(4C-8) to predict (q„ax>°^ variations relating to the standard case = 
Pi/po=U incompressible flow), and the results are shown in Figures 4C.1 
and 4C.2. Figure 4C.1 shows reasonable agreement between the predicted qmax 
variation (for constant-density mixing) and the measurements of Yule (refs. 
125, 126). In Figure 4C.2 experimental spreading-rate ratios Oq/o (where Oq 
refers to the standard case above) are shown from the results of Sabin (ref. 
37), Baker and Weinstein (ref. 127), Abramovich et al) (ref. 12o) , Johnson 
(ref. 129), and Brown and Roshko (ref. 4l). The scatter is considerable, 
but the predicted a variation is roughly consistent with the observed trends. 


JET NOISE PREDICTION 


Combining (4C~7) and (4C-8) with (4C-5) gives the prediction that 

I - (Vi -Vds)® Vds2 , 

I - Vi^ (1 -Ads)^ ^ds^- ^ 


The frequency parameter Sj defined by (4C-6) may be replaced, 

(4C-7), by 

c - fd 

^ds Vjd - Ads) ■ 


in view of 


(4C-11) 
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Figure 4C.l Variation of turbulence intensity with velocity ratio 
in two-stream mixing (low Mach no., Pq/p: =1). 

Data from Yule (ref. 125, 126). 
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(A = 0, .2, .4, .6) PREDICTION BASED ON (4C-8) 


Figure 4C.2 Spreading rate variation with velocity and density 
ratio, in incompressible two-stream mixing 



EFFECT OF INCREASING JET TEMPERATURE 
AT ZERO VELOCITY RATIO 


The effect on jet noise intensity of raising the jet temperature can be 
predicted from (AC-IO); as noted at the outset, what is being predicted is a 
source alteration effect, and any flow-acoustic interaction effects must be 
allowed for subsequently. We denote the isothermal -jet intensity at the 

same jet velocity by (in the corresponding normalized frequency band). 
Then gives 


I 


e.g. 

for 



for 

Tj/To 


^ ~ ^ds)^ ^ds^^hot 

^ds)^ ^ds^ ^ i sotherma 1 

=2, I/Iq = 1 -5 dB; 

= 10, I/Iq = 3,8 dB.^ 


(4C-12) 


These increases in effective source strength are in fact small compared with 
the reductions in intensity expected from flow acoustic interaction effects. 


EFFECT OF INCREASING VELOCITY RATIO 
FOR GIVEN JET CONDITIONS 


The intensity radiated from a jet in a surrounding stream can be 
re a ted, using (4C-10) to the intensity radiated by a jet of the same exit 
velocity and temperature mixing with fluid at rest. The result, shown in 
Figure AC. 3, is a prediction of the forward flight effect on jet noise, 
purely from the source alteration viewpoint. It should be noted that Figure 
AC.3 IS calculated for jet exit Mach numbers of order 1, but by specifying 
the remaining parameters as A and T^oj ./Tq the effect of Mach number varia- 
tions IS kept small. Moreover, inclusion of flow-acoustic interaction 
effects at 90° to the jet axis [by inserting a factor (T„/Tj)3 on the 
radiated intensity] makes a negligible difference (less than i dB) to the 
predicted intensity reduction factor, over the range of parameters covered 
t>y Figure ^C.3. 


^Caloulat^ from tl^ values of Korst and Chow (ref. 67), using a main 
stream Croaoo number - 0, There is no significant change in these 
f^gure8 tf = 0.2 is used. 
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90° TO JET AXIS 
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Figure 4C.3 Predicted effect of forward flight on jet mixing noise at 90 



SOURCE NON-COMPACTNESS AND CONVECTION EFFECTS 


A reasonable hypothesis for the source convection velocity under various 
conditions of jet density and velocity ratio would be to assume 

Vc “ Vds (4C-13) 

For the standard jet (X = 0, Pg/pj “1), for example, this gives V = 0.62 V| 
using Korst and Chow's value for Ajj. 

The scaling velocities Vgj which appear In the source wavenumber- 
frequency spectrum model may be taken as proportional to q__„, and hence to 
(Vi -Vds). 
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APPENDIX AD 


VELOCITY SAMPLING CALCULATION 
AND LV PROCESSING EXAMPLE 




In reference 1 it is argued that if there is no dependence of particle 
number on particle velocity, that is, there is no uniform particle distri- 
bution m the flow, then the true mean velocity is given by the reciprocal 
of the arithmetic mean of the particle periods (the time taken for the 
particle to cross a fixed number of fringes), multiplied by a proportion- 
ality constant. 


K w particle velocity sample Vj is related to the period sample T; 

oy Vj -K/Tj, then the true mean velocity V is given by 

N 

V = KN/ j; T| . 
i = 1 

This IS equivalent to the harmonic mean of the velocity samples since that 
I s g I ven by 

1 = i y _L 

i/'^ i=i ''i ■ 

It can also be shown that the biasing effect leads to a weighting of the 
ve ocity probability density function by a factor proportional to the 
velocity. Thus, by dividing the measured probability density function by 
the velocity and normalizing, the true velocity probability density function 
can be obtained. The measured velocity p.d.f. is denoted by p„ , in a 
discrete form. The mean value of Vg is given by ^ 

K = I Vg p^ . 


If a new p.d.f. is generated in the manner outlined above, the new p.d.f. 
can be given by 


and the mean value is 


= 


n n 


V = 


1 


y ^ 

^ V 
n n 


. "s Pvg 
^ u 


This last term is unity and the remaining expression 
mean of V, the true velocity. A simple example will 
strate this equivalence. 


represents the harmonic 
now be given to demon- 
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Consider the following samples of velocity; 11, 10, 12, 15t 11, 12, 13, 
13, 12, 13, 1^, 1^, 15, 1^, 1^, 16, 1^, 15, 13, 17* The arithmetic mean is 
found to be 13*^0 and the harmonic mean is 13*18. The discrete probability 
density function is, 


< 

o 

12 

13 

14 

15 

16 

17 

Pv5 -05 .1 

.15 

.2 

.25 

.15 

.05 

.05 

Then following the analys 

is above: 






Vs 10 11 

12 

13 

14 

15 

16 

17 

Ps .0659 .1198 

.16^7 

.2027 

.2353 

.1318 

.0412 

.0388 


Using the modified probability density function, the mean value of Vg 
is given by 


V = I VsPs 
s 


and is found to be V = 13*18. 

The procedure for processing the recorded LV data will not be given. 

The unprocessed data is printed out in the form shown in Figure 4D,1. The 
histogram, taken from the computer print out, is shown in Figure 4D. 2. The 
validated points for axial velocities of 225-^ m/s and 709*3 m/s have a 
considerable effect on the higher order moments. In view of their high 
values, the jet exit velocity being l6l.5 m/s, they may be Justifiably 
regarded as incorrectly validated noise. The velocity data is therefore 
limited to the range ^5*7 ”173*7 m/s. Having bounded the data in this manner, 
it is reprocessed and weighted by the inverse of the velocity. The resulting 
histogram is shown in Figure 4D.2 with the corresponding computer printout 
given in Figure 4D.3* It can be seen that the true velocity distribution 
found in this manner Is very nearly Gaussian, having a skewness of -.0101 and 
kurtosis of 2.7956. (The Gaussian distribution has correspond i ng values of 0 
and 3*) 

If the flow velocity changes very rapidly the particles may lag behind 
the flow. This would lead to low instantaneous velocity reading for rapid 
flow acceleration and a high reading for rapid flow deceleration. Some 
slight particle lag was expected and observed close to the jet exit at the 
highest Mach numbers where the flow has been rapidly accelerated. For a 
particle already accelerated to the local mean flow velocity, the frequency 
of fluctuations in flow velocity it observes will be lower, in its moving 
frame, than that seen by the stationary observer, such as at the measurement 
volume. For the particle size distribution used in the present measurements, 
there is no observable amplitude or phase lag for frequencies in the moving 
frame less than 1 KHz (ref. 58). Since the expected dominant observed fre- 
quencies at the stationary frame were never more than an order of magnitude 
greater than this particle size distribution was expected to have a negligible 
effect on the measurements. Thus, no correction for particle size was applied 
to the data. 
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Figure 4D.1 Unprocessed velocity data. 
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CORRECTED HISTOGRAMS 
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Figure 4D.3 Processed velocity data 



APPENDIX 5 
LIST OF SYMBOLS 


a. ajj 
aj, aj 

by 

B 

Bij 

Cn 

Cp 

Cp 

CA 

d 

dA 

D 

D„ 

m 

f 

fc 

p(v) 

H 


I 

k 

k 

m 


local, ambient speed of sound 

speed of sound in primary jet, tunnel flow 

primary jet, plenum cross-sectional area 

velocity profile constant 

f i 1 ter bandwidth 

volume displacement quadrupole source strength 
effective nozzle coefficient 
specific heat at constant pressure 
free-jet data correction factor for flow effects 
free-jet data correction factor for distance effects 
eddy convective amplification 
primary Jet nozzle diameter 
, area of ray tube cross-section 
Doppler factor 
modified Doppler factor 
frequency 

center frequency of one-third octave filter 
flow factor 

stagnation enthalpy; or normal height or distance of 
aircraft flight path from microphone 

acoust i c i ntens i ty 

turbulent kinetic energy 

wavenumber vector 

exponent of relative velocity; or 

total mass flow per unit area [Hill (ref, hj)] 


386 



LIST OF SYMBOLS (conf d) 


M 

Me 

Ma 

Mj 

Mr 

Mt 

Mto 

n 

OASPL 

P 

P* 

Po 

Pi 

PR 

PT 

P 

P(o)) 

Pr 

Pr 

Pt 

q 

q(x, t) 


mass concentration on jet centerline 

nozzle Mach number; or axial flow Mach number, V^/a; or 
twice the average sum of momentum and pressure forces per 
unit area [Hill (ref. ^7)] 

axial eddy convection Mach number, V^/a^ 

aircraft Mach number, 

jet Mach number, Vj/aj 

reservoir or plenum airflow Mach number 

tunnel or free-jet Mach number, Vj/a^ 

minimum tunnel or free-jet Mach number 

wavenormal unit vector 

overall sound pressure level 

static pressure 

acoustic pressure 

ambient pressure in anechoic room 

free-jet intake static pressure 

reservoir or plenum static pressure 

tunnel or free-jet test section static pressure 

total pressure 

power spectral density of acoustic pressure 
Prandtl number 

reservoir or plenum total pressure 

tunnel or free-jet test section total pressure 

total turbulence intensity 

acoustic source distribution 
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LIST OF SYMBOLS (coafd) 


Q(k, w) 

r 

r 




'■.5 

R 




RN 


s 

S 


m 


•^1 j 

SPL 


t 

t 

T 


cross-power spectral density of source distribution 
radial distance; or (cylindrical) radial coordinate 
ray path position vector 

(cylindrical) radius at which axial mean velocity is 66.3% 
of centerline value 

ha 1 f -ve 1 oc i ty rad i us 

potential core radius 

primary jet nozzle radius 

tunnel equivalent radius 

normalized (cylindrical) radial coordinate, r/rj 
microphone, measurement radius (spherical) 

(spherical) radial coordinate 
source-observer separation at reception time 
radius of duct 

run number in free-jet experimental program 
distance along ray path 
Strouhal number 
modified Strouhal number 

volume acceleration quadrupole source strength 
sound pressure level 
t i me 

pulse travel time, Rp/a^ 

temperature; or flyover noise record length 
free-jet intake temperature 
temperature on jet centerline 
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LIST OF SYMBOLS (confd) 


Tj 

Tr 

Tt 

V 

Vj, V2, V3 

y 

Vi. V2 
Va 

Vb 

Vc. Vci 

Vds 

Ve 

^ei > ^e2 > 
Vet 
Vr 
Vs 

Vx» V^,V^ 

Vx 

Va 

Vj 

VrEL 

Vt 

VtO 


primary jet exit temperature 
reservoir or plenum flow temperature 
tunnel or free-jet test section temperature 
acoustic particle velocity 

axialy radial y azimuthal velocity fluctuation 

mean flow velocity vector 

axial, radial mean velocity 

axial mean velocity on jet centerline; or 
lower stream velocity in 2-D mixing layer 

upper stream velocity in 2-D mixing layer 

eddy convection axial mean velocity 

dividing streamline velocity 

eddy volume 

wavenumber scaling velocity components 
transverse wavenumber scaling velocity, * 
ray propagation speed 

axial mean velocity of fluid at source position 
mean flow velocity components in x, r, directions 
velocity profile, e*g. V^/Vj 
ai rcraf t speed 

primary jet exit axial mean velocity 
relative veloci ty, (Vj - Vj) or 
tunnel or free-jet velocity 
minimum tunnel or free-jet velocity 
axial distance or axial coordinate 


X 



LIST OF SYMBOLS (conf d) 


Xq 

xn 


X 


I 

s 


X 

a 

Of 


Y 


A0e 

At 


e 


0 

9c 

00 


K 

X 

u 


virtual origin of mixing 
potential core length 

axial protrusion of primary nozzle downstream of 
free-jet exit plane 

axial distance between effective source location and 
primary nozzle 

normalized axial coordinate, x/rj 

scaling parameter in a modified Doppler factor 

wavenormal unit vector 

axial wavenumber scaling velocity constant > ®l^a,rel * 

transverse wavenumber scaling velocity constant, 

^a,rel ” ^et 

ratio of specific heats 

vorticity thickness of shear layer, Va,re|/M''x^^*"lmax 

axial derivative of vorticity thickness, spreading rate 

angular resolution in aircraft flyover tests 

noise record length corresponding to AS^ 

turbulence energy dissipation rate; or 
kinematic eddy viscosity (Appendix AB) 

wavenormal angle (except 0^^) relative to downstream jet axis 
cone of silence angle 

polar angle, relative todownstream Jet axis, at emission 
time 

measured polar angle 
radial wavenumber 
velocity ratio, Vj/Vj 
index of refraction 
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LIST OF SYMBOLS (conf d) 


p 

PQ 

PJ 

o 


•I- 

$(k, u) 
X 


0 ) 


eddy vi scosi ty 
fluid mean densi ty 
density in surrounding flow 
primary jet exit density 

jet flow spreading parameter; or spread rate defined by Gortler 

[ref. (69)] such that similarity of mean profiles were obtained 
by plotting against (oy/x) 

a for A = 0 

(cylindrical) azimuthal coordinate 
eikonal function 

source function cross power spectral density 

similarity coordinate for velocity profile 

ray angle relative to downstream jet axis; or polar angle, 
relative to downstream jet axis, at reception time; or 
stream function 

source frequency; or radian frequency 


Subscripts 

a 

e 

' . j 

J 

m 

rel, REL 
R 
s 
T 


jet centerl ine value 
at emission time 

denote Cartesian coordinate directions (i, j=1,2,3) 

primary jet exit value 

measured value; or modified value 

relative to velocity or Vj 

reservoir or plenum value 

refers to source conditions within the primary jet flow 
wind tunnel or free-jet uniform test section value 
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LIST OF SYMBOLS (conf d) 


Overbars 


ambient value; or stationary observer reference frame in 
flyover case; or external flow conditions (Appendices 3A 
and 3B) 


denote non-dimensional i zed quantities 
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